Plasma Analogy and Non-Abelian Statistics for Ising-type Quantum Hall States 



Parsa Bonderson, 1 Victor Gurarie, 2 and Chetan Nayak 1 ' 3 

'Microsoft Station Q, Elings Hall, University of California at Santa Barbara, Santa Barbara CA 93106, USA 
^Department of Physics, CB390, University of Colorado, Boulder CO 80309, USA 
Deparment of Physics, University of California at Santa Barbara, Santa Barbara CA 93106, USA 

(Dated: February 3, 2011) 

We study the non-Abelian statistics of quasiparticles in the Ising-type quantum Hall states which 
are likely candidates to explain the observed Hall conductivity plateaus in the second Landau level, 
most notably the one at filling fraction v = 5/2. We complete the program started in Nucl. Phys. B 
506, 685 (1997) and show that the degenerate four-quasihole and six-quasihole wavefunctions of the 
Moore-Read Pfaffian state are orthogonal with equal constant norms in the basis given by conformal 
blocks in a c = 1 + | conformal field theory. As a consequence, this proves that the non-Abelian 
. statistics of the excitations in this state are given by the explicit analytic continuation of these 

wavefunctions. Our proof is based on a plasma analogy derived from the Coulomb gas construction 
of Ising model correlation functions involving both order and (at most two) disorder operators. 
We show how this computation also determines the non-Abelian statistics of collections of more 
than six quasiholes and give an explicit expression for the corresponding conformal block-derived 
O^l ' wavefunctions for an arbitrary number of quasiholes. Our method also applies to the anti-Pfaffian 

wavefunction and to Bonderson-Slingerland hierarchy states constructed over the Moore-Read and 
anti-Pfaffian states. 
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I. INTRODUCTION 



Non-Abelian braiding statistical— is currently the subject of intense study, partly because the experimental ob- 
servation of a non-Abelian anyon would be a remarkable milestone in fundamental science and partly because of its 
potential application to topologically fault-tolerant quantum information processing 8 - - — . At present, the state which 
is the best candidate to support quasiparticles with non-Abelian braiding statistics is the experimentally-observed 
v = 5/2 fractional quantum Hall stated—. Efforts to observe non-Abelian anyons in this stated— and harness 
them for quantum computatio n 23 ' 30 " — are predicated entirely on two assumptions: (i) The observed state is in the 
same universality class as either the Moore-Read (MR) Pfaffian stated or the anti-Pfaffian stat o 35 ' 36 , an assumption 
which is supported by numerical studies 37-40 . (There is another non-Abelian candidate, the so-called SU(2)2 NAF 
state 4 ^, for this plateau, but it is not supported by numerics.) (ii) Quasiparticle excitations above these ground 
states are non-Abelian anyons. In order for this assumption to hold, it is necessary for there to be a degenerate set 
of n-quasiparticle states and for quasiparticle braiding to transform these states into each other in such a way that 
i-H different braiding transformations do not commute. 

' Moore and Read^i conjectured that the MR Pfaffian state is non-Abelian while Greiter, Wen, and Wilczek— argued 
that it is Abelian. It was subsequently shown by Nayak and Wilczek — and by Read and Rezayi 4 ^ that there is a 
2Lf J-i-fold degenerate set of n quasiparticle states. To show that assumption (ii) is correct, it is further necessary 
to show that these degenerate states are transformed into each other by non-commuting transformations enacted 
by quasiparticle braiding. Several different argument o 43 ' 45 " — strongly support this hypothesis, but a proof has been 
\ missing until now. By "proof," we mean an argument that relies on no unproven assumptions beyond the existence 
of an excitation gap and the existence of a screening phase for particular classical two-dimensional (2D) plasmas at a 
particular temperature and, therefore, is at the same level of rigor as the Berry's phase calculation for quasiparticles 
in the v = 1/M Laughlin states^ 5 -. In this paper, we supply such a proof by mapping matrix elements of the MR 
Pfaffian state to the partition function of a classical multi-component 2D plasma, possibly with magnetic charges. 
Our derivation extends and completes a partial result obtained in Ref. |45L Numerical studies provide very strong 
evidence that the plasmas corresponding to the v — 1/M Laughlin states with M < 70 are in the screening phased. 
Similar numerical evidence confirming that the plasma (described in our paper) corresponding to the v = 1/2 MR 
state is in the screening phase has recently also been obtained 5 -^. 

One approach to the calculation of the braiding statistics of quasiparticles in fractional quantum Hall states is 
based on an idea due to Moore and Read 3 ^. These authors proposed to use the conformal blocks of conformal field 
theorie s 58 ' 59 (CFTs) as trial wavefunctions for fractional quantum Hall effect states. The conformal blocks are the 
holomorphic parts of correlation functions. Unlike correlation functions, conformal blocks are not single-valued. The 
conformal blocks which are used as trial wavefunctions for fractional quantum Hall effect states are single-valued 
in electron coordinates but are not single-valued in the coordinates of the quasiparticles, and it was conjectured 
that the properties of the conformal blocks under analytic continuation of the quasiparticle coordinates define their 
non-Abelian statistics. 
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However, as emphasized by Blok and Wen£°-, the analytic continuation properties of wavefunctions are only part of 
the story. An additional contribution to the statistics is given by the Berry's matri x 55 ' 60 " —. Wavefunction analytic 
continuation only gives the quasiparticle statistics correctly if the conformal blocks, as electron wavefunctions, have 
matrix elements which are independent of the quasiparticle coordinates (when they are well-separated). This includes, 
but is not limited to, the diagonal matrix elements, which are the wavefunctions' norms. When this condition is 
satisfied, the Berry's matrix is trivial, apart from a term which accounts for the Aharonov-Bohm phase due to the 
(charged) quasiparticles' motion in the magnetic field. This is because the wavefunctions are holomorphic in the 
quasihole coordinates, except for the Gaussian factors (which give rise to the resulting Aharonov-Bohm terms). 

The effective field theory of a fractional quantum Hall state is expected to be a Chern-Simons theory. Chern-Simons 
theories are related to conformal field theories 6 ^: the Hilbert space of a Chern-Simons theory with fixed non-dynamical 
charges at points 771 , . . . , r\ n is equal to the vector space of conformal blocks in an associated CFT with primary fields at 
771 , . . . ,r) n . Thus, if the multi-quasiparticle wavefunctions of a fractional quantum Hall state can be identified with the 
conformal blocks of a CFT, it is very natural to conclude that this fractional quantum Hall state is in the universality 
class of the associated Chern-Simons theory. In fact, one can hardly imagine any other possibility. However, this 
identification is only correct if the braiding properties of the multi-quasiparticle wavefunctions are equal to those of 
the Chern-Simons theory. This, in turn, requires the Berry matrix (in the basis given by the conformal blocks) to be 
trivial. 

Thus, the logic may be summarized as follow s 43 ' 52 . Let us suppose that the quasiparticles of some quantum Hall 
state have the special property that when n quasiparticles are present at arbitrary fixed positions 771 , . . . ,r] n , there is a 
g-dimensional space V n of degenerate states of the system. Now let us suppose that ^ a (Vn'i z i) with a = 0, 1, . . . , q — 1 
are the q conformal blocks of a correlation function in a CFT, where Zi, . . . ,Zn are coordinates of the N electrons. 
(We choose the CFT and the operators in the conformal block so that they are single- valued in the Zi, but possibly 
multi- valued in the 77^.) If the v t , Q (r7 M ; z{) form a basis for V n , then we wish to show that the overlap integrals 

f n 

J k=l 

are proportional to 5 a p diagonal and independent of the quasiparticle positions 77^, in the limit where the 77^ are far 
apart from each other. If we can show this, then the braiding properties of the n quasiparticles are determined by the 
analytic continuation properties of the wavefuctions 4' Q (r7 M ;zi). 

There is significant previous literature which addresses this problem by analytic or numerical method o 43 ' 45 ' 48 ' 51 " — . 
In Section [XIII1 we will discuss these previous results and clarify their relation to the result of this paper. 

In this paper, we prove, for the MR Pfaffian state, that the overlap integrals of Eq. ([1]) are diagonal and independent 
of the quasiparticle positions 77^, in the limit in which the 77^ are far apart. Specifically, we show that 

G a ,ff(v^) = C5 a p + O ( e -l"«-^l/ £ ) , (2) 
which allows us to define orthonormal states (77^)) by dividing by the common normalization constant 

te|#a M = G-^in^rj^^-z,). (3) 

We obtain Eq. @ by expressing the desired matrix elements in the form of the partition function of a classical plasma 
and relying on the screening property of a plasma, thereby extending Laughlin's plasma analogy 65 arguments to these 
non-Abelian states. Our derivation completes the program started in Ref. HH, where such a plasma representation 
was used to prove that the diagonal sum of norms in Eq. (fTJ), ^ Q G QjQ is a constant independent of the quasiparticle 
positions (so long as they are well-separated). The methods used there did not, however, allow one to prove that their 
norms are independently constant and equal, nor that off-diagonal matrix elements G a ^ are zero. We accomplish 
this by extending and elaborating on the methods proposed in Ref. l45l One of the important steps in our approach 
is the explicit construction, via the Coulomb gas formalism 6 ^—, of Ising model correlation functions including both 
order and disorder operators, shown in Eq. (|156l) . This equation is one of the significant results of our paper and is 
interesting in its own right in the context of the Ising CFT. 

Although we can directly calculate the Berry's matrix only for the two-, four-, and six-quasiparticle wavefunctions 
in this way, our results determine the braiding properties of arbitrary numbers of quasiparticles. We show that the 
enumeration of multi-quasihole states 43 [which can be done without computing the integrals in Eq. (fT])] allows us to 
compute the braiding statistics of an arbitrary number of quasiparticles, given a mild assumption of locality. This 
derivation uses special properties of the MR Pfaffian state and works in a particular basis (the "qubit basis"), but 
does not need any further assumptions beyond the existence of a gap in the energy spectrum. 

We can also utilize similar locality assumptions in the form of the more refined formalism of anyon models, which 
describes a topological phase with a braided tensor category. For this, the topological structure is specified by the 
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number of topologically-distinct quasiparticle species, their fusion algebra, the F-symbols (which encode associativity 
of fusion), and the i?-symbols (which encode braiding). As we discuss, the F- and i?-symbols can be determined 
merely from the two- and four-quasihole wavefunctions. Thus, the underlying structure of a topological phase allows 
us to bootstrap from the four-quasiparticle case to an arbitrary number of quasiparticles. In contrast to the previous 
derivation in the qubit basis, this derivation economizes on the necessary input, i.e. not requiring six-quasiparticle 
wavefunctions, because it allows (in fact incorporates) changes of basis, in the form of the F-symbol transformations. 

The results of our paper also apply to the anti-Pfaffian wavefunction, constructed as the particle-hole conjugate of 
the MR Pfaffian stat e 35 ' 36 . They similarly apply to the Bonderson-Slingerland (BS) hierarchical states 6 -^ constructed 
over these, which provide candidates for all the (other) observed quantum Hall plateaus in the second Landau level. 
In particular, this includes BS candidate states for v = 12/5, for which there is also some numerical evidence^. 

The methods we develop here should also be generalizable to other quantum Hall states, most importantly to the 
Read-Rezayi (RR) series of parafermion states? 1 . Doing this in practice requires a careful development of the Coulomb 
gas construction for these states, which has not yet been accomplished, and overcoming additional obstacles that do 
not exist for the Ising-type states analyzed in this paper. This will remain the subject of future work. 

This paper is organized as follows. In Section UH we review the derivation of the Berry's matrix for adiabatic 
processes involving degenerate states. In Section Hill we discuss adiabatic transport of quasiparticles in the MR state, 
and describe the problem to be solved. In Section ITVl we discuss Laughlin's plasma analogy arguments. In Section [Vl 
we review the Coulomb gas construction of the Ising CFT, following Ref. 67. In Section IVlIl we reproduce the result 
of Ref. 45 on the sum of the norms of multi-quasiparticle wavefunctions. In Section [Villi we extend this Coulomb gas 
representation to arbitrary matrix elements of the four-quasihole and six-quasihole wavefunctions, thus proving that 
they are orthogonal with equal norms. In Section HXl we show how these results determine the non-Abelian statistics 
for an arbitrary number of quasiparticles. In Section [Xl we use the plasma analogy to show that two wavefunctions 
(with quasiparticles) are orthogonal if they do not have matching types of quasiparticles at the same coordinates. In 
Section KB we use the previous results to determine the statistics of quasiparticles in the anti-Pfaffian state and BS 
states. In Section IXIH we briefly discuss the application of the methods we have developed to other candidate states 
based on other CFTs. Finally, in Section IXIIIl we discuss previous works that have made progress toward establishing 
the braiding statistics of the MR state. In Appendix [XJ we specify the normalization conventions that we use for free 
bosons. In Appendix [Bl we review Mathur's procedure^ for relating products of contour integrals to 2D integrals 
in the Coulomb gas representation of CFTs. This relation plays a crucial role in our analysis. In Appendix [D] we 
use the Coulomb gas representation to compute the (multi- valued) correlation function of two order and two disorder 
operators in the Ising model. In Appendix [D] we discuss the behavior of electric and magnetic operators in the 
plasma phase of a two-component Coulomb gas. In Appendix [El we review and generalize the Debye-Hiickel theory 
for application to the plasmas that arise in this paper. In Appendix IF1 we compute the 2"2 _1 conformal blocks of 
n o fields (where n is even) and an arbitrary number N of ip fields in the Ising model; this gives a preferred basis 
for the q = 2^ _1 degenerate states of n quasiholes in the MR Pfaffian state. The Berry's matrix is trivial in this 
basis and braiding properties are given explicitly by the analytic continuation properties of these wavefunctions. In 
Appendix [Gj we give an incomplete argument that would allow one to compute the braiding statistics for an arbitrary 
number of quasiparticles directly from the wavefunctions with arbitrary numbers of quasiparticles. Although, as we 
show in Section IIX1 this is not necessary, it would nevertheless be a particularly simple and elegant route to deriving 
quasiparticle statistics, if it could be completed. In Appendix [HI we provide two explicit examples demonstrating the 
orthogonality of wavefunctions that do not have matching quasiparticle types at the same positions. 



II. BERRY'S MATRIX 



In this section, we review the derivation of Berry's matrij* 6 ^— for an adiabatic process when there are energy 
degeneracies. We consider the Hamiltonian H (R±(t), . . . ,R n (t)), which depends on a set of parameters R/j,(t) that 
are varied in time t. All states evolve according to the Schrodinger equation 

ih±mt))=H(t)m))- (4) 

One can define orthonormal energy eigenstates \a(R±, . . . ,R n )) f° r the Hamiltonian at particular values of the pa- 
rameters R^, such that 

H (RJ \a(R^)) = E a (Rfj,) \a{R^)) (5) 

and {a{Ra)\ (3(R^)) = 5 a p. When the parameters R^ are varied with f, we will leave the R^ dependence of quantities 
implicit, e.g. writing H (t) and \a(t)). We consider a Hamiltonian such that the Hilbert space splits into subspaces 
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of degenerate energies H(t) — Y^Ett) U-E(t)- We now focus on one of these subspaces HE (t) ( e -g- the subspace of 
ground-states), and assume that the energy gap between it and the other subspaces does not close during the adiabatic 
process. The adiabatic theorem tells us that if we start at t = with a basis state \ip a (Q)) = l a (0))) G %_e (o)> then 
the time evolved state \ip a (t)) will be in the T-LE {t) subspace, and can thus be written in the form 



|Va*)) = e-i^o(0*V (t)|a(t)) 



(6) 



where Uq is the Berry's matrix, which is a generalization of Berry's phase. It is a unitary transformation in the Eq 
subspace, i.e. Uo(t) : %£ (t) — > TJ-E (t), such that Uq(0) — 1, and the dynamical phase exp — j? J* E (t')dt' has been 
separated from the Berry's matrix term. Since it is a matrix, the Berry's matrix can potentially exhibit non-Abelian 
behavior. Taking the time-derivative of Eq. [6] and taking an inner product with another time-evolved state in HE (t)i 
we have: 



ih{Mt)\ | \Mt)) = E (t) (a(t)\ P{t)) +ih{a(t)\ Uo\t)^P- |0(t)) +ih(a{t)\ | \/3(t)) 
Re- writing the left-hand-side by using Eq. (|4]), one finds 

ih^ a (t)\j t \Mt)) = (Mt)\H(t)\Mt))=E (t) (a(t)\ 0(t)). 



Combining Eqs. ([7]) and ([8]), we obtain 



{a{t)lU -i {t) ^p. m) 



Ht)\ Jt \m) 



(9) 



Solving this expression for Uq, one finds 



U (t) = Pexp 



i / A{t')dt' 



t + y^i n I dt n dtn-x... dt 1 A{t 1 )...A{t n ) 
Jo Jo Jo 



(10) 



where V stands for path-ordering (putting operators to the right of those with smaller t and to the left of those with 
larger t), and we have defined the Berry's connection for the HeqU) subspace 



A^(t) ee i(a(R l7 ... 7 R n )\^ r \f3(R 1 ,...,R n )) 



(11) 
(12) 



Defined this way, A is Hermitian. 

The term Uo(t) only has a gauge- invariant meaning if the Hilbert space is the same as the original one. For 
this, one must make a closed circuit in configuration space. Let us consider an adiabatic process with t running 
from to tf (where tf is large enough compared to the inverse of the energy gap that the process is adiabatic), 
where %E {t f ) — ^e (o) an d the path in configuration space is a closed loop, which includes processes that exchange 
identical (quasi-)particles. Even though HE (t f ) = He (o)i it is possible to have \a(tf)) ^ |a(0)), e.g. if we have 
defined \a (i?i, . . . , R n )} which is multi-valued as a function of the R^. However, they must be related through a 



transformation B : % 



Eq(0) 



Ke ( ), defined by 



B a ,p = (a(fl)\j3(t f )) 



(13) 



so that \a(tf)} = B |a(0)). For such an adiabatic process, we can now write the time-evolved state at t ~ tf in terms 
of operators acting on the initial state \ip a (0)) = |a(0)) £ T~Le {o) 



\ip a (t f )) =e-*io' / E (t')dt' p 



exp 



A(t')dt' 



(14) 
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and thus, it may be applied to an arbitrary initial state |\P(0)) = J2a c <* IV'a(Q)) = J2a c a l a (0)) m tne subspace We (o) 



exp 



A(t')dt' 



B\9(p)). 



(15) 



If we never consider states outside the %E {t) subspace, we can obviously ignore the common dynamical phase. Thus, 
we see that the evolution of the initial state in the He (o) subspace under an adiabatic process is (apart from the 
common dynamical phase) composed of the Berry's matrix and the wavefunction transformation B. 



III. QUASIHOLE WAVEFUNCTIONS AND NON-ABELIAN STATISTICS 

In this paper, we will be discussing a set of wave functions and their braiding properties, i.e. the evolution under 
adiabatic exchange of quasiparticles in 2D systems. We will make little reference to the Hamiltonian of the system, 
other than to assume that the Hamiltonian has a gap above its ground state(s). The wavefunctions which we discuss 
can be regarded as trial wavefunctions for the Hamiltonian of electrons in a magnetic field interacting through the 
Coulomb interaction. Alternatively, they can be viewed as exact eigenstates of electrons in the lowest Landau level 
at filling fraction v = l/M interacting through a special model Hamiltonian with three-body interactions, 

H = H^ 1 . (16) 

For the case of bosons at v = 1, the Hamiltonian has the form 

N 

Hl = X J2 <P(* **)■ (17) 

i<j<k 

where A > 0. For the case of fermions at v — 1/2, Fermi statistics dictates a more complicated for m 38 ' 42 : 

JV 

Hi = X ]T S ijk {dtd*}5 2 (zi- Zj )6 2 (zj - z k ). (18) 

i<j<k 

where Sijk is a symmetrizer. Our focus in this paper will be wavefunctions with an even number n of quasiholes. For 
the model Hamiltonians in Eq. (|16|) . the n quasihole wavefunctions which we will discuss are zero-energy eigenstates. 
(This is typical for such ultra-local Hamiltonians; quasiparticles cost finite energy, so there is a finite energy cost for 
a quasiparticle-quasihole pair.) As we will see, when we fix the positions r)i,...,r) n of these quasiholes, we will still 
have a 2t _1 -fold degenerate space of states spanned by ^> a , a = 0, 1, ... , 2? _1 — 1. For the sake of precision, let us 
momentarily assume that the system is on a sphere of fixed area and that the number of electrons is fixed (and that 
the magnetic field is tuned to accommodate n quasiparticles). Then the only assumption that we will need about the 
spectrum of the Hamiltonian of Eq. (|16p is that all other states with quasiholes at 771 , . . ., r) n will be separated from 
span^a) by a finite energy gap. 

When we consider states with quasiholes, we will need to augment this Hamiltonian with a potential which pins 
the quasiholes at fixed positions: 

H = H 3 + -ffpinning- (19) 

This is necessary to guarantee that there is a gap in the multi-quasihole case; otherwise, it would cost no energy to 
move the quasiholes to other positions. An elegant choice of pinning potentials is constructed in Ref. [Hi]. However, the 
Berry's matrix is computed solely from a set of wavefunctions, with no explicit reference to the Hamiltonian, apart 
from the assumption that it provides a gap. Thus, the pinning potential, though important as a matter of principle, 
is not, as a practical matter, important in its details for our calculation. 

The MR Pfaffian ground state wavefunction for an even number N of particles is given by 3 - 4 -: 

N N 

Hi(z 1 ,..., ZN )=Pi(—Lj\ TT( Zi _ z .)^ e - 4 £ W2 . (20) 

M is a positive integer, taking odd values if the particles are bosons (which may occur, e.g. for neutral bosons 
in a rapidly rotating trap 72 i 73 ) and even values if they are fermions (e.g. electrons in the quantum Hall effect). 
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Throughout most of the paper, we set the magnetic length Ib — y^hc/eB to 1, as we have done in Eq. (f20|) . and will 
only reconstitute it when it provides necessary clarification. The symbol Pf stands for Pfaffian: 

N/2 

Pf E jyjj E S S n ( CT ) II A °(2k-l)M2k) (21) 

cr£Sjv k— 1 

where A is an antisymmetric N x N matrix (where N is even) . The square of the Pfaffian of an antisymmetric matrix 
is equivalent to the determinant, i.e. [Pf (Aj.j)] = det(Aj.j). This wavefunction has the same form as the BCS 
wavefunction in real-spac o 34 i 42 multiplied by a Laughlin-Jastrow factor. 

The wavefunction in Eq. ([20]) is the unique exact ground state of the Hamiltonian in Eq. (fl6| . The M = 2 case is 
an approximate ground state for electrons with Coulomb interactions at v — 5/2 (assuming that the lowest Landau 
level of both spins is filled and the wavefunction in Eq. (|20p is transposed from the lowest Landau level to the second 
Landau level)2Zr— . The M = 1 case is an approximate ground state for neutral ultra-cold bosons in a rotating 
trapZ^iS. 

The wavefunction in Eq. (|2T)|) can be written as a conformal block in a CFT, as was first proposed in Ref. [34. The 
relevant CFT is a (restricted) product of two theories, one at central charge c = 1/2 describing the Pfaffian part of 
the wavefunction, and the other at c = 1 describing the Jastrow factor Oi<j( z i — z o) °f the- wavefunction, as well 
as the Gaussian factor. Specifically, one writes 



( v 2 - v {z N ) e -i 1I ^ m id 2 z V (z) 



N N 

= Pf(^) xn^/." 1 ^ 1 ' (22) 

Here tp represents the holomorphic free Majorana fermion (the operator with conformal dimension — 1/2) of the 
c = 1/2 Ising CFT, and if is the free boson of a U(l) CFT. Various conventions can be used to describe the free 
boson. We adopt the one presented in Appendix fXl with Eq. (IA13|) and g = 1/4. 

For future reference, let us note that the c = 1 correlator is charge neutral, that is, it is invariant under the change 
ip —> <y2+const. Indeed, under such a change, the exponential factor acquires a term N^M/2 — Aj (2iry/2M), where 
A is the total area. However, M — A/2irN is the inverse filling fraction of the quantum Hall state, since A/2ir is the 
total number of available states in a Landau level which we fill with N particles, and so NyjM/2 - A/(2nV2M) = 0. 

An excited state wavefunction depends on the positions Zi of the electrons as well as the positions of the 
quasiparticles. It is important to recognize that the quasiparticles' coordinates are simply parameters of the electrons' 
wavefunction (and underlying Hamiltonian), not to be treated on the same footing as the electrons' coordinates. These 
wavefunctions were constructed as eigenstates of Eq. (fi"6j) in Refs. 14311441 . Given that the ground state can be expressed 
as a conformal block in the c = i + 1 CFT, it is natural to try to construct wavefunctions with n (fundamental) 
quasiholes in the same CFT. The natural guess^ is that they are given by: 

*q (ryi,...,?7„;zi,...,ZAr) = (a(rjx) . . . <r(r) n ) <i/)(zi) . . . V>(zjv)) q 

75S7^) . . . e l ^7m^) e ^^) . . . e'Vf^^lf,-'^ / d ^ *>M \ . (23) 

Here a are the holomorphic spin operators of the Ising CFT, with conformal dimension h a = 1/16. The bosonic part 
of the correlation function is chosen in such a way that the wavefunction is a polynomial function of the Zj. 

Notice the index a in Eq. (1231) . The holomorphic spin operators of the Ising CFT have many conformal blocks, 
which we label by the index a. In fact, it is well-known that the total number of conformal blocks is 22" -1 , thus 

a = 0,l,...,2*- 1 -l. (24) 

The wavefunctions f a represent the set of degenerate wavefunctions at fixed positions of the quasiholes and form the 
basis for their non-Abelian statistics. 

To find the wavefunctions of Eq. (|2"3"]) explicitly, we need to evaluate the appropriate conformal blocks of the CFT. 
For n = 2, there is only a single conformal block for Eq. (|23[) ; evaluating it for N even, we obtain the two-quasihole 
wavefunction: 

TV 

/III, — " I I 1 1. , — " i — i I ; i \ 

^{rji,r]2;z 1 ,...,z N ) = (771 -772)417-8 Pf 



V Zi — Zj J ,^ , 

(25) 



*<3 
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This wavefunction is, indeed, a zero-energy eigenstate of the Hamiltonian in Eq. (|T6|) (see Refs.[43jj44jfor details). Since 
there is only a single generator for the two-particle braid group, a counterclockwise exchange of the two particles, non- 
Abelian effects cannot be seen - they require at least two different braids which do not commute with each other. The 
effect of braiding can, therefore, only be a phase which is acquired by the wavefunction in Eq. (I25[) . This wavefunction 
is single-valued in electron coordinates, as it must be, but is multi-valued in the quasihole coordinates. Taking the 
analytic continuation of this wavefunction at face value, we would conclude that the effect of a counterclockwise 
exchange of two quasiholes in this state is a phase exp [z7r (^jjy — g)] ■ However, this conclusion is premature, because 
we must also take into account the Berry's matrix (which, in this case, is simply a phase). 

Before discussing the Berry's matrix, let us consider the four-quasihole wavefunctions and, briefly, the general n 
quasihole wavefunctions (with n even). In the four-quasihole case, we are faced with evaluating Eq. (f2"3")l for n = 4. 
This calculation is more difficult, but was accomplished in Ref. |43|. For N even, it results in the following two 
wavefunctions 



-A- 1 I (?7n??24) 4 / / s ^OT £ l^cl 2 

*0(Th.,m,ria,W,ZU'--,ZN) = [[V™ 8 ) niU ^' ■ (fr(l 3 )(24) + Vl ~ X fr(14)(23)) e - =1 , 

„o y/l + y/l-x 

4 1 1 4 2 

*i(r)i,V2,V3,r h ;zi,---,ZN) = f[ V$~ k - ^- 13??2 ^_ (tt (13)(24) - %/l~^ ( i 4 )(23)) e OT -=^ (26) 



where the so-called anharmonic ratio x, well-known in CFT, is given by 

_ 7?i2?734 _ (m - m) (m - m) 



(27) 



Here, wg have introduced ttiG notation tj^u — Tj^ — 77^, and the shorthand ^ ( a b){cd) f° r 

N N 

* (ab)icd) = Pf ( iVa - ^ - ^ ~ - Zj) + {l " j) ) nfe-^^e-S" 12 (28) 

V z% Zj J . . 

The wavefunctions v I'(i3)(24) and 4 r (i4)(23) are zero-energy eigenstates of Eq. (IT5|) and they form a basis of the two- 
dimensional space of states with four quasiholes at fixed positions^. The state x I , (i2)(34) is not linearly-independent 
of these two because of the identity^: 

*(12)(34) - *(13)(24) = (1 - X) ( V I'(12)(34) ~ *(14)(23)) • (29) 

Even though ^m)(24) and v I / (i4)(23) form a basis of the four quasihole Hilbert space, they do not provide an orthonormal 
basis. In this paper, we demonstrate that the linear combinations VPo and defined in Eq. (|26|) are, in fact, an 
orthogonal basis. Moreover, we show that V&o and \&i have the same norms (though we do not compute the precise 
value of their overall normalization constant), and thus can provide an orthonormal basis by dividing by a common 
normalization constant. 

It has been argued since Ref. that using the wavefunctions in Eq. (|26|) allows us to read off the non-Abelian 
statistics of the quasiparticles in a straightforward manner. Indeed, if the quasihole at r\\ is exchanged with the 
quasihole at 772 in a counterclockwise fashion (or, equivalently, if the quasiholes at 773 and 774 undergo counterclockwise 
exchange), a straightforward analytic continuation of the wavefunctions leads to the transformation rules: 

771 r\ 2 or 773 774 : * ^ e^^^^^o, *i H> e"^™"-;?)^. (30) 

To see this, we note that 1 — x = 77147723/77137724 ^ and ^(i3)(24) x f'(i4)(23) under this exchange. We see that 
the phase exp [iir \ — \ )] acquired by ^0 is the same as that acquired from counterclockwise exchange of the two 
quasiholes in the n — 2, N even case. 

On the other hand, if the quasiparticles at 772 and 773 undergo counterclockwise exchange (or if the ones at 771 and 
774 are exchanged), then we get 

r) 2 <=> 773 or 771 ±=i 774 : \P e "Um + sJ u — ; ^ ^ ^Wa/ + S J ' . (31) 

v 2 V 2 

Finally, if the quasiparticles at 771 and 773 undergo counterclockwise exchange (or if the ones at 772 and 774 are exchanged), 
then we get 

771 ±5 773 or 772 ±3 774 : *o >-> e^A+i) 2°+^- , *i h> e™( A+l) + ^ . (32) 

v 2 v 2 
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These exchange transformations are more difficult to show, but can be checked using algebraic manipulations as in 
Ref. |H. These three exchange operations Eqs. (I3TJ1) . (f3"Tj) . and (f3"2"j) constitute the building blocks of the non-Abelian 
statistics of states with four quasiholes. 

The explicit form of the conformal block wavefunctions for n > 4 was not previously calculated. In Appendix [F] 
we show that they have the following form: 



/n/2 

II V2i-l,2j-l %»,2j 
i<j 



(Pi,P 2 > 



\ 




r ■ p 



n/2 

n 

k<l 



-1/2 



> k -r,\/2 
J k,l 



K-r,\/2 



(1+^,3+^, ...)(2-r 1 ,4-r 2 , 



(33) 



The indices take the values p, = 0, 1, with the constraint that Y^i=iP i * a even i so there are 22 1 such wavefunctions. 
(If we were to consider the case where the number of electrons N was odd instead of even, then we would instead 
require Y]™- 2 Pi to be odd.) Here, Xki = !| 2fc ~ 1,2 ' r,2l ~ 1 ' 2k and "Fn+r 3+r )(2-r 4-r 1 i s a generalization of the 

r /2fc-l,2I-l ^2kJ2.l 1» 1 2> M 13 2' i 

notation in Eq. ([28]) which was introduced in Ref. |43| and is explained in Appendix [Fj As discussed there, the 
wavefunctions of Eq. (1331) form a basis of the 2"? _1 -dimensional space of zero-energy n quasihole eigenstates of the 
Hamiltonian in Eq. (fT6]). For the special case n = 4, Eq. (|33]) is identical to Eq. (|26[) . The analytic continuation 



properties of wavefunctions with an arbitrary number of quasiholes can be read off from Eq. 

However, calculating the explicit analytic continuation of the wavefunctions is not, in principle, sufficient to establish 
the statistics of quasiholes. One also needs to calculate the Berry's connection. It is defined as 



■4a,/?(*) 



E 

i"=i 



dt 



^ dt 



A 'If - 

where the overlap matrix is defined by: 



N 



if 

k=l 
N 

IF 

fe=l 



Zk 



Zk 



*q d 

* a d 

Gi /2 a 9% 



G 



1/2 



^13 



N 

I II ^ Zk *«(^! Zi^piVnl Zi). 
J k=l 



(34) 



(35) 



(36) 



(37) 



We have allowed for the wavefunctions in Eqs. (|35|) and ()36|) to be un- normalized, since we will not determine the overall 
normalization constant of the wavefunctions we work with in this paper. When the quasiparticles are adiabatically 
transported along the coordinate paths Tj^it), forming a closed circuit in parameter space as t goes from to it/, an 
arbitrary state 'J in the (2 "2 _1 -dimensional) degenerate ground-state space is transformed under the following unitary 
evolution, combining the explicit transformation of the wavefunctions resulting from analytic continuation with the 
Berry's matrix transformation resulting from the Berry's connection (see Section [IT] for more details) 



l*M*/)]>=Pexp 



A(t)dt 



BMV»(0)]) 



(38) 



where V stands for path-ordering, and B is the unitary transformation describing the analytic continuation of or- 
thonormal states 



l*«M*/)]> 



9-1 

E 

/3=0 



B/i.al^MO)]) 



(39) 



(We have dropped the overall dynamical phase, since it is the same for all states in the ground-state space.) For 
example, the analytic continuation matrices corresponding to the exchanges in Eqs. (1301) . (1311) . and (|32[) (assuming 
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the wavefunctions have equal norms) are, respectively, given by 



r ( AM 



1 

i 



V2 



V2 



1 -1 
1 1 



(40) 



We will show that the wavefunctions in Eq. (121))) are orthogonal for large separations — rj v \ —> oo, such that 

G a ,f, = C5 afi + 0(e-b»-n*\/*), (41) 

where C and I are ^-independent constants. This implies that the Berry's connection is zero, up to terms that give 
the Abelian Aharonov-Bohm phase, as may be seen from the following calculation: 



N 



IF 

k=l 



Zk 



d 







GlildVv KG 1 /}, 



2 */3 



. d 

dip 
■ % x 

l 0„R 

8M ap 



J £Ji d ^ \gUZJ g 
r N 



r l/2 r l/2 



1/2 

13,13 



8M 



0(e 



-\v„-v»\/e\ 



(42) 



We have integrated by parts to go from the first line to the second. Similarly, we have: 

N 



J fe=l u oj 



/3 



__ - r -l/2 r -l/2 
= "'IS** + ^ 



V 8M 



AT 



2fc 



fe=l 



-|«7m-^I/^ 



(43) 



In Eqs. (|42|) and (|43|) - we have used the fact that the dependence of $ a on r)^ and on r/^ comes only through 
the Gaussian factor exp(— giy |?/ M | 2 ) (?/ M and fj^ are considered independent of each other for these purposes). The 
resulting Berry's connection is diagonal in the space of wavefunctions, giving rise to the Berry's matrix 



"Pexp 



dtA(t) 



exp 



— E 



8M 



/i=l 



tf dt fx _ „ *k 
at I dt ^ * 



dt 



i+ o( e - | ^-"" |/£ ), 



(44) 



which is the same for all wavefunctions in the degenerate subspace. When the quasiparticle coordinates are taken 
around a closed loop (or the exchange paths of identical quasiparticles form a closed loop) , this term is equal to the 



phase exp I 



J 2M/ 



which is proportional to the total enclosed area A encircled by the quasiparticles in the counter- 



clockwise sense (area encircled in the clockwise sense contributes negatively to A). This is unlike particle braiding 
statistics which depends only on the enclosed particles and not on the area. By reconstituting the the magnetic length 
£b = y/Hc/eB (which we set equal to 1) in this expression, we see that this phase is simply the Aharonov-Bohm phase 



exp 



1 



E 



8M& 



,, , _ d% 



exp 



A 



'2Ml% 



exp 



e BA 



2M he 



*7777^- = ex P *T7 («) 



he J 



acquired by a charge q = e/2M particle encircling a total flux $ = — BA due to the background magnetic field 
B = —Bz. (We use the convention where Zj = Xj +iyj, which corresponds to holomorphic wavefunctions for electrons 

of charge — e in a background magnetic field B = —Bz.) This reconfirms the interpretation of the given wavefunctions 
as corresponding to charge e/2M quasiholes. As long as Eq. (|4ip is fulfilled, however, no other contributions arise in 
the Berry's matrix. In particular, it docs not affect the non- Abelian statistics, which comes from the explicit analytic 
continuation of quasiparticle coordinates in the wavefunctions. 

There are other length scales one should be aware of when considering non- Abelian quasiparticles. In general, 
topologically non-trivial excitations can tunnel between non-Abelian quasiparticles, which has the effect of splitting 
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the degeneracy of their states^. Such tunneling is exponentially suppressed with separation distance, and thus 
introduces correlation length scales associated with the tunneling of topological excitations and determined by the 
(non- universal) microscopic physics of the system. As long as the quasiparticles are farther apart than these correlation 
lengths, the topological degeneracies are preserved (up to exponentially suppressed corrections), but otherwise the 
notion of the non-Abelian state space and braiding statistics transformations upon it breaks down. For a non-Abelian 
quasiparticles in Ising-type topologically ordered systems, the relevant correlation length is which corresponds 
to tunneling of the tp excitation, i.e. a Majorana fermion. For p-wave superconductors, £^ is identified as the 
superconducting coherence length 7 - 5 -. For the MR state, corresponds to tunneling of the neutral fermion (ifjQ in the 
notation of Section [IX Bp . Numerical studie d 53 ' 76 provide an estimate of £0 « 0.8£b — 2.3£b for the v = 5/2 MR state. 

The wavefunctions in Eq. (|2l)|) were derived as correlators of some CFT. There is no reason a priori to expect that 
they will form an orthogonal basis obeying Eq. (|4"Tj) with respect to the inner product of nonrelativistic electrons in 
a magnetic field. It is the goal of this paper to show that this is indeed so. 



IV. LAUGHLIN'S PLASMA ARGUMENT 



We proceed by first recalling an argument due to Laughlin 7 - 7 - which he used to deduce the normalization of the 
Laughlin wavefunction with N electrons and n quasiholes in the v = l/M quantum Hall effect. Such a wavefunction 
has the form 



,Z N ) 



n N 



N 



i— 1 i<j 



M. 



E l^l 2 



E Nr 

i = l 



(46) 



Note that the prefactor Yiu<v( r ln ~ depends only on the quasihole coordinates r/^ and is independent of the 

electron coordinates Zi&. Therefore, it can be regarded as part of the normalization of the wavefunction. By including 
it explicitly in the definition of the wavefunction, we are anticipating that it will result in a norm of the wavefunction 
that is independent of the quasihole positions. Laughlin proved that 



N 



nd 2 z k * 1 # 1 =C\ 



O e 



-\V»-Vu\/ti 



(47) 



fc=i 



where C\ and £\ are constants independent of r/^ . (We use the subscript 1 here to indicate quantities that correspond 
to the one-component plasma and to differentiate them from similar quantities occurring elsewhere in the paper.) 
The proof proceeds as follows. One observes that the normalization integral Eq. (I47|) can be rewritten as 



. N , N 

J k=l J fe=l 



2 Zfee -* 1 /T = e -F 1 /T 



(48) 



JY 



Yj t?2 lo § \vm - vA - Yl 1? lo § 1^ - «ii ~Yq 2 1o § 



fj,= l »=1 



N 



(JL=1 2—1 



(49) 



where T — Q 2 /2M. We note that the 2D Coulomb interaction between two charges q\ and qi separated by a distance 
R is — qi q2 log R. Thus, $1 can be interpreted as the 2D Coulomb-interaction potential energy for N charge Q 
particles at and n charge Q/M particles at 77^, together with a uniform neutralizing background of charge density 
p = — 2 ^M£2 [which is the uniformly negatively charged disk, represented by the Gaussian terms in Eq. (|4"6")) ] 125 . 

Consequently, F\ can be interpreted as the free energy of a classical 2D one-component plasma at temperature T of 
N charge Q particles in the presence of n additional test particles of charge Q/M at the fixed positions ?y M and a 
uniform neutralizing background. Clearly, one can ascribe different charge values Q to the plasma particles, as long 
as one similarly alters the test charges and temperature in a compensating manner. One convenient choice is to take 
T = g and Q = \J2Mg. (Another typical choice is T = M /2 and Q = M, which gives the test particles unit charge.) 
In any case, the coupling constant T = Q 2 /T = 2M remains invariant under such redefinitions, and it is known from 
Monte Carlo simulations 56 that the freezing point of such a classical 2D plasma is at T Cl w 140 (i.e. T Cl w Q 2 /140). 
Hence, the plasma is a screening fluid for M < 70, whereas it freezes into a crystal for M > 70. 
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It is important to distinguish this M Cl w 70 transition point between the fluid and crystal phases of the analogous 
2D one-component plasma from the M c ss 9 transition point between the quantum Hall fluid and Wigner-crystal 
phases of the physical electron systems^. The later determines the physical range where quantum Hall states exist, 
while the former indicates that the plasma analogy argument indeed applies to the Laughlin wavefunctions for all the 
physically relevant filling fractions. 

When the plasma screens, the free energy F% in Eq. (|48|) cannot depend on the positions rj^ of the Q/M test charges, 
so long as 1^ — r\ v \ 3> l\, where l\ is the Debye length of this plasma, since they are screened by the elementary 
charges. The Debye length can be estimated using Debye-Huckel theory (see Appendix [Ej to be l\ = Ib/V%, where 
£b = y/Hc/ 'eB is the magnetic length (which we have set to 1) of the quantum Hall system. Thus, the overlap integral 
is indeed a constant, as long as the test charges are sufficiently far away from each other. 

It follows that the Berry's connection for adiabatically transporting Laughlin quasiholes using the wavefunction as 
normalized in Eq. (I46[) is given by 



A 7h ' 



j fc=i 1 



* 3 / * 



«£(M-/n**f)-«/n* 



Z k 



9% 



1*1 



(50) 



and 



r N t 

a* = i n w 

j i — i ii ■ ■ 



W d ( * 



k=l 



d% VII* 
4M v ; 



(51) 



This gives a Berry's phase of 



V exp 



i I dtA(t) 
o 



= exp 



exp 



1 

AM 



E 



dt 



dt 



dt 



0(e 



-\v»-v*\/ei) 



M he 1 



\Vn-Vu\/ti\ 



(52) 



where A is the area encircled by the quasiholes in the counter-clockwise sense. This contributes only the Aharonov- 
Bohm phase exp(iq&/hc) acquired by charge q = e/M encircling an area A containing flux $ = —BA from the 

background magnetic field B = —Bz. The remaining contribution to the unitary evolution resulting from adiabatic 
transport of the quasiparticles comes from explicit analytic continuation of the wavefunction, which is thus the braiding 
statistics of the quasiparticles. This proves that the Laughlin quasiholes are anyons that accumulate a statistical phase 
9 = tt/M as the positions of two of them are exchanged in a counterclockwise fashion, as can be explicitly seen from 
analytic continuation of the term (ry^ — rj v )^i in the wavefunction of Eq. (J46j) . 



V. THE COULOMB GAS CONSTRUCTION 



A. Intuitive Approach 



In the previous section, we saw that, although we could not explicitly evaluate the norm of the Laughlin wavefunc- 
tion, we could make a strong statement about its dependence on quasihole coordinate by appealing to the screening 
property of a Coulomb plasma. We would now like to construct such an argument to prove Eq. (|41|) . but we must 
first note that, taken at face value, the overlap integrals of *S> a , defined in Eq. (|37|) . have little to do with the partition 
function of a plasma. Indeed, the plasma argument seems to be custom tailored for wavefunctions which can be 
written as products of differences, such as Eq. (|4l)|) . The MR Pfaffian ground state wavefunction Eq. (f2"U|) . or the 
wavefunctions with quasihole excitations, such as Eqs. (|23|) . (|25|) , and (f26|). are in fact sums of products and cannot be 
written as exponentials of logarithms. Nevertheless, there exists an approach, called the Coulomb gas construction, 
which allows one to represent conformal blocks in terms of a plasma. Let us review this approach, in its particular 
application for the c = 1/2 Ising CFT of interest here. 
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In the next subsection, we will follow the logic and notation of Feigin and Fuchs 7 - 9 -, Dotsenko and Fatee v 66 ' 80 , 
Felder— , and Mathur— , which is essentially an algebraic approach to the Coulomb gas construction of the minimal 
models 58 . However, there is another approach to the Coulomb gas construction which is more intuitive; we briefly 
describe it here (see Refs. l8ll - [83l) . 

The basic question which we answer in this subsection is: Why does the Ising model, which does not have a conserved 
U(l) charge, have anything at all to do with a gas of electric and magnetic charges interacting logarithmically? One 
way of answering this question lies in the following steps. 

1. Write the Ising model partition function on the honeycomb lattice (we choose this lattice for convenience) in 
the form 

Z= e ~ PH = n cos hf3J[l +tanh(/3J)er i o- i ] (53) 

{o-i=±l} {cn=±l} i,j 

where the Ising Hamiltonian has the form H = —J^l<i j> a i (J j- 

2. In any term in the expansion of the product in Eq. (|53[) . each lattice bond can either receive a 1 or an xo~iO~j, 
where x = tanh/3J. Notice that any term in the expansion vanishes upon summation over {er^ = ±1} unless 
every spin <Ji appears either zero times or twice. Consequently, the bonds which receive an xctiOj form closed 
loops, and the partition function takes the form: 

Z(x) = ^ b{a) (54) 
M 

where {a} is a configuration of loops on the honeycomb lattice and 6(a) is the total length of all of the loops in 
the configuration {a}. The critical point of the Ising model occurs at x — 1/V3. 

3. Observe that this partition function can be obtained from the following local rules for the Boltzmann weights. 
When a loop turns left, it acquires a factor xe lx and when it turns right, it acquires a factor xe~ lx . Since the 
number of left turns minus the number of right turns is ±6 for any closed loop, every such loop receives a factor 
2x b cos 6x after summing over both orientations of the loop, where b is the number of bonds in the loop. We 
obtain the Ising model partition function in Eq. (1541) provided 2 cos 6% = 1. 

4. Write the critical partition function (as defined by these local weights) as a height model on the honeycomb 
lattice. A height model is a model of a fluctuating interface which is specified by its local height z — <f)(x, y). The 
loops are interpreted as domain walls between regions with different heights (the heights live on the plaquettes 
and the domain walls on the links). In the continuum limit, the energetic penalty for domain walls between 
different heights becomes a gradient energy (V</>) 2 so that the partition function for the interface can be viewed 
as a quantum field theory for a scalar field </>. 

5. Write the height model as a free bosonic field </> (i.e. as a Coulomb gas) with stiffness g = 1 — ^ = 4/3 together 
with a coupling (1 — g) = —1/3 to the scalar curvature R. (Note that we had to take — 7r < 6% < in order to 
obtain the critical point of the Ising model; taking < 6% < tt would give us the low-temperature fixed point.) 
The stiffness gives the correct energy penalty for a domain wall. A background charge is necessary to correctly 
describe the coupling of the bosonic field to the curvature, because the number of left turns minus the number 
of right turns will be different from ±6 around a point of non-zero curvature. For the honeycomb lattice on the 
plane, this reduces to a background charge at infinity 2(1 — g) = —2/3. 

6. In terms of the bosonic field, </>, the effective action of the height model has a marginal operator uie~ 21 ^ which 
enforces the fact that the heights take values that are integral multiples of ir (which would otherwise be lost 
in the passage to the continuum limit). In fact, the term which enforces the integrality of the heights is more 
complicated. We have kept only the most relevant term in its Fourier expansion, which is marginal; the other 
terms are irrelevant. Thus, the effective action takes the form: 

9 f 4*„(<VAtf J- l S± Z £j / ^2„p , , ,„ / J2 -2« 



S = J <* MW + y J d 2 xR<j) + wj tfxe-^f + ... (55) 

The . . . denotes other (irrelevant) terms in the Fourier expansion of the potential term which enforces integral 
heights; we have kept only the marginal term. Rescaling the field <j> 2 v /g0, the effective action takes the form: 

S= / d 2 x{\74>) 2 + ^ ~jL ] [ d 2 xR<t> + w[ d 2 xe~^^+... (56) 



16tt J SiiT^/g 
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7. When we compute an Ising correlation function with n spin fields in this model, only the term of order w n 
is non-zero, i.e. this correlation function has n insertions of the marginal operator e ~' l 't>V3/2 ^ ner6j we nave 
substituted g = 4/3), which we call a screening operator. 

The preceding logic makes it seem natural for correlation functions in the Ising model (and, in fact, a large class 
of models which have a height model representation) to have a Coulomb gas representation. It is, thus, helpful 
for understanding our results intuitively. However, it is not the most convenient way to derive the Coulomb gas 
representation for the conformal blocks which we need. For that, we use a more technical approach, described in 
the next subsection. We note that, although the two approaches are very similar, there is not really a one-to-one 
correspondence between them, although the results which we find in this paper strengthen the connection. 



B. Algebraic Approach 

The algebraic approach to the Coulomb gas takes, as its starting point, the action for a free boson with a background 
charge a a at infinity from step 5 of the previous subsection. This can be re-written as a total derivative term with 
imaginary coefficient iao- This total derivative term changes the energy-momentum tensor, thus shifting the central 
charge from its free boson value, c = 1, to c = 1 — 24«g. Since the added term is imaginary, the theory is not unitary. 
However, for certain values of ao, including the one relevant to the Ising model, the theory has a unitary subspace. 

This approach was introduced by Feigin and FuchsZ^, and developed for the minimal models^ by Dotsenko and 
Fateev in Refs. 66.8(J. The method was subsequently refined by Felder— , who both elucidated its BRST cohomological 
structure and extended the results to the torus. The advantage of Felder's approach is that it holds at the operator 
level, not merely at the level of correlation functions, allowing a more systematic description. This leads to a simple 
prescription which can be applied in a uniform manner. Thus, we adopt Felder's notation. The next few paragraphs 
are a short review of the procedure, whose full details can be found in Ref. l67l . 

The approach consists of taking a holomorphic free boson field (p(z), whose two-point correlation function is given 

I'.v 



( P^iM^a) ) = -2 log (z\ - z 2 ) . 



(57) 



This corresponds to Eq. (|A5I) with g — 1/4. This field can be used to construct the vertex operators of charge a, 
e 2ay(z) jf the ctj satisfy charge neutrality, YljLi a j = Oj then a collection of vertex operators has correlation function 



iaiip(zi) e ia 2 v(z2) 



N 



*jV¥?(ziv) \ _ 



IK* 

i<j 



(58) 



If 7^ 0> then this correlation function vanishes. This matches up precisely with Laughlin's plasma analogy, 

where the vertex operators in the CFT Coulomb gas formalism correspond to the particles that comprise the plasma, 
and charge neutrality must be obeyed (though, for the Laughlin states' plasmas, charge neutrality is achieved through 
a uniform background charge density). However, one can deviate from this simple Coulomb gas in two ways: (i) If a 
set of vertex operators violates charge neutrality, one can place an additional compensating vertex operator of charge 
2aa = —^2jOtj at oo to obtain a non- vanishing correlation function, and (ii) one can introduce "screening charges" 
which modify the vertex operators. As mentioned before, this will produce a unitary theory only for special values of 
ao, ctj, and screening charges. 

The minimal model CFTs with central charge c = 1 — 6 are denoted M[p,p') for positive integers p and p', 

where p > p' and the unitary theories are given by p = p' + 1. (The Coulomb gas formulation also applies to the 
non-unitary minimal models with p ^ p' + 1, so we leave these integers arbitrary in the following expressions.) The 
Ai(p,p') CFT has two screening charges given by 




(59) 



a possible charge at infinity of 



2ao = a_ + a + 



P-P 

Vpp 7 



(60) 
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and allowed vertex operator charges 



a nm = \ i (l-n)a- + \{l-m)a + , (61) 



for n = 1, . . . ,p — 1 and m = 1, . . . ,p' — 1. The vertex operator e lanmV ^ for these charges are used to describe the 
minimal model's primary field 4>( n ,m) with conformal scaling dimension 

U 1 I , \2 2 {jTVp - Tip') 2 - (P - P') 2 / r0 N 

ra ra ,m = 7 »a- +ma + ) - a = — , (62) 

4 4pp' 

which can be arranged into the conventional Kac table. These fields obey the identification <j>( n .m) = 4 l (p~n,p'~m)i 
so there are (p — l)(p' — l)/2 distinct fields in the theory, each of which can be represented by two distinct vertex 
operators. Of course, these vertex operators by themselves cannot correctly represent the minimal model fields, so 
screening operators 

Q± = j dze la±{p{z) (63) 

must also be introduced. 

In this way, one can generate the correlation functions of the minimal model CFTs as a Coulomb gas with screening 
charges. However, the specification of the contour integrals of the screening operators of Eq. is a crucial matter. 
For this, we follow Felder's prescription of combining the screening operators with the vertex operators to form 
screened vertex operators^ 

V^(z) = f[<f dw k f[I d Ul e la ^ {z) e la - v{wk) e la + v( - u '\ (64) 

k= l-!C k 1=1-1 Si 

where the screening charges' integration contours in Eq. (|64[) are taken to be concentric circles of radius \z\ centered 
at the origin (with the a+ contours Si inside the a_ contours Cfe), as shown in Fig. [TJ These contour integrals have 
divergences that must be regularized in some manner, i.e. either by an appropriate point-splitting at z or through 
analytic continuation. 

The full conformal block of the CFT operators is represented by the correlation function 

T a {z u ...,z N ) = (<j> (numi) ( Zl ) . . . (f> (nN , mN) (z N ) ) a = ( V n ™(zi) . . . V^ N {z N ) ) (65) 

for a set of screened vertex operators V^^ k {zk), where the indices n k and are chosen to represent the minimal 
model field 4>(n k ,m k ) a t ^k- The conformal blocks are labeled by the intermediate states in the fusion sequence of fields 



a = {(>i,Mi) (vn,^n)} 



(66) 
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which must obey the fusion algebra of the CFT, i.e. 



("fc-l,Mfc-l) 



(n k ,m k ) 



This means that 



l + \n k - v k -i \ <v k < min{n fe + z/ fc _i - l,2p- n k - v k -i - 1} , 
1 + \m k - < fi k < mm{m k + u k -i - 1, 2p' -m k - /ife-i - 1} . 

(n k + v k -i + v k ) mod 2 = 1, 
(m k + u k -i + fi k ) mod 2=1. 



(67) 



Here, we have (fo,//o) — (^n,^n) — (1, 1), indicating that the entire collection of fields must fuse to vacuum [the 
identity field / = 0(i,i)]- These different conformal blocks correspond to different choices of the numbers of screening 
charges r k and s k assigned to each screened vertex operator, which must therefore satisfy the conditions 67 



fj, k -l = m k + fi k - 2s k - 1 
V k -\ = n k + v k - 2r k - 1. 



(68) 
(69) 



When these rules are not satisfied, the correlation functions will evaluate to zero. It is important to recognize that 
these rules require that the sum of all Coulomb gas representation charges (vertex operator charges a nm and screening 
operator charges a±) is zero. 

One can alternatively represent the same conformal block as in Eq. (|65p by using the fact that 4>( n ^ m ) = 4>( p - n ,p'-m) 
and, in particular, that the identity field can also be represented by / = 0( P -i,p'-i)- With this in mind, it becomes clear 



that we should obtain the same conformal block if we replace V^ 1 ^ (z\) with 



l± (zi), require {v N ,n N ) = (1, 1) 



and (^o, /io) = (p—l,p'—l), and place a charge 2ao at infinity. In this case, the sum of all Coulomb gas representation 
charges together with the charge 2ao at infinity is zero. 

We now focus on the Coulomb gas representation of the c = 1/2 Ising CFT, which is the -M(4, 3) minimal model 
and has three primary fields /, ip, and a. The two screening charges for this CFT are given by 



V3 
2 ' 



2 



The six vertex operator charges are constructed out of a_ and a+ , according to 



1 



m) a_|_, 

It is convenient to put these together into the "Kac table" : 



1,2,3, m = 1,2. 



(70) 



(71) 



2 


l 


i 


1 


1 





4 


V3 
2 


m/n 


1 


2 


3 



(72) 



The columns of the table are labeled by the index n = 1,2,3, and the rows by m — 1,2. The entries of the table 
are the charges a nm of the vertex operators. They represent the operators of the c = 1/2 Ising CFT according to the 
identification 



2 


4> 


a 




1 


i 


a 




m/n 


l 


2 


3 



(73) 



Here / is the identity field, and just as in Eqs. (j2"2"j) and (|2"3")l . a is the dimension 1/16 operator, and tp is the dimension 
1/2 operator of the c = 1/2 Ising CFT. 

We can now examine in detail a concrete example of how the same conformal block can have several different but 
equivalent representations using the screened vertex operators. Consider, for example 



>(*)#>)> = -. 

z 



Since both V12 and V31 correspond to "0, this correlation function can be represented in three different ways: 

1 

z 



(iffW^o)) 



(74) 



(75) 
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< V$(z) V™(0) )=Uu Z ' 3 4 , (76) 

J (Z-1()3M3 

and 

< vffco <(o) > = hw, ld W2 zl ^" W2)§ 3 ( 77 ) 

^ ^ (z — W±) 2 (z — W 2 ) 2 U>2 

It should be clear that all three methods give the same answer, 1/z, up to an overall unimportant constant (as may 
be verified by making the change of variables: £ = uj z and = Wk/z). The correlation function in Eq. (|75j) has total 
charge — 2ao — 1/2-^/3 (which is canceled by a charge at infinity) while the other two correlation functions, Eqs. ([76)) 
and ([77]) have the total charge 0. 

Now we can use these techniques to represent the conformal block corresponding to the Pfaffian as 

Pf (-^—) = ( ^1) • ■ • i>{z N ) ) = < Vi?{z x )V$?{z 2 ) . . . V£Wi)^TM ) , (78) 

(where N is even). This is not the only way to construct this conformal block, but it is the most convenient for 
subsequent generalizations. 

Now consider a conformal block with four a operators (which correspond to four quasiholes). There are two such 
conformal blocks [as we saw, for instance, in Eq. (I26|) ]. which we denote by: 

Takn^zi) = MmM^M^M^Mzi) . . . ^{z N )) a (79) 

where a = 0, 1 corresponds to the block in which the first two a fields fuse to I or ip, respectively. We can represent 
7-0(77^; Zi) in the following way: 

7b(^; Zi) = < V^( Vl )V^(V2)Vinm)V^°(m) W(*0 . . . V^(zn-i)V^(z n ) ) . (80) 

This representation mirrors that of Eq. (175)) in that it only uses vertex operators from the m = 1 row of the Kac table 
in Eq. ([72"]) . The total charge of all the operators involved in Eq. ([50")) is equal to zero. Furthermore, the total charge 
of the first two screened vertex operators is also zero, V^ V^ ~ /, which is the reason for the identification of this 
Coulomb gas correlation function with the c = 1/2 Ising conformal block J r o(Vni z i)- If we wish, instead, to compute 
.7-1(77^; Zi), then we need a Coulomb gas correlation function in which the first two screened vertex operators have 
total charge \/3/2 corresponding to the ip field: 

^( Vfl ; Zi) = < V^irn^^V^imW^M Vi 1 °(z 1 )V 3 ° 1 °(z 2 ) . . . VifizN-JV™^) ) . (81) 

Since the screening operators are attached to the first two vertex operators, rather than the first and third, the 
construction Eq. (|81[) can be interpreted as simply a different choice of contour for one of the screening operators in 
Eq. 

We note, for later use, that we can also represent J-"o and T\ in an alternative way: 

Jb(^; Zi) = { V 22 (vi ) V21 (V2 ) V21 (»?3 ) V21 (174 ) V™( Zl )Vf?(z 2 ) . . . Vi?{z N ^)V™{z N ) ) , (82) 

F^Zi) = { VgMVjfMVgMVnfa) V™( Zl )V™(z 2 ) . . . V^{z N . 1 )V^{z N ) ) . (83) 

Unlike in Eqs. (150)) and (f5T)) , the total charge of the vertex operators involved in Eqs. (152)) and ([53")) is equal to 1/ (2 v3) 
(which is another representation of the identity) . 

Finally, we can also construct conformal blocks with any even number n of a fields (corresponding to wavefunctions 
with n quasiholes), e.g. 

^ (0 ,o,..,o)(^; «») = < vj?(m)vff(v») . . . v a VWi)v&V) v^C* Wte) ■ ■ • vfiWi^TM ) (84) 

= < VgMVgiTriVtfMVgfa) . . . V2\°(Vn-i)V2 1 °(r 1 n) V™ ( Zl )V™ (z 2 ) . . . V^{z N . l )V^{z N ) ) (85) 

The subscript (0,0,..., 0) denotes that this is the conformal block in which the first and second a fields fuse to /, the 
third and fourth a fields fuse to /, . . ., the (n — l) th and n th a fields fuse to /. 
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Similarly, the Coulomb gas construction gives the general conformal block 

JW,..,ir n/a )(»fc; *i) = ( II [^"^'Wi W'V)] Vi 1 °(z 1 )V 3 a 1 °(z 2 ) . . . V™{z N ^)V™{z N ) \ (86) 



,/2 



faO^faO II [^i"^ 1 ' Wi)^2i 'V')] t&VMteO • • . t$W-i)V£W) ) , (87) 



i=2 



in which the 1 st through 2j th er fields collectively fuse to J if Wj = and to t/> if 7Tj = 1, and where 7To = 7r n / 2 = 0, 
indicating the overall parity constraint that the a fields must collectively fuse to / since there are an even number 
N of ijj fields. This is presented in the "standard basis," where fusion channels are specified by fusing in the anyons 
one at a time from left to rights. For the Ising CFT, we can trivially transform between the standard basis and the 
qubit basis 

-7-qubit _ x-standard 

• A (pi,p 2 ,...,P„/2) ~~ • / "(^1^2,---,7T Il /2) \ 00 ) 

in which the (2j — l) th and 2j th a fields fuse to / if pj — and to tp if pj — 1, by simply using the conversions 

^PA^mod2 (89) 

Pj = (ttj — 7Tj_i) mod 2. (90) 

Since this is a trivial change of basis (i.e. it is just a different way of presenting the subscript label), we can interchange 
between the two freely. For the purposes of describing conformal blocks using the Coulomb gas formalism, the standard 
basis is more natural. For describing the explicit evaluation of the conformal blocks using bosonization methods, the 
qubit basis is more natural. Henceforth, we differentiate the use of these two bases through context. 

A similar expression can also be used if the number TV of electrons is odd. Specifically, for N odd, one would use 

/ n/2 \ 

-W 2 ,...,.„ /2 )(^) = (jl [^~*'- 1 'Wi)t£ ,,< W] V^(z 1 )V^\z 2 )V^(z 2 )...Vi 1 (z N ^)Vl (z N )J(91) 

= ( v™{ m )v^\ m ) |j [vir^imi-iWS' ^)] v^(z 1 )vi 1 °(z2)v 3 ° 1 Q (z 2 ) . . . v$t(z N -i)vg?(z N ) \ (92) 

with 7To = and n n / 2 — 1, which indicates that the n as have overall fusion channel We note that the number of 
screening charges in both Eqs. (|86|) and ([9T]) is TV + j. 

The explicit expressions for correlation functions such as Eq. (|84p involve products of powers of differences of 
coordinates, and integrals over some of them, as in the simple examples of Eqs. ([TBI and ([77|) . This has a reasonably 
similar structure to the Laughlin states, such as Eq. (|46|) . so it brings us closer to the goal of constructing an effective 
plasma describing Eq. (jTTj) . 



VI. PLASMA REPRESENTATION FOR THE NORM OF THE GROUND STATE WAVEFUNCTION 

Using the preceding expressions for the conformal blocks to construct the overlap integrals Eq. (j^TI) . we see that 
they do appear superficially similar to the plasma construction of Eq. (jlS)) . The difference is that the screening 
operators need to be integrated over their holomorphic and antiholomorphic coordinates along some specially chosen 
contours. As a result, Eq. (|37[) no longer takes the form of the partition function of a classical plasma. In what 
follows, we construct the overlap integrals in a slightly different way which leads to an expression which does take the 
form of a classical plasma's partition function. For this, we crucially utilize the method invented by Mathur in Ref.l68l 
of relating expressions involving products of holomorphic and antiholomorphic screening charge contour integrals to 
expressions involving 2D integrals over screening charge positions. (We review this method in Appendix |B|) 

We will begin by considering the case with no quasiholes, ie. the ground-state wavcfunction. We will construct 
a representation of the norm of the ground-state wavefunction which takes the form of the partition function for a 
classical plasma. We begin by ignoring the charge part of the wavefunction and focusing on the Pfaffian: 

»(-!_) 



(93) 
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In order to represent this as a plasma, we take the conformal block represented by Eq. ([78]) . and multiply it by its 
complex conjugate. Then, instead of integrating the screening operators over the contours in the complex plane of 
their respective holomorphic and antiholomorphic coordinates, we integrate the screening operators over the entire 2D 
plane. 

To see why this procedure is valid, we first consider the expression 



Pf 



1 



Z% Zj 



= < Vg^VSte) . . . V^(zn-i)V^(z n ) ) 



dwi <p dw 2 e la3iv{zi) e la - v{wi) e ta - v{w2) e la3iv(z2) x... 
C Z1 Jc zi 



. . x 



dwN-l 



q] Wn e ia3iip(z N ^ 1 ) e ia-.<p(w N -i) e iu-ip(w N ) e ia 31 ip(z N ) \ ^g^\ 



where C x is used to indicate a contour of radius \x\ centered on the origin (with appropriate regularization, i.e. taking 
contours at the same radius to be infinitesimally concentric and point-split at the Zi coordinates). To obtain the norm 
squared of this wavefunction, we multiply this expression by its complex conjugate: 



Pf 



1 



dwi w (iir-> 



dyj 2e ia3iv(zi) e ia^v(w 1 ) e ia^ip(w2) e ia 31 ip(z2) x 

dWN e ia 3iv(zN-l) e ia-<p(w N - 1 ) e ia-tp(w N ) e ia 3lV (z N ) \ _ ^gg-j 



dw 



N-l 



Evaluating the correlation functions of vertex operators, noting that 031 = — a- = V3/2, we obtain 
1 



Pf 



2 N/2 

=n 



dw 2 k-i f dw2k-i 



j> dw 2 k f dw 2 k Yl 

•' C '2k-1 JC *2k-\ i<j 



N N N 

Wi - wA 3 TT \wi - zA~ 



n 

hi 



(96) 



It is important to emphasize that Wi and Wi are independent variables in this expression, so terms such as \u>i — Wj\ 

3/2 3/2 

should really be understood as shorthand for (wi — Wj) (wi — Wj) ■ Retracing Mathur's steps, as explained i 
Appendix [Bl we rewrite the product of Wi and Wi contour integrals in Eq. (|96p in terms of 2D integrals: 



jirf 2 Wfc Yi k - w /\ 3 n \ wi ~ z i\ 3 n \ zi ~ 2 ji 3 

k i<j i,j i<j 

N/2 N N N 

= Y\f dw 2k-i f dw 2 k-i f dw 2 k f dw 2 k Y\_\ w i - w jf Y\.\ Wi ~ z ^~ 3 Yi 



(97) 



Therefore, we can write the square of the Pfafhan in the form: 



Pf 



N 



N 



n^n 



Wi — Wi 



N 

n 



N 



Wi — z 



1<J 



Yi. \ Zi z i 

i<j 



(98) 



Note that the right- hand-side of this equation is divergent as any Wi approaches any Zj. It can be made well-defined 
by analytic continuation. In other words, we define this expression by evaluating the integral 



„ N N N N 

/ yi d2w k n ~ w 3 i 3a n \ m ~ z i\~ 3a n 

k i<j i,j i<j 



3a 



(99) 



for a < 2/3, where the integral is convergent, and analytically continuing to a = 1. This analytic continuation gives 
the right- hand-side of Eq. (|98|) . As we will discuss, the associated plasma does not go through a phase transition as 
a is varied from a < 2/3 to 1, so the right-hand-side of Eq. (|98l) is a useful representation of the left-hand-side. 

If, instead, we modify the right-hand-side at short distances by, for instance, introducing a short-ranged repulsion 
(e.g. a hard-core cutoff), then the right-hand side will be modified for — > Zj but will be unchanged at long-distances. 
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This will produce a wavefunction in the same universality class as the Pfafhan. However, rather than introduce a 
cutoff and work with a modified wavefunction, we prefer to define Eq. (1981) by analytic continuation, as described 
earlier. 

Now we can interpret the norm of the Pfaffian in terms of a two-component plasma. Specifically, we can write 



Pf 



1 



N 



IP 

k=l 
N 



Zk 



Pf 



1 



N 



= / Y[d 2 z k d 2 w k e'^ T = e~ F ^ T 



k=l 



N 



N 



$2 = -^Q 2 logK-^| + ^Q 2 logK-^|-^Q 2 log| 



*j i ) 



(100) 
(101) 



1<J 



i<3 



where T — Q 2 /3. Now $2 is the 2D Coulomb-interaction potential energy for N charge Q particles at Zi and N 
charge — Q particles at Wi. Thus, F 2 is the free energy of a classical 2D two-component plasma of charges ±Q at 
temperature T. (We use the subscript 2 to indicate the two-component plasma.) Again, we can let T take any value 
as long as Q is adjusted accordingly. One convenient choice is to take T = g and Q — \/3g. 

It is known that such a two-component plasma with coupling constant T — Q 2 JT is a screening fluid for T < 4, 
i.e. the condensation temperature is T C2 = Q 2 /4, but that it needs a short-ranged repulsive interaction, such as a 
hard-core cutoff, in order to be stable against collapse into neutral bound pairs for T > 2 (i.e. T < 2T C2 )^~—. For 
r > 4 all particles are bounded into neutral pairs, while for V < 2 all pairs are broken. The Pfaffian wavefunction's 
corresponding plasma is precisely in the range 2 < T < 4 where it is a screening fluid as long as a short-ranged 
repulsion is introduced. This fits with the preceding discussion regarding the need for a short-distance repulsion or 

analytic continuation, and is intuitively clear from the fact that Pf^ -^j^- j diverges as Zi — > Zj. We discuss the 

screening properties of this plasma in more detail using field theoretic methods in Appendix [Dj Its Debye screening 
length can be estimated (see Appendix [E| to be £ 2 — (12tt«/)~ 1 ^ 2 , where rif is the electron density. 
Adding the charge part of the MR ground-state, we have: 



\*(zi,...,z N )\ 2 = 



Pf 



1 



2 N 

m 



-h EM 2 



N 



N 



|2M 



N 



n 



i<3 



N 

n 

i<j 



z i\ 



2M+3 



-3 T, N.I 2 



(102) 



This expression is antisymmetric under exchange of z% with Zj while holding Zi and Zj fixed. It is also of degree 
M(N — 1) — 1 in any of the ZiS and degree M(N — 1) — 1 in any of the ZiS. Indeed, there is a unique polynomial 
satisfying these properties, so it is clear that once the right-hand-side is computed by analytic continuation, it will 
give the squared modulus of the MR Pfaffian ground-state wavefunction. 

Now we can write the norm of the MR ground-state wavefunction in terms of a classical plasma by writing 



4>i 



4>9 



r N f N 

/ Ip 2 ^ \*(zu...,z N )\ 2 = \{d 2 z k d 2 w k e-^+^l T 
J fc=i ^ fe=i 



-F/T 



k=l 
N 

5Z<9? log |^ 

i<j 
N 

Ql lo s \ w -> 



z A + ^Ly \ Zl \ 2 
31 AM 11 

i=l 



N 

E 

hi 



N 



Ql log \wi -Zj\-y^Ql log \zi - zj\ , 



(103) 
(104) 
(105) 



where T — g, $1 corresponds to the 2D Coulomb potential for N charge Qi = ^/2Mg particles at Zi in a uniform 
neutralizing background of charge density p\ = — 27r ^^ , and $2 corresponds to the 2D Coulomb potential for N 

charge Q 2 = V^g particles at Zi and N charge —Q2 = —\/3g particles at Wi (and no neutralizing background charge 
density). Thus, F is the free energy of a classical 2D plasma (at temperature T) in which its particles can carry 
charges corresponding to two independent types of Coulomb interactions, differentiated using the subscripts 1 and 2. 
In particular, the plasma described here consists of N particles at Wi carrying charge — Q2, N particles at z.j carrying 
charge Q\ and Q2, and a uniform background of charge density p\ = — 2 Ji,ii' 2 ano - P2 = that neutralizes the charges 
of type 1. 
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When plasmas 1 and 2 are independently in the screening liquid phase for the corresponding values of T, Qx, and 
Q2, we expect that the combined plasma should also be in the screening liquid phase (except, perhaps, as these 
parameters become close to their critical values for either of the two plasmas) . A recent numerical study 5 -^ of this 
combined plasma for M — 2 found that the behavior is very similar to the M = case (i.e. the two-component 
plasma) in that it is in the screening phase for Q\jT < 4. In particular, this study verifies the MR state's analogous 
plasma (which has Q\jT = 3) is in the screening phase for the most important case M = 2. (For additional results 
on similar plasmas in the context of vortices in multi-component superconductors, which support the idea that the 
combined plasma is likely to be in the screening liquid phase, see Refs. I88l490l .) We estimate the Debye screening 
length of this plasma (see Appendix [Ej to be 

/ M \ 1/2 
e=[ ; Ib. (106) 

For M = 2, this gives i « 1.2i B - 



VII. PLASMA REPRESENTATION FOR THE TRACE OF THE OVERLAP MATRIX 

The situation gets more complicated when we turn to wavefunctions with multiple quasiholes. We would again like 
to be able to treat quasiparticles as test charges in the analogous plasma. However, this is not as straightforward to do 
as for the Laughlin case. There are multiple degenerate wavefunctions (corresponding to multiple conformal blocks) 
in such a case. These different conformal blocks are distinguished in the Coulomb gas formalism by the location of 
the screening charge operators' contours. Thus, if we exchange a pair of screening contour integrals, one holomorphic 
and one antiholomorphic, for a 2D integral too naively, we would elide the distinction between the different conformal 
blocks, which would clearly be incorrect. 

Thus, we must proceed with greater caution. To do this, it is useful to recall that the c = 1/2 Ising CFT with its 
conformal blocks is but a mathematical tool to construct the correlation functions of the Ising model at its critical 
poin t 58 ' 91 . These correlation functions are real, not complex, and they depend on the 2D coordinates of the operators 
of the Ising model, not just on the holomorphic part of these coordinates. In particular, consider a correlation function 
of four Ising spins (order operators) a, as well as N Ising energy operators e: 

( (7(771, 771 )ct(?72, 772)^(773, m)cr(m, fj4,)e(zi,zi) . ..e(z N , z N )) . (107) 

Note that these are non-chiral operators. For instance, e = ijjip, where ip is the chiral Majorana fermion field 
introduced earlier and ip is its antiholomorphic counterpart. This correlation function can be written in terms of the 
two conformal blocks, JF and T\. These conformal blocks are the chiral part(s) of the correlation function Eq. (|107p . 
which we denoted in the previous section as: 

Faiv^Zi) = (< :r (?7l)cr(772)f7(?73)cr(?74)l/'(zi)...-0(2Ar)) a (108) 

Note that these are now chiral operators cr(r]^) and ip(zi). The subscript a = 0, 1 denotes whether the first two a 
fields fuse to I or ip, respectively. The explicit forms of To and T\ are: 

= / 77 13 7724 \ 8 1^ / f ± ^— - % \ (109) 

Vf?12??23»734774l/ y/l±y/l-X V 

where ^(13) (24), ^(14) (23) are defined in Eq. (|F11[) . Note that T$ and T\ are clearly multi-valued functions; they 
transform under the braiding of coordinates in exactly the same way as the functions ^0 and of Eq. (|26[) . up to 
an overall phase [which is due to the c = 1 CFT present in Eq. (J5SJ)]. Indeed, ^o^i were constructed by multiplying 
in Eq. (|109[) by a Laughlin wavefunction-like factor coming from the c = 1 CFT. 
The antiholomorphic part of the correlation function is similarly given by !Fq and J-±. However, the non-chiral corre- 
lation function ( 17(771, fji) . . . 17(774, 774) e(z\, z\) . . . e(zpj, Zn) ) must combine holomorphic and antiholomorphic sectors 
in such a way as to be single- valued. There is a unique way to do this, which is the trace: 

(ct(t7i,t7i)ct(?72, 772)^(773, 77 3 )cj(774, 774) e(zi,z~x) . . . e(z N , z N )) = Tbfe; Zi)Fo(Vn'> *i) + FiiVni z t)^i(v^ ( 110 ) 

Indeed, this is the only combination of the conformal blocks which is single valued as 771 and 771 = 77^ are taken all over 
the complex plane, and similarly for the other 77^ and 77^,. This may be checked by using the analytic continuation 
properties of Jq and T\, which are exactly the same as Eqs. (|3U| . (|3"Tj) . and (|3"2")l for >f and \I/i (up to the overall 
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phase, which obviously cancels between holomorphic and antiholomorphic terms anyway). This expression is also 
real, as expected for a real correlation function of the Ising model. 

Since the sum of the squares of the four-quasihole conformal blocks is single-valued, we can form a plasma repre- 
sentation for the sum of overlap integrals: 

f N 

/ Y]_ d 2 z k [Jb(r? M ; 2j)Jb(^; ^i) + J"i(r?^; z^Ti^; 2*)] . (Ill) 

^ k=l 

We cannot do this for each of the individual terms in this sum. In order to express Eq. 1 1 1 II in terms of a classical 
plasma, we begin with the Coulomb gas representation for the conformal blocks of Eqs. (|80p and (|8ip [we could 
have equally well chosen the representations in Eqs. (|82p and (|83p . but this choice is more suitable for subsequent 
generalizations, as we will see later], and multiply them by their complex conjugates. The conformal block .Fq(?7^; zi) 
is precisely the expression which we defined in Eq. (|80| : the conformal block in which the first two as fuse to /. 
Written explicitly in terms of the vertex operators, this is: 

Mv„ Zi) = { Vi?( Vl )V 2 ?(V2)V™( m )V™(VA) V™( Zl )V™(z 2 ) . . . t$W_i)t$W) ) 

d,Wi e ia av(vi) e ia-f(wi) e ia2iip(m) A e ia2itp(r)3) e ia-ip{w 2 ) e i(X2itp(ri4) 

C r!l J 

x l dw* £ dwAe ia31 ^e ia -^ w ^e ia -^ Wi h io131 ^ x ... 



3 

C„. JC 



...xi dw N+1 J> dwN+2e i a 2i-p{zN-i) e ia- V {wN+i) e ia. V (w N+2 ) e ia 3 Mz N ) \ ( 112 ) 

We multiply this expression by its complex conjugate 

To{%',zi) = (<f dw 1 e ta2lvini) e ta - ip(lSl) e ia21 ^ rl2) I dw 2 e iajl ^ s ) e ia -^ a V aai ^ 

\JC fjl JCfj 3 

x I dw 3 I dme iasMSl) e^- v ^ 3) e ia -^ iE4) e ia31 ' p ^ 2) X ... 

...X 1 dw N+1 <f dffi ^ +ae fasi*>(«»r-i) e to-p(«w+i) e *a-^(ffl w+ ,) e iasi^)\ ( U3 ) 

In Eqs. (|112|) and (|113p . there are N coordinates (electrons), four coordinates rj^ (quasiholes) , and N + 2 coordinates 
w a (screening charges) . The correlation functions of vertex operators in Eq. (|112l) and (|1 13[) can be evaluated using 
Eq. (58]): 

N+2 , . N+2 N+2 4 N+2 N 



Fairi^Zii^^^Zi) = Y[ j> dw c j> dw c (w a -w 6 ) 5 JJ J\ {w a - VnV 11 JJ JJ K - z») 

c— 1 * a<6 a— 1 fi— 1 a— 1 z— 1 

JV 4 iV 4 

x jj (zi - Zj )% jj j] - 2l ) f n ^ - 

2<j /.i— 1 2= 1 /^<z/ 

7V+2 iV+2 4 iV+2 iV 

n k n n 1 n n ( 



N+2 JV+2 4 N+2 N 

xn^-^nn - ^r 1 n n 

a<b a— 1 fi— 1 a— 1 z— 1 

TV 4 TV 4 

* n> - %) f n n - *) f n ^ - ■ ( m ) 

j<J (U=l i=l f4<f 

In these expressions, the appropriate choice of integration contours (which we left implicit here) tells us that we are 
computing J-'oi'Hu'j zi)To{fi tl \ zi). However, by choosing a different contour for one of the screening charges in Eqs. (|112p 
and (|113p . as per Eq. (l8Tj) (specifically, if the contour C 2 corresponding to w 2 was a circle of radius \r] 2 \ rather than 
radius 1 773 1 ) , we would obtain J-\(r]^ zi)fi{f]^]zi) instead. To obtain the non-chiral correlation function, we should 
add the right-hand-side of Eq. (II 14[) to the corresponding expression for .7-1(77^; zi)T\ {fj^; zi), with these different 
integration contours. 
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Instead, following Mathur 68 once again, we replace the integrations over pairs of contours by integrations over the 
plane as described in Appendix [Bl This replacement gives us neither J-cife; Zi)To(fj^; zi) nor Fife; 2i)F"ife; zi) but, 
rather, the combination Fbfe; 2i)Fofe; zi) + Fife; Zi)Fife; zi). Thus, we obtain 



Fo(Vri z i)Fo(Vrf%) +Fl(Vri^)Fl(V»':Zi) = / II d2yj c W_ \w a ~ W b \ 3 JJ J]j U>a ~ ?^ |~ 5 JJ J^[|w a -Zi|" 3 

c— 1 a<6 a— 1 fi— 1 a— 1 i— 1 

iV 4 W 4 

xn^-^/nnj^-^ 1 n fe-^i f • ( n5 ) 



The reason that the particular combination Fo(?7^; z^Fofe^ zi) + Fifej; Zi)J 7 i(fj tl ; zi) appears on the right-hand-side 
is, as shown by Mathur—, that when the contour integrals are replaced by 2D integrals, as described in Appendix IBl 
this has the effect of computing a sum of holomorphic and antiholomorphic conformal blocks, such that the entire 
combination is single- valued as a function of all variables. 
We now define 

f N 

G a,p(v^,V^)^ W^ZkFairi^z^Fp^^Zi). (116) 

k=l 

We denote this overlap matrix as G^ to distinguish it from the closely related G defined in Eq. (f3"T|). which is the 
overlap matrix of the MR wavefunctions If we take the integral of Eq. (II 15|) over the coordinates Zi, we obtain 

TrG-*^" = Gq q + G^i 



/ II d2zk [^oiv^i z;)F fe; 7 -i) + Fife; «i)Fj,fe; zi)] 
^ fe=i 

/N+2 N+2 N+2 4 N+2 N 
n d 2 w c n \ Wa - Wb \ 3 n n k ^r f n n 
1 ^- ^ ,-, 1, 1 „ 1 A 1 



\w„-zi\ 3 



a<6 a— 1 /i— 1 a— 1 i=l 

N 4 JV 4 



Comparing with Eq. (|A4|) , we can rewrite this in terms of the partition function of a plasma 

„ N N+2 

TrG^ = / TT d 2 z k TT d 2 w c e~^-l T = e~ F ^ T (118 



. iV N+2 

n rf2z fe n d2wc e_ 

fc=l c=l 


-*a/T _ e -FilT 




Af+2 

^2 Q 2 \og\w a -Wb\ 

a<b 


N+2 4 Q 2 

+ EEy lo sK" 

a— 1 ju= 1 


iV+2 TV 
a— 1 2—1 


N 
i<j 


/i— 1 i= 1 


4 o 2 



-Vu\, (119) 

where, for temperature T = g, $2 is the 2D Coulomb-interaction potential for N particles of charge Q — ^/3g at 
positions z% and N + 2 particles of charge — Q at positions u> a , in the presence of four fixed test particles of charge 
Q/2 at positions 77^. This plasma obeys overall charge neutrality, as can be seen by adding up all the charges. As 
previously mentioned, it is known that the 2D two-component classical plasma comprised of particles of opposite 
charge +Q is in the screening fluid phase for T — Q 2 /T = 3, though a short-ranged repulsion (e.g. a hard-core cutoff) 
is needed. (We discuss the screening properties of this plasma in more detail in Appendix |D|) . Since this plasma 
screens, the free energy F2 in Eq. (|118|) is independent of the positions rj^, as long as they are farther apart than the 
screening length I2 of the plasma. This proves, for the case of n — 4 quasiholes, that 

TrG* = 2C 2 + Oie-^-^/* 2 ), (120) 

where C2 is a constant independent of ry M . 
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We now turn to the overlap matrix G of the MR wavefunctions The wavefunctions an d *i differ from 
and JF\ by an additional c = 1 correlation function, as is clear from Eqs. (|26p and (|23[) . This additional correlation 
function is straightforward to calculate [as before, use Eq. (|A6| with g = 1/4] 



.3 



4 4 N N _ _j_ A | ,2_ i ■+ , ,2 

n, . _i_ TT TT / x I TT , x M 8M £-r I^Ml 4 A» 1**1 

- »*o 4M n n - n (* - ^) e m ^ * =i ■ ( m ) 

H<v fi— 1 z— 1 i<j 

Consequently, the analog of Eq. (|1 1 7[) for these functions is given by 
TrG = Go,o + Gi,i 

N 

Yl d 2 z k [*ofe; Zi)^o(v^, Zi) + 2i)*i(?7 M ; z,)] 

fc=i 

n ^ 2zfc / 1 1 n ~ Wb ' 3 n n ~ 

fe=l c=l a<b a=l ju=l 

AT+2 AT JV 4 iV 4 

n n k *r 3 n i* - ^i 3 n n I*?*. - n - ^i 1 

a=l i=l »<j fi = l i=l /i<i/ 

W 4 N 4 j_ A , |2 _i " , ,2 

ki-^-i nni^-^i n i^-^i 2M e *=* «=* . (122) 

i<j /x=l i=l 

Similarly, this can be interpreted in terms of a plasma for which there are two independent Coulomb interactions, 
denoted using subscripts 1 and 2, by rewriting Eq. (|122[) as 

„ iV JV+2 

TrG = / [] d 2 z k [] d 2 w c e -(*i+*0A = e"^ (123) 

fe=l c=l 

4 q2 4 AT Q2 TV 

fi<v fi—1 i—1 i<j 

^.2 4 n 2 n 



X 



8M 2 ^ " Ml 4M 

JV+2 W+2 4 n 2 N+2 N 

$2 = l0 S l Wa _ + E E ~2 l0g \ W * ~ ^1 + E E Q 2 log ha " *| 

a<6 a— 1 a— 1 i—1 

^ 4 ^ 2 4 „ 2 

- £ lQg ^ _ ^ _ log |^ - Zi \-J2^ log IfM - . ( l25 ) 

i<j i—1 n<v 



where, for T — g, we have $i corresponding to a 2D Coulomb potential for JV charge Qi = \J2Mg particles at Zi 
and four fixed test particles with charge Q\/2M = \J g/2M at rj^, in a uniform neutralizing background of charge 
density pi = — t^mp? ' an( ^ corres P onal ing to a 2D Coulomb potential for JV charge G2 = V^ff particles at 2^, JV + 2 
charge — Q2 particles at w a , and four fixed test particles with charge Q2/2 at 77^. Hence, this plasma consists of JV 
particles (corresponding to the electrons) at Z{ which carry charges Q\ and Q2, JV + 2 particles (screening operators) 
at w a which carry charges — Q2, four fixed test charges (quasiholes) at 77^ which carry charges Q\/2M and Q2/2, 
and a uniform neutralizing background of charge density p\ = — 2 ^mP ( an d P2 =0). As previously mentioned, we 

expect such a plasma to be in the screening phase for roughly T > T Cl ,T C2 , where T Cl = Q^/140 and T C2 = Q|/4 are 
the critical temperatures above which plasmas 1 and 2 are individually in their screening fluid phase. Therefore, this 
plasma at temperature T — g with Qi — y/2Mg and Q2 — y/Sg should be in the screening phase for M < 70. This 
has been numerically confirmed^ for M = 2. When the plasma is in the screening phase, the free energy F will not 
depend on the positions of the test charges 77^ , as long as their separations are larger than the screening length t of 
the combined plasma. 
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Thus, for sufficiently small M, we have proved that the trace of the overlap matrix G [defined in Eq. (|37|) ] for n — 4 
quasiholes is an ^-independent constant for large separations, or 

N 

TrG= £ / Hd^k^^z^^z,,) = 2G + 0( e -l'"^). (126) 

a=0,l"' fc=l 

Hence, we have established that both the trace of G, which includes the charge sector as in Eq. (|122j) , as well as the 
trace of G^ , without the charge sector, as in Eq. (|117l) . are constants. 

The preceding derivation can be generalized to an arbitrary even number n of quasiholes, for which the formulas 
analogous to Eqs. (fTT8| . (|T23| . ([P24]) . (|T25|) are 

N q-l „ JV w +f 



TrG* = f f[d 2 z k Y, Mv^Fa^Zi) = f f[d 2 z k [J d 2 w c e~* 2 / T = e - F ^ T (127) 

•* k=l q=0 k=l c=l 

r N q-l . N N+% 

TrG = / \{d 2 z k ^ a {%;-z^ a {ri»,z,)= / \{d 2 z k [] <?w c e "(*x + * 2 )/T = & -F,T (12g) 

7 i ~ — n *> l — i 1 



fc=l a=0 " fe=l c=l 

™ Q 2 " W Q 2 N 

$1 = " £ 4^2 lo § - ^1 - S S 2M l0g ~ Zj| ~ ^ °' bg |Zi ~ Zj 

n 2 n ra N 



■ MP E 1*1' + SEW (129) 

/l — 1 2—1 

$2 = - £ Ql^g\w a -W b \+ ^ £ "^logl^a ~ ^1 + £ XI ^2 log ko - Zi\ 

a<b a—1 j.L—1 a—1 i—1 

N n N n 2 U (I 2 

-J2 qI \zi - zj\ - J2 £ t lo § - *i - £ t log ^ _ ^ ' (130) 

where q = 2§ _1 . The sum over a can be replaced by a sum over 7ij = 0, 1 or pj = 0, 1, for j = 1, . . . , n/2 with 
the parity constraint n n /2 = (^2jPjj mod 2 = (for N even). Summing the diagonal product of holomorphic and 
antiholomorphic conformal blocks over all conformal blocks, we obtain the single- valued expression on the right-hand- 
side. 

The arguments discussed thus far are a carefully worked-out version of the arguments presented in Ref. fiM Their 
main drawback is that they do not prove Eq. (|41|) . They only prove the weaker statements 

9-1 „ N 



TrG^ = £ / f[d 2 z k T a (^;z l )T a (r lfl ;z z ) = qC 2 + 0{eT^^/ k ) (131) 

Q = 0"^ fe=l 

9-1 , N 

TrG = £ / n**k9 a (T}S,Zi)* a (VtfZi) = <lC + 0(e-^-^\/t) (132) 

Q=0^ fc=l 

which is necessary but not sufficient, with one exception, for the (nontrivial part of) Berry's connections Eq. (|34[) to 
vanish. In the next section we extend the proof to show that the stronger statement Eq. (|4Tj) is true. 

The noted exception is the case of two quasiholes. Since there is only a single conformal block in this case, equal 
to either of the Coulomb gas expressions 

Jo(^5 *) = { V^)V£°M V 3 \°( Zl )V™(z 2 ) . . . V^{z N -!)V^(z N ) ) 

= ( V£V)V£°M Vi?{ Zl )V™{z 2 ) . . . V^{z N . l )V^{z N ) ) , (133) 

the overlap matrix is a 1 x 1 matrix which is equal to its trace, and so, for the n = 2 quasiholes case, we have 

N 



G o,o = [f[d 2 z k T (fj^Zi)Mv^Zi) = C 2 + 0(e-l""-"-l/ £a ) (134) 
J k=i 

N 

Go,o - J II d2z ^o(v», Zi)Mw Zi) = C+ 0(e-l""-^l^). (135) 
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Thus, as concluded in Ref. |45| by the same logic, when there are only two quasiholes, the effect of a counterclockwise 
exchange is the accrual of a statistical phase exp [in — |)] and an Aharonov-Bohm phase exp (— i^jlj^r)- In 
fact, we should be more precise: This is the phase that is accrued when there is an even number of electrons in the 
system. When there is an odd number of electrons, one can repeat these steps using 

MV„1 Zi) = < V™{ m )V™{ m ) V™( Zl )Vi?(z 2 )V™{z 3 ) . . . V^{z N ^)V^{z N ) ) 

= ( V 2 ° 2 (vi)V 2 \ (v2) V^(z 1 )Vinz2)V 3 °^(z 3 ) . . . t£W_i)V£W) ) , (136) 

which has the two quasiholes fusing into the ip channel. The overlap matrix is again a 1 x 1 matrix which is equal to 
its trace, but in this case it gives: 

N 

G M = / n<Pzk?i(fi»,Zi)ri{VrfZi) = C2 + 0( e -l^-^l/ fe ) (137) 
J fe=l 

N 

g m = / II d2 ^ # ife; (»/*.; zd = c + o( e -i'"^). (138) 
fe=i 

(We note that the number of screening charge coordinates w a is N + 1 for both N even and odd.) Thus, when there 
are only two quasiholes and an odd number of electrons, the effect of a counterclockwise exchange is the accrual of 
a statistical phase exp [in (^jy + |)1 and an Aharonov-Bohm phase exp (— i^jj^-)- This difference in the resulting 
phase is an indication of non-Abelian braiding statistics, specifically due to the fact that the two quasiholes must be 
in different fusion channels I and ip when N is even and odd, respectively. This is is discussed further in Section llXl 



VIII. PLASMA REPRESENTATION OF THE OVERLAP MATRIX FOR FOUR QUASIHOLES 

To represent all the entries of the overlap matrix as a plasma, we need to find Coulomb gas representations for 
arbitrary products of conformal blocks. While it does not seem possible to find such a representation for an arbitrary 
product TaTp (see Appedix [G] for an incomplete approach), it turns out to be possible to do so for particular linear 
combinations. These combinations, in turn, are nothing but the correlation functions of the order and disorder 
operators in the Ising model. The disorder operator (j,(r},fj) in the Ising model has the same scaling properties as 
the order operator (7(77,77), but it changes sign as it is taken around the order operator. The analog of Eq. (|1 10|1 for 
correlation functions of two disorder operators and two order operators is: 

(^1*^1)^2,772)^(773, 773)^(774,774)6(21,21) . ..e(z N , z N )) = Jb z l )T (f] ll ; z^ - Ji^; z l )T 1 (f] ll ; (139) 
(<r(m>m)v(V2,n 2 )a{ri3,fi3)v(v4,fj4)e(zi,zi) ■ -.e(z N ,z N )) = J^oiVn't z i)^i(n^ z~i) ~ Fiivtf Zi)^b(ty*;»i) 5 ( 14 0) 
( ^(771,771)^(772,772)^3, »73)m(?74, Vi) e(zi,z\) . . . e(z N , z N ) ) = Jb^; z i )J'i(fj IJ ,; Zi) + Ji^;^)^)^; Zi). (141) 

The expressions for these correlation functions without the energy operators e were given in Ref. [58l Since the 
transformation laws of the conformal blocks T a are the same as those of the wavefunctions ^ a (up to an irrelevant 
phase), we can use Eqs. (|30|) . (|3Tj) . (|32|) to verify that the expressions in Eqs. (|139j) . (|140j) . and (|14ip indeed change 
sign when an order operator is taken around a disorder operator. We can similarly verify that the right-hand-sides 
of Eqs. (I139p . (|140j) . and (|141[) remain invariant if an order operator is taken around the other order operator or a 
disorder operator is taken around the other disorder operator. 

If we can prove that the integrals fYlkd 2z k °f tne three expressions in Eqs. (|139l) . (|140|) . and (|141l) are equal to 
zero, then we will have proved that the overlap matrix defined in Eq. (Ill 6[) is proportional to the identity matrix, 
since we would know that 

Go,o — Gi,i = 0, Gq }1 — Gi Q = 0, G 0>1 + Gi Q = 0. (142) 

Combined with the already proven identity, ^ a G^ a = 2C 2 , this would prove that G^p = C 2 b a ^. 

To do this, let us construct the Coulomb gas representation for Eqs. (|139[) . (|140[) , and (|14ip . To the best of 
our knowledge, such a representation has not previously been constructed in the literature. It is straightforward to 
construct it, however, using what we have learned so far. First, let us consider Eq. (|139|) . We take Eqs. (|82jl. (|83|) for 
the holomorphic part of our representation [unlike Eq. (|115|) where we used the alternative Eqs. (|80]). (|81|)]: 

Mv„ zi) = { v™{ m )v™{ m )v™{m)v 2 °nm) t&'toOTte) • • • ^T(^-i)^T(^) ) , (143) 

7-1(77^; z^ = ( V^(m)V 2 \°( V2 )V^( V s)V^( m ) V 3 \°( Zl )V™(z 2 ) . . . V™(z N -i)V™(zn) ) • (144) 
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For the antiholomorphic part, we take Eq. (11431) but with the first two operators exchanged, which gives 

Jbfo; zi) = < v™{ni)v^{m)v™{m)v^{m) vl°(z 1 )v^(z 2 ) . . . ) ■ (145) 

For the antiholomorphic sector's other conformal block, we use this same string of vertex operators, except with the 
screening charge moved from the 773 vertex operator to the 772 vertex operator. This gives 

- Fifazi) = ( V^(fk)V^(fh)V^(fk)V^(fk) V£?{ti)Vg?{z2) . . . V™{z N . x )V™{z N ) ) , (146) 

so that combining these different representations gives the correlation function in Eq. (|139[) . 

In other words, we use a representation for the Ising disorder operators in which we use the vertex operators, not 
including possible screening charges, given by 

Mum) = V^(m)V^{fk) = e-*^**™-"^*™ = e^^Oifafk) (147) 

2V3g 

for the first disorder operator (at the position 771,771) and 

Km,m) = V™(v2)V™(fk) = e l ^ M - l ^^ 2) = j^^G^ (772,772) (148) 

2\/3g 

for the second disorder operator (at the position 772,772)- The two charges appearing in Eqs. (I14T[) and (|148[) in the 
holomorphic and antiholomorphic parts of the vertex operators are the two allowed charges which can represent the 
Ising spin, as can be seen from the table in Eq. (|72")l . Such objects (whose left and right charges are distinct from each 
other) can be viewed as particles in the 2D plasma which carry not only electric charge, but also magnetic charge, in 
the sense of Eqs. (fAij) . (fA7|) . ([A8]l in Appendix El 

Meanwhile, the Ising order operators without their screening charges are represented by purely electric operators 

-7(773,4,773,4) = V£W) V 3 T(»Jm) = e i ^^ +i ^^ (149) 
as are the Ising energy operators (without screening charges) 

e(zi, Zi) = V^(zi)V^(zi) = e ^M+^m) (150 ) 

and the screening operators 

e ia^ip(w) e ia-tp(w) _ e -i^ip(w)-i^ip(w) (151) 

In this construction we exclusively use the a_ = —y/3/2 screening operators. The crucial part of the proposed 
construction is that almost all of the operators used here are mutually local, i.e. are single-valued when any one is 
taken around any other. The exception is when an order operator is taken around a disorder operator (or vice versa), 
which results in a — 1 . This is easy to check if one uses 

Thus, when one is taken around the other, a phase 47r(a r /3 r — aifii) is acquired. For example, if a disorder operator 
represented either by a r — — , ai = or by a r — , ai — — is taken around the order operator represented 

by Pr = 4 j Pi = 4 i this phase is 7r. (Similarly, when a disorder operator is taken around an energy operator, 

represented by (3 r — (3i = it produces a phase of 2ir.) 

Of course, we would have obtained the same analytic continuation properties if we had switched the representations 
of the first and second disorder operators, i.e. had taken: 

Kvum) = v™( m )v™( m ) = e ^^)-^(m) (153) 

for the first disorder operator (at the position 771,771) and 

= V™{ m )V™{m) = e-'i^J+^to) (154) 

In fact, since this correlation function must be the same if we exchange the two disorder operators (or exchange the 
two order operators), we must take an equal linear combination of both representations for the disorder operators. 
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Thus, we conclude that we obtain the correct analytic continuation properties when operators are taken around other 
operators or when identical operators are exchanged if we write 

Fair) n~> z i) ^(v fi',Zi) - FiiVn] Zi)^i(fj^, z~i) = (^(771, 771)^(772, 772)^(773, 773)^(774, 774) e(zi,Zi) ■ ■ ■ e(z N ,z N )) 

= ( v^Mv^imW^imW^M v3?(*i)t&°(*a) ■ • • v™(zn-i)v™(z n ) ) 

x ( v^{n l )v^{n2)v^{nz)v^{f li )v^{z l )v^{-z 2 ) . . . v^-i)^ ^) ) 

+ < V^( m )V^( m )V^( m )V^( m ) Vg(zi)Vg(z*) ■ ■ • Vg(zN-i)V$?(z N ) ) 

x ( v^(m)v^(m)v^(fj 3 )v^(m)v^ (^)y 3 T(^) . . . v^-i^TM ) 
+ < v™(m)v^(m)v£(m)v^{m) vi?{z x )v™{z 2 ) . . . vi?{z N - X )v™{z N ) ) 

x ( ^(m)^Tfe)v 2 \^%)v 2 T(^)v 3 2 i (^)^T(^) . . • vi?(zN-i)v™(z N ) ) 
+ < v™(m)v£lm)v^(m)v^{m) vi?{z x )v™{z 2 ) . . . v£?{z N -i)v™{z N ) ) 

x ( Vg(fk)Vj?(fh)Vg(fk)Vg(fk)V£? (zi)V 3 °?(z 2 ) ■ ■ ■ V^(z N -i)Vl°(z N ) ) . (155) 

If we again pursue Mathur's strategy and replace the integrations over pairs of contours, such as § dw c § dw c , by 
integrations over the plane J d 2 w c , we obtain the following expression for Eq. (|155[) . 

Fo(Viji\ Zi)Fo(f)M, Zi) ~ Fiiv^, Zi)Ti{f)ij,\ Zi) = < m(?7i, 771)^(772, 772)^(773, 77 3 )cr(774, 774) e(zi, zi) . . . e(z N ,z N ) ) 

. N+l 4 

no —1 TT r — - - - - - - - — — 1 - 

d w c \r)i -772] 4 11 |_( r /i — ^m) 8 KVi ~ Vn) s (V2 - V^) 8 (V2 - %) 8 J|??3-??4| 4 
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AT+1 W+l iV TV 

x n ka-whi 3 n ni wa " z ^ 3 ni zj ~^i 3+ c - c - ( i5g ) 

a<6 a— 1 z— 1 i<j 

We emphasize that Mathur's procedure can be applied to this case because the screening charges have trivial mon- 
odromy with all other vertex operators, including, in particular, those of the disorder operators. Eq. (|156[) is one of 
the main results of this work. It is the correlation function of two order and two disorder operators in the Ising model 
obtained via the Coulomb gas approach. In order to understand why this expression is correct, it is helpful to observe, 
first of all, that the holomorphic and the antiholomorphic parts of Eq. (|156l) indeed reduce to the second through fifth 
lines of Eq. (|155[) while the holomorphic and the antiholomorphic parts of the complex conjugate [the "c.c." at the 
end of Eq. (|156[) ] reduce to the sixth through ninth lines of Eq. (|155[) . (The complex conjugate is also necessary to 
make the correlation function symmetric.) More importantly, one should note that the expression inside the integral 
in Eq. (I150[) is single valued if any of w a ,w a is taken around any other variable, or if any Zi,z~i is taken around any 
other variable, so that the integrals over these variables are well-defined. However, it changes sign if 771,771 is taken 
around 773,773, as well as if 771,771 is taken around 774,7/4, if 772,772 is taken around 773,773, or if 772,7/2 is taken around 
774,774. This is exactly as we would expect for the correlation function of Eq. (|139[) . Thus, when the d 2 w c integrals are 
decomposed into sums of products of conformal blocks, following Ref. 68 as outlined in Appendix [B] and the discussion 
following it, the analytic continuation properties automatically select the correct combination of conformal blocks. 

Similarly, the two correlation functions in Eqs. (|140|) and (|141[) can be obtained by a simple permutation of the 
variables 77^,77^. Thus, all three correlation functions from Eqs. (|139[) . (|140[) . and (|141l) can be constructed in this 
way. 

In fact, the analytic continuation properties noted previously are sufficient to conclude that if Eq. (|156|) is nonzero, 
it is equal to the correlation function we need to compute. Thus, if one were simply handed Eq. (|156[) . one could 
verify it without the arguments of Ref. |(38| (although, of course, one would probably not discover this equation without 
Ref. |68T ) by appealing to these analytic continuation properties and showing that the expression Eq. (|156l) is non-zero. 
In Appendix ICl we explicitly evaluate Eq. f|156[) in the absence of the energy operators. In this case, there is only one 
screening operator involved, and only one integral over w, which we calculate. The result of the evaluation, given in 
Eq. (|C21[) . explicitly produces the correct combination of blocks as given in Eq. (I156p . 

It further follows that when the energy operators are included, the representation Eq. (I156[) cannot simply vanish. 
Indeed, we can always take the four Ising order and disorder operators far away from the energy operators, and 
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the correlation function factorizes. It is clearly non-zero in this limit; by analyticity it will remain nonzero at finite 
separation between them. 

Thus, we have the following expression for the overlap integral of the difference of the product of conformal blocks, 
a generalization of Eq. (|117[) , 
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This particular expression gives us the integral over d 2 z k of the correlation function in Eq. (I139p . The correlation 
functions in Eqs. (|140l) and (| 141[) can be obtained by the simple permutations of the rj^ and 77^ coordinates. 

Now we reinterpret Eq. (jl5T[) as a partition function of a plasma by matching it against Eq. (|A4I) and rewriting it 

as 
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where T = 5, Q = an d m — !/%/%• We see that $2 is the 2D Coulomb interaction potential for TV particles 

with electric charge Q at positions Zi, N + l particles of electric charge — Q at positions w a , two particles of electric 
charge Q/2 at positions 773 and 774, a particle with electric charge Q/6 and magnetic charge to at position 771, and a 
particle with electric charge Q/6 and magnetic charge —to at position 772. Additionally, there is (implicitly) an electric 
charge of — Q/3 at infinity, corresponding to the background charge needed to maintain charge neutrality. Thus, we 
can interpret F2 as the free energy of a plasma at temperature T, comprised of TV particles of electric charge Q and 
N + l particles of electric charge — Q, in the presence of two test particles with electric charge Q/2 at positions 773 
and 774, and two "special" test particles which carry electric charge Q/6 and magnetic charge ±m at positions 771 and 
772, respectively (as well as a test particle of electric charge — Q/3 at infinity). 

We observe that this plasma screens, just as the plasma described in Eqs. (|118|) . (|119|) screened, due to a large 
number of ±Q electric charges at T > T C2 = Q 2 /4 present in it. However, when there are magnetic, in addition to 
electric, test charges that are at at sufficiently large distances from each other in the plasma, the free energy F2 will 
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diverge due to the confinement of magnetic charges, as explained in Appendix [D] and in Eq. (|Df 0[) . As a result, we 
find that Gq — l = (up to corrections that vanish exponentially as separations between r/^ become larger than 
the screening length), just as we expected. 

Constructing a similar plasma mapping for Gq^ ± Gf also results in a corresponding plasma that has test particles 
carrying magnetic charges. Thus, we can similarly show that G^ ± Gf — 0, and, consequently, we can conclude 
that the overlap integrals for the case of n — 4 satisfy 

G£„ = C 2 5 a , p + 0(e-l'"l/ £2 ). (160) 

This result is relevant to 2D chiralp-wave superconductors, whose real space wavefunction is the Pfaffian without the 
charge sector (Laughlin type factors)^, up to short-range modifications. Thus, Eq. (1 16 Oil presents another approach 



to computing the non-Abelian statistics of this state, distinct from that of Refs. 46 47 49.50.52. 



By a simple extension, we can now prove that the matrix of overlap integrals for the full wavefunctions are also 
proportional to the identity. Indeed, all we need to do is to multiply the integrand of Eq. (|157|) by the charge sector 
terms from Eq. (|12ip to obtain 
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with T = g, Q\ = \J2Mg, Q2 = \/3g, and 777.2 — l/\/3g. Again, F can be interpreted as the free energy of a 2D two- 
component plasma with two independent types of Coulomb interactions (denoted by subscripts 1 and 2) comprised of 
N particles (corresponding to the electrons) at Z{ which carry electric charges Qi and Q2, N + l particles (screening 
operators) at w a which carry electric charge —Q2, two fixed test charges (quasiholes) at 773 and 774 which carry electric 
charges Q\/2M and Q2/2, two fixed test charge (quasiholes) at 771 and 772 which carry electric charges Q\/2M and 
Q2/6 and magnetic charge ±7712, and a uniform neutralizing background of charge density p\ = — ymP (and p 2 = 0, 

but an electric charge of —Q2/3 at infinity). As previously mentioned, we expect such a plasma to be in the screening 
phase for T > T C1 ,T C2 , where T C1 — Q 2 /140 and T C2 — Q|/4 are the critical temperatures above which plasmas 1 
and 2 are individually in their screening fluid phase, so this plasma at temperature T = g with Q\ = ^/2Mg and 
Q2 = \/3g should be in the screening phase for M < 70. This has been numerically confirmed 57 for M = 2. When 
the plasma is in the screening phase, the free energy F will diverge due to confinement of magnetic charges, as long 
as the separation between the magnetic charges (at 771 and 772) are larger than the screening length £ of the plasma. 
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Thus, for sufficiently small M, we have now proved that Go,o — Gi.i = 0. By the suitable permutation of 77^, 77^,, we 
similarly obtain Go 1 ± G1.0 = 0. Combining this with the already proven relation Go,o + G1.1 = 2G, we have finally 
proven Eq. (j4Tj) 

G a ,p = C5 aJj + 0(e-^-^/ e ) (164) 

for the four-quasihole wavefunctions. 



IX. n-QUASIPARTICLE FUSION AND BRAIDING 

We would now like to go beyond the four quasihole case considered so far. At first glance, it might appear that 
going beyond four quasiholes is easy. All we need is to consider correlation functions of a larger number of order and 
disorder operators. Indeed, consider a wavefunction with N electrons and n quasiholes. As we know, there are 2 
such wavefunctions. Therefore, there are 22"" 1 • 2~~ 1 = 2™~ 2 overlap integrals [the number of entries in the G a ,/3 
overlap matrix of Eq. (|37[)]. On the other hand, we can imagine computing the correlation function of n Ising spins 
(as well as N Ising energy operators), as well as correlations with n — 2 order and 2 disorder operators, n — 4 order 
and 4 disorder operators and so on. If n/2 is odd, we should stop at § — 1 disorder and f + 1 order operators. If n/2 
is even, then we should stop at n/2 disorder and n/2 order operators. The total number of such correlation functions 
(which also depend on which operators are chosen to be order, and which are chosen to be disorder) is 
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if n/2 is even. Either way, the total number of combinations of conformal blocks one can get in this way (a generaliza- 
tion of Eq. (|139l) - (|14ip to n er-operators) is exactly equal to the number of entries in the G a ^ matrix, so computing 
the integrals over cPzk via the plasma construction would allow us to deduce what every entry of G a ^ is, for the 
general case of n quasiholes. 

Suppose we pursue the strategy of the previous sections to compute G Qj/ 3 for an arbitrary number of quasiholes. 
Eq. (|115p is easy to generalize to an arbitrary number n of quasiholes. Replacing the four quasihole operators in that 
expression with n of these operators resulted in Eqs. (|127|) . (I128|) . (I129[) . (I130p . which led us to conclude that 

TrG^ = qC 2 + Oie-^-nvUti) (167) 
TrG = qC + 0{e-^-^ /e ), (168) 

as discussed in Ref. l45l. Likewise, Eqs. (|157|) . (|158[) . (|159[) . (|16ip . (|162p . (|163p can also be generalized to n quasiholes 
by extending the product or sum over 3<A*<^<4to3</i<f<n and the products or sums over k = 3, 4 to 
k = 3, . . . , n. However, this only corresponds to the case of n — 2 order and 2 disorder operators. There are (™) such 
correlation functions in addition to the correlation function with n order operators. However, for n = 2 quasiholes, the 
correlation function with 2 order operators is equal to the correlation function with 2 disorder operators by Kramers- 
Wannier duality. Similarly, for n — 4 quasiholes, Kramers- Wannier duality reduces the number of distinct correlation 
functions with 2 order and 2 disorder operators from (2) = 6 to 3, as we have seen in the previous section. For higher 
n, Kramers- Wannier duality relates distinct correlation functions and therefore does not lead to any such reduction. 

For n — 2, we have a plasma representation for the only non-trivial correlation function and G a » is a 1 x 1 matrix, 
so we can compute it. For n = 4, we have a plasma representation for the 4 non-trivial correlation functions and G a ^ 
is a 2 x 2 matrix, so we can compute all of its entries. For n = 6, we have a plasma representation for the 1+ ( 2 ) = 16 
non-trivial correlation functions and G a ^ is a 4 x 4 matrix, so we can compute all of its entries. However, for n = 8, 
we have a plasma representation for the 1 + u) = 29 non-trivial correlation functions but G aj p is an 8 x 8 matrix, 
so we can compute fewer than half of its entries. The situation gets worse with increasing n since 2" -2 grows much 
faster than 1 + . 

In order to compute all of the entries in G a _p, we clearly need to be able to compute correlation functions in which 
there are more than two order and more than two disorder operators. Unfortunately, we do not know how to generalize 
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Eq. (|157[) to more than two Ising disorder operators, so we cannot compute by this strategy. The problem is 

that we use operators with opposite magnetic charges, Eqs. (|147[) and (|148p . to represent two disorder operators. If 
we have more than two disorder operators, then we need more than one operator of each type, Eq. (|147j) as well as 
(|148l) . However, Eq. (I147P is not local with respect to itself, and neither is Eq. (I148p . This prevents us from using 
Eq. (|147|) or (I148|) more than once in any correlation function. 

Another approach is to separate the screening operators into ones that are applied to tj) vertex operators and ones 
that are applied to a vertex operators, and then attempt to apply a procedure like Mathur's to change holomorphic- 
antiholomorphic pairs of tp screening operator contour integrals into 2D integrals over the complex plane. If this 
could be done, then the result before performing the a screening operator contour integrations would be a plasma 
with test particles. Among these test particles are ones corresponding to the holomorphic and antiholomorphic a 
vertex operators' screening operators, which, in addition to carrying electric charge — Q 2 /2 = —^3g/4, respectively 
carry magnetic charge =F3m2/2 = =fy/3/4g. Thus, if the holomorphic and antiholomorphic a screening operators 
are not paired up properly, they will give a vanishing result because of confinement of magnetic charge. This would 
give the sought after orthogonality (exponential in the separation of quasiparticles) for arbitrary n-quasiparticle 
wavefunctions, since the different conformal blocks (degenerate wavefunctions) correspond to assigning the screening 
operators to different a vertex operators. Unfortunately, there is a barrier in this approach to applying a Mathur-style 
argument, which is that there is a branch cut which prevents the J-terms from canceling in a simple way, and we 
have not managed to push the argument past this barrier. We provide more details on this incomplete approach in 
Appendix [Gj 

Fortunately, neither of these approaches is really needed. From the basic underlying structure of a topological phase, 
we know that, given a few basic assumptions (which rest on the assumption of an energy gap), we can deduce the 
braiding statistics of arbitrary numbers of quasiparticles, given much less information about quasiparticle statistics. 
This will be made precise in the following subsections. 

A. Braiding in the "Qubit Basis" 

In this section, we will use some features special to the MR Pfaffian state to deduce the braiding properties of an 
arbitrary number of quasiholes, given the Berry's matrices which can be computed from the 2-, 4-, and 6-quasihole 
wavefunctions by the methods described earlier. We will assume that the system is governed by a 3-body Hamiltonian 
with pinning potentials as in Eq. (TT9"]) . The only assumption we make is that there is a gap between the degenerate set 
of ground states with n quasiparticles at fixed positions (determined by the pinning potentials) and all higher excited 
states. So long as the gap remains open, the braiding properties that we discuss cannot change if the Hamiltonian is 
modified from Eq. (fl~9|) to something more realistic. 

From Ref. HH, we know that when there are n fundamental (charge e/2M) quasiholes at fixed positions, there is a 
degenerate set of states, rather than a unique state, and the following is a basis of these states: 



Zj ) M e- k ^ i? (169) 

where Ta =0,1. This double-counts the number of states, since these wavefunctions are invariant under the interchange 
of indices: (1+r - !, 3+r 2 , • • • , n— l+r a ) O (2— 7*1,4— r 2 ,...,n— r»). Thus, there are 2^ _1 such states. One can think of 
these states as each pair (1, 2), (3, 4), . . ., (n— 1, n) of quasiholes being a two-state system, i.e. a qubit, with an overall 
parity constraint on the n/2 qubits. For the sake of concreteness, we will fix this parity constraint by taking r» = 0. 
Note that this is a special feature of this particular topological phase; in a generic non-Abelian topological phase, the 
n-particle Hilbert space will not decompose into such a tensor product of two-state systems. Each of these two-state 
systems can be measured in the following way when the corresponding pair of quasiparticles is far from all of the others. 

Suppose you want to know if rj is or 1. Take 772.7-1 to 7/ n _i. If ^ (i+r 1 .3+r 2 ,....n-i)(2-r 1 A-r 2 n) now vanishes when 

any Zi approaches r?2j_i = Tfo-i, then rj = 1. If, instead, we take 7/ 2 j to r] n -i and ^(1+^,3+r ) ..., n -i)(2-r 11 4-r ,...,n) 
now vanishes when any Zj approaches 772 j = 7/ n _i, then rj =0. 

Consider, in this basis, the effect of a braid group generator T2i— 1, which executes a counterclockwise exchange 
of quasiholes 2i — 1 and 2i. This can be done with all of the other particles far away. Since the state in each of 
those other two-level systems can be measured while keeping those pairs far away (as described earlier), exchanging 
quasiparticle 2i — 1 and 2i must, by locality (which follows from the assumption of a gap), act as the identity within 
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each of those two-dimensional vector spaces. Therefore, it must be of the form: 

Tli-\ = &2 ® . . . ® 1 2 ® #2 



la 



Is 



(170) 



so that it only acts non-trivially on the i pair. Thus, the computation of r 2 j-i reduces to the computation of B?- 
As we discuss in the next subsection, the eigenvalues of B2 are the numbers i?j 2 1<Tl and R^" 1 ■ 

The braid group generator r 2 i affects pairs (2i — 1, 2i) and (2i + 1, 2i + 2) and must, therefore, take the form 
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so that it only acts non-trivially on the i th and (i + l) th pair. Once again, by locality, we can ignore all of the other 
quasiparticle pairs. 

Indeed, locality further guarantees that B% and B4 cannot depend on the number of other quasiparticle pairs in 
the system (since all of the other pairs can be taken far away), so long as there are enough so that there is, indeed, 
a two-state system on which B 2 can act and a four-state system on which B4 can act. If there is only a single pair 
of quasiholes, then there is a unique state, so at least four quasiparticles are needed in order to compute B%. The 
cognoscenti may object at this point by noting that there is a second state of two quasiholes, namely the state with 
an odd number of electrons. However, these states do not lie within the same Hilbert space, since they require their 
wavefunctions to have different electron number. In order to show that B2 can be computed from the combination 
of a N even electron number computation and a N odd electron number computation, we need to use more powerful 
consistency arguments, which are discussed within the next section. Put slightly differently, unless we know in advance 
that the MR Pfaffian state is a state of Ising-type anyons (i.e. unless we assume the answer), there is no reason to 
assert that the two eigenvalues of B% are given by the Berry's matrices of two quasiholes with an even or odd number 
of electrons. Fortunately, we do not need to make any such assumptions and can compute B2 from the four-quasihole 
case (or the six-quasihole case, though this is overkill). Similarly, B4 can be computed in a system with six quasiholes, 
for which there are the needed four degenerate states and which is the largest number of quasiholes for which our 
order-disorder operator strategy allows us to compute the full Berry's matrix. 

Thus, the representation matrices for all of the braid group generators can be obtained from the wavefunctions with 
four and six quasiholes and, therefore, the representation of the entire braid group can be obtained. (As mentioned 
earlier, if we use consistency conditions more effectively, we can reduce this to two- and four-quasihole wavefunctions, 
as discussed in the next section.) 

In order to actually compute the desired matrices B 2 and B4, we need to go into the qubit basis defined by the 
appropriate conformal blocks in the c = 1/2+1 theory, which are computed in Appendix IF1 
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On the left-hand-side of this expression, the indices take the values pi — 0, 1, respectively, and obey the overall parity 
constraint that X^P' be even. The preceding arguments hold for this basis as well, for which, given the plasma 
screening arguments for the four- and six-quasiparticle wavefunctions, the braiding matrices can be computed from 
analytic continuation to be: 
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Alternatively, if one used the standard basis, one would similarly have (for N even) 
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with 
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where again we can use the four-quasiparticle wavefunctions to compute B 2 and B' 2 (which are, respectively, i?( 1£ ^ 2 ) 
and B( 2 ^ 3 ) of Eq. (|3D|l). but must use the six-quasiparticle wavefunctions to compute B' A . 

Comparing the braiding generators in the qubit and standard bases, we see that t 2 .;_i acts on a two-dimensional 
subspace in the qubit basis, while T 2 i (as well as t\ and r„_i) acts on a two-dimensional subspace in the standard 
basis. Thus, we can imagine that if we had a way of ensuring a priori that braiding operations are consistent through 
changes of basis (which we can, in fact, show, but only after the six-quasiparticle computation of this subsection), 
then we would only need to compute braiding transformations on two-dimensional subspaces, and hence could avoid 
the need to compute anything using wavefunctions with more than four quasiparticles. 



B. General Considerations 



In this section, we use more powerful and general arguments to show that the braiding properties of anyons can 
be deduced from the two- and four- quasiparticle wavefunctions, provided we can compute the corresponding Berry's 
matrices. Although these arguments are more general in the sense that they do not use any special properties of the 
MR Pfaffian state, the assumed locality properties, which follow from the existence of an energy gap, are equivalent 
to those used in the previous argument. 

The long-distance, low-energy properties of quasiparticles (e.g. their braiding statistics) in a topological phase are 
assumed to be completely specified by an "anyon model," a.k.a. a unitary braided tensor categor y 14 ' 92 " — . An anyon 
model is characterized by: 

1. A finite set C of "topological charges" a, b, c, . . ., which are conserved quantum numbers specifying the different 
types of quasiparticle excitations. 

2. Fusion rules specifying how these topological charges may combine and split, as described by a commutative 
and associative fusion algebra 

ax6 = ^JV c 6 c, (179) 
cec 

where the integer N° b indicates the distinct number of ways the charges a and b can combine to produce charge 
c. For simplicity, we will restrict our attention to the case where there are no fusion multiplicities, i.e. N° b = 0, 1, 
since this is all that is needed for this paper, but there are more general cases 1 ^. There is a vector space V° & 
with dimV" & = N° b assigned to each fusion product, and one can represent the basis states of these spaces 
diagrammatically as 

b 

= \a,b;c) eVf. (180) 



The anyonic states describing a collection of quasiparticles can be represented by fusion diagrams such as 
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where the topological charge a,j of the jih quasiparticle is assigned to the jth endpoint at the top of the diagram. 
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3. Associativity of fusion within the fusion state spaces, which is specified by unitary isomorphisms written dia- 
grammatically as 

a b c a b c 

^^=E^ a6 1e/ (182) 
a d 

The i^-symbols are analogous to 6j-symbols, providing a change of basis between the basis in which topological 
charges a and b are first fused and then their result is fused with c to the basis in which topological charges b 
and c are first fused and then their result is fused with a. 

4. Braiding of topological charges enacting unitary transformations on the state space, which are written diagram- 
matically as 




(183) 



It is also worth mentioning that there is a unique "vacuum" charge, denoted or /, for which fusion and braiding 
is trivial. Furthermore, each topological charge a has a unique conjugate a with which it is allowed to fuse to vacuum 
(in a unique way), i.e. N° b = <5 ab . 

In summary, C, N£ b , [-P^ c ] e /' anc ^ comprise the basic data which completely specifies an anyon model. Given 
this basic data, one can describe the operation representing an arbitrary fusion and braiding process by using a series 
of applications of associativity and braiding (F and R). Hence, once this basic data of a system's anyon model is 
obtained, they can be used to compute the complete fusion and braiding statistics of (arbitrary configurations and 
exchanges of) its quasiparticles. 

In order to provide a coherent description of the fusion and braiding, this basic data must satisfy certain consistency 
relations known as the "polynomial equations" (i.e. the pentagon and hexagon equations)^, which ensure that any two 
series of applications of F and R starting in the same state space and ending in the same state space are equivalent^!. 
These consistency relations put strong constraints on the basic data (up to gauge transformations of the trivalent 
basis states). Consequently, it is often possible to start with an incomplete subset of the basic data and derive the 
rest using nothing more than these consistency relations. 

There are several important gauge invariant quantities worth describing here. The first is the quantum dimension 
of charge a, given by 

d a = d n =\[Ff aa ]J- 1 . (184) 
Quantum dimensions can be shown to satisfy the relation 

d a d b = Y. N abdc- (185) 
The second is the topological spin of charge a, given by 

9a = 9-a=T, = ^ [ ^ aa] 00 ( R a T ■ ( 186 ) 

cGC a 

These two (gauge invariant) quantities are particularly important, because it is often the case that the fusion algebra, 
together with the quantum dimensions and topological spins (or even a subset of them) can uniquely identify an anyon 
model. 

Once it is known that the braiding statistics of quasiparticles is given by analytic continuation of the wavefunctions 
given by CFT conformal blocks, we do not need to explicitly perform analytic continuation directly on the wavefunc- 
tions. Rather, one can instead obtain all the basic data directly from known properties of the CFT. For example, 
the topological spin is simply given by a = e l27Tha , where h a is the (holomorphic) conformal scaling dimension of the 
primary field corresponding to topological charge a, while the F- and i?-symbols can be obtained from the operator 
product expansions (the i?-symbols can simply be read off). For a complementary discussion of how the structure of 
topological phases can be decoded from wavefunctions - in particular, from their pattern of zeros - see Refs.[98||9^ 
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1. Isvng 

As an example of an application of this consistency, we consider starting with nothing more than the Ising fusion 
algebra 

IxI = I, I X ip = tp, I x a = a, ^ 187 ^ 
ip x ip — I , x a = a, a x a = I + 0. 

This fusion algebra is believed to describe the v = 1 bosonic MR state (i.e. M = 1, for which one has the identifications 
■02 = lo = I , &i = cr, and I2 = "00 = 0) and also the closely related SU(2)2 Chern-Simons theory (equating 1 = 0, 
cr = 5, and = 1). In the next subsection, we will show that this identification is correct. Given this fusion algebra, 
one can solve the consistency conditions to find exactly eight different anyon models (up to gauge transformations) 94 > 95 . 
These eight different anyon models are completely distinguished by their values of 6 a = e l27r ( 2 - J+1 )/ 16 , where j = 0, . . . , 7 
for the eight different theories [e.g. Ising anyons have j = and SU(2)2 has j = 1]. Hence, knowing only the fusion 
algebra one can use consistency to determine the theory up to this eight-fold degeneracy of theories. 

Since we only need to supplement the Ising fusion algebra with the value of 6 a in order to completely identify 
the anyon model describing such a system, we should examine it more closely to determine what is left to compute. 
However, we first note that one can easily determine the quantum dimensions to be di = = 1 and d a = v2 from 
Eq. (|185p . Now we can write out the definition 

B a = -L [R°° + R°°) = V2 [F™°] u i R T)* ■ (188) 



From this we see that all we need is R° a and either RZf or [F° aa ] 
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2. Putative Anyon Model for the Moore-Read Pfaffian State 



One can similarly analyze the anyon model which is generally assumed to correspond to the MR Pfaffian state. In 
the next subsection, we will compute the F and R matrices to show that this identification is correct. We define the 
topological charges through the corresponding CFT operators of the quasiparticles, which are the product of an Ising 
sector operator with a U(l) bosonic charge sector vertex operator 

I- = e l 75M v , 0j = vbe^ mV , Oj = aev^ v , (189) 

where j must be even for Ij and ybj, and odd for t?j. From the construction of wavefunctions from CFT correlation 
functions, we can identify the charge Iq as corresponding to vacuum (an insertion of this operator leaves the wave- 
functions unchanged), while 02M corresponds to the underlying particle of the system, e.g. the electron or atom. 
Consequently, when M is odd, we identify the bosonic underlying particle with vacuum, i.e. 02Af = ^o- When M 
is even, the underlying particle (the electron) is a fermion, so we cannot identify it with vacuum (though one can 
introduce a Z2 grading and put it in a doublet with vacuum). Instead, for M even we identify a pair of electrons with 
vacuum, i.e. 2 m x ip 2 M = hn = h- 

Furthermore, we can identify <j\ as the label corresponding to the charge e/2M fundamental quasihole. One can 
think of I2 as the e/M Laughlin quasihole obtained by inserting one flux. Alternatively, on can think of I2 and 0^ as 
the quasiparticles obtained at 771 when one takes r) 2 — > r\\ in the four quasihole wavefunctions and respectively. 
The labels /2fe+2, 02fe+2, and o~2k+i can then be thought of as describing quasiparticles obtained by similarly combining 
k charge e/M I2 Laughlin quasiholes with quasiparticles of type I2, 02, and a\, respectively. Thus, we have the fusion 
rules 

Ij Xl k = I[ J+ k] , Ij X fe = 4> []+k] , Ij X (T fe = a y+k] , , igQ . 
1pj X 0fe = I[j+k] , Ipj X (Tfe = CTy+k] j X (Tfe = iy +fc] + 0y +fe] . 

where we have defined the short-hand [j] = j'mod4M. For M odd, one has the additional identifications Ij = ipu+2M] > 
ipj = Iy+2M] i and aj = cry +2 M] • I n this way, there are 6M distinct topological charge types for M even and 
3M topological charge types for M odd. Given these fusion rules, the quantum dimensions can be found to be 
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1 (for j even) and d a , — v2 (for j odd) from Eq. f|185[) . The topological spins are 
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(191) 
(192) 
(193) 



As mentioned earlier, the eight consistent anyon models permitted for the M = 1 MR fusion algebra can be uniquely 
distinguish ed b y their value of 9 ai . For M = 2 (corresponding to the v = 1/2 fermionic MR state) , it was found 
m Refs. l95lfT00l that given the corresponding fusion algebra, there are 32 different possible anyon models that satisfy 
the consistency conditions and that these 32 different anyon models can be uniquely identified by their topological 
spins. In fact, closer inspection reveals that they can be uniquely identified merely by their values of 6< Jl and 9i 2 . 
Moreover, it can be showniffi that this is true for general M, i.e. that 9 ai and 9i 2 uniquely identifies the anyon models 
corresponding to the MR fusion algebra. 



C. Identifying the anyon model 



Once the anyon model corresponding to a state is known, one can use it to compute the braiding statistics of an 
arbitrary n-quasiparticle wavefunction. Thus, we now turn to the wavefunctions of the MR Pfaffian state, to extract 
the quantities necessary to identify the corresponding anyon model. 

We now give the charge e/2M quasiholes in this system the label o\. At the moment, this labeling is completely 
innocent, but we note that <j\ has the same quantum dimension as the similarly-named quasiparticle of the previous 
subsection since d ai = \[2 and, as discussed in Section IIX Al there are (v^)" states for n even quasiparticles. 
Now consider the wavefunctions vJ/ and "I^ defined in Eq. (2U These wavefunctions are linear combinations of the 
wavefunctions in Eq. [25] Suppose that, in fy , we bring the o\ quasiholes at 771 and 7y 2 close together and take the 
quasiholes at 773 and 774 far away. Then there is a localized excitation at 771 « 772. Let us call this localized excitation 
12- Suppose we do the same thing with <3/i; we call the resulting excitation ip2- A priori, we do not know if I2 and 
ip2 are distinct excitations or are topologically-equivalent (or, perhaps, are different superpositions of excitations). 
However, from the four-quasiparticle braiding matrices computed in Section [Villi through plasma analogy arguments, 
we know that if the quasihole at 773 is taken around the excitation at 771 « 772, the resulting phase in ^1 differs from 
that in by — 1. Thus, it is correct to assign two different labels to the corresponding excitations; the labels I2 and 
-02 are as good a choice as any. 

We now compute the F-symbols and the fusion algebra. This can be done with four-quasihole wavefunctions. 
The i^-symbol is a unitary transformation between two different bases of the two-dimensional Hilbert space of four 
quasiparticle states. One basis is given by ^0 and ^1, in which the ct\ quasiparticles at 77! and 772 fuse to I2 and i\)2-, 
respectively. From the computation of matrix elements in Section I Villi we see that these wavefunctions provide an 
orthonormal basis: 



<*o|*o) = = 1. 



(*o|*i) = (*i|*o> =0. 



(194) 



(Strictly speaking, we have shown that "Jo and vl/i are orthogonal and have the same norm, which is simply an 
77^ -independent constant in the limit that the 77^ are all far apart. Thus, to normalize them, we need to divide both 
wavefunctions by their common norm.) The second orthonormal basis is given by I^q ^ , in which the o\ quasiparticles 
at 771 and 774 fuse to I2, ip2- Since 1, 2, 3, and 4 are an arbitrary labeling of the quasiparticles, the states in which 1 
and 4 fuse to I21 1P2 are given by changing the labels to obtain: 
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Then, using the identity Eq. (|29|) . we have: 
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Therefore, 
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Thus, from Eq. (|197[) . we conclude that: 
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where a, b = 7a, ip2- One can similarly compute all the F-symbols directly from wavefunctions with no more than four 
quasiparticles. 

We note that we computed the F-symbols directly from the wavefunctions, with no appeal to orthogonality or 
the plasma analogy, so it might at first appear that orthogonality played no role here. However, it is important to 
establish that the associativity encoded in the F-symbols is unitary, since we are describing quantum mechanical 
systems. For this, the orthogonality result (obtained from the plasma analogy) is crucial, since it establishes that the 
wavefunctions in question provide orthonormal bases that are related by this f-symbol, which can thus be interpreted 
as a unitary change of basis transformation. 

We now compute the i?-symbols. With the definitions of I2 and "02 given earlier, we can read these off from the four- 
quasihole wavefunctions "Jo and "Ji, since our plasma analogy has shown that the braiding statistics is simply given by 
analytic continuation of the wavefunctions. Thus, as found in Eq. (|30l) for the counterclockwise exchange 771 ^772, we 
have i?j 2 1<Tl = e^VOT - ^ R^ 1 = e J7r (™"+§). We no t e that these i?-symbols can also be obtained from the analytic 
continuation of the wavefunctions ty' and in Eq. (I196P corresponding the counterclockwise exchange 772 ^> ?73- 
However, once we know that the braiding statistics is given by explicit analytic continuation of the wavefunctions, we 
can save ourselves the trouble of explicitly computing the i?-symbols from wavefunctions, since we know the analytic 
continuation of the wavefunctions is determined by the CFT through its operator product expansion. 

Generalizing this to quasiparticles of type <jj by appropriately modifying the U(l) charge sector vertex operators, 
we similarly obtain 
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Similarly, one can compute the remaining i?-symbols from analytic continuation by additionally introducing quasi- 
particles carrying topological charge Ij and/or tp^ into the wavefunction, since the screening properties can easily be 
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shown to still apply (as long as M is not too large). The results are 



R\ j * h = e m ™ (201) 
R+j** = -e™& (202) 
R T f" = = e ™& (203) 

R a 3 I k = R l k a 3 = / 204 x 

R ^k = jriw _ e w(A-4). (205) 

As previously mentioned, there is some ambiguity in these quantities, since they are not all gauge invariant. However, 
we can use them to obtain 

0j. = e' i27r OT (206) 
e. 4>] = e i2 *(&* + ?) (207) 

This can be used to specify the corresponding anyon model. In fact, as previously mentioned, only a subset of 
this information is absolutely necessary to uniquely identify the anyon model. Hence, the fusion and braiding of the 
MR state can be completely determined through consistency using only the fusion algebra and the plasma argument 
applied to the wavef unctions with at most four a quasiparticles. 

In fact, wavefunctions with at most two a quasiparticles are sufficient to compute the i?-symbols, provided that 
they are allowed to also include quasiparticles that can carry either of the Ising charges I or ip (which will not change 
any of the plasma arguments). In particular, consider the wavefunctions (for N even electrons) with two o\ quasiholes 
at 7/i and 772 and either quasiparticle of type 1% or ip2 at 773. These two wavefunctions can, respectively, be obtained 
from the four o\ quasihole wavefunctions ^0 and ^1 by taking 7/4 — > 7/3. Clearly, since 773 and 774 can be taken to be 
far away from 771 and 772, locality dictates that the braiding statistics factors R^ ai and R^f 1 obtained by exchanging 
771 772 in the wavefunctions with two o\ quasiholes will be the same as those obtained from the wavefunctions with 
four o\ quasiholes. 

We saw an effect similar to this at the end of Section lVlIl where we considered wavefunctions with two o\ quasiholes 
(and no other quasiparticles) when the number of electrons N was either even or odd. There we found the two different 
braiding factors R^ ai and RZ^ 1 for N even and odd, respectively. This is, of course, related to the same properties 
discussed here. In particular, each electron has topological charge if>2M, so changing between even and odd N for 
wavefunctions with two o\ quasiholes has a similar effect as changing between the type I2 and ip2 (at fixed N) of an 
additional quasiparticle, as described in the previous paragraph. However, the cases of N even and odd necessarily 
belong to different Hilbert spaces, since the number of electrons in the system is different. Hence, one could not have 
simply taken those results as proof of the non-Abelian statistics, since doing so would have required making additional 
assumptions about the nature of the state, essentially equivalent to assuming the answer. 



X. ORTHOGONALITY FOR UNMATCHED QUASIPARTICLES 



It is often assumed that the overlap between two wavefunctions that do not have the same types of quasiparticles 
at (nearly) the same locations should vanish. For example, such an orthogonality postulate is usecU°-i (sometimes 
implicitly) in the determination of the braiding statistics and other properties of hierarchical states, such as the Abelian 
Haldane-Halperin (HH) hierarchy state s 102 ' 103 , which are built on the v = 1/m Laughlin states, and the Bonderson- 
Slingerland (BS) hierarchy states^, which can be built on arbitrary states, notably providing a non-Abelian hierarchy 
over the v = 1/2 MR and anti-Pfaffian states. In fact, this orthogonality is of paramount importance, since it is 
necessary to establish the interpretation of the wavefunctions as describing anyonic quasiparticles. Without it, we 
would be missing the notion of a specific, distinguishable localized conserved quantum number (topological charge) 
associated with the quasiparticle coordinate, i.e. property 1 of anyon models in Section IPC Bl and thus the rest of 
the properties (fusion, braiding, etc.) that we can derive would lack a proper interpretation beyond some algebraic 
relations between certain special wavefunctions. In this section, we will prove this orthogonality using the plasma 
analogy. 

First we consider the wavefunctions in terms of their Coulomb gas CFT formalism. Next we recognize that, as 
described in Appendices [X] and |D] the chiral vertex operators can be expressed in terms of electric and magnetic 
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Coulomb charges. Specifically, 

e i av (z) = e »$#(*,s) ( Z; z) (209) 
corresponds to a particle carrying electric charge q = a/2 and magnetic charge m = a jig. Similarly, 

corresponds to a particle carrying electric charge q = a/2 and magnetic charge m = —a/2g. If the corre- 
sponding holomorphic and antiholomorphic vertex operators have coinciding positions, they produce an operator 
e iatp(z) e iaip(z) _ e ta<t>{z,z) correS p 0nc } m g to a particle carrying electric charge a and no magnetic charge. If quasiparti- 
cle coordinates are not matched up appropriately between two wavefunctions, then, in the plasma analogy, they will 
leave stray (uncanceled) magnetic charges in the plasma. Since magnetic charge is confined in screening plasmas, this 
makes it clear that the overlap between wavefunctions described by a plasma analogy will vanish unless they have 
matching quasiparticles at nearly coinciding positions. More specifically, the overlaps will be zero up to O fe - *"^) 
corrections, where r is the largest separation between unmatched quasiparticles (i.e. between magnetic charges in the 
plasma). 

With this orthogonality established, one might wonder whether a stronger result can be established. Indeed, one can 
attain stronger results for holes in the v = 1 filled Landau level and for single Laughlin-type quasihole wavefunctions. 
We explain this in more detail in Appendix [HI but note the results here that 

n n 

G(ffc,<)=Ci £ (-irJJc-i(l^l a +^l a -^M^)=C 1 (27r)" £ (-1)^ <*lll (w - <) (211) 

7TgS„ JU=1 7TGS„ f*=l 

for n holes in the v = 1 filled Landau level, and 

G (fj, rf) = C.e-M^+Wf^m') = Ci {2itM) ^ (j? _ ^ (212) 

for a Laughlin-type quasihole in any /y = 1/Af quantum Hall state with plasma analogy (for example, this holds for 
the I2 excitation of the MR state). 

Based on these examples, we conjecture that the overlap of wavefunctions that are not properly matched up in the 
U(l) boson charge sector will actually vanish with Gaussian, rather than exponential, falloff. This is expected to result 
from the neutralizing background that gives rise to the Gaussian terms in the wavefunctions. We also expect that 
the overlaps will behave effectively as delta-functions, projected into the appropriate subspace, while keeping track 
of the braiding statistics due to analytic continuation of the wavefunctions. This leads us to a general orthogonality 
postulate for quantum Hall states of the form 

G^(7K,<J~^(rw,<JE E IPM'M/O-^J* ( 213 ) 

aeC 7rGS„ a f*=l 

where there are n a quasiparticles of type (topological charge) a at coordinates r]^ a in one wavefunction and at 77^ 
in the other wavefunction (the delta- functions are only between quasiparticles of matching topological charge), and 
B a ,/3 {VnaTV'fia) ls a term that only keeps track of the braiding statistics due to exchange of the quasiparticles. For 
example, Laughlin-type quasiholes would have 

1 

n 

Eq. (1213[) clearly cannot be exact, since the right-hand side does not obey the necessary analytic properties mandated 
by the wavefunctions one is taking overlaps between, but we expect the answer with the correct analytic properties 
to have this effective form with respect to quasiparticle wavefunctions, up to exponentially suppressed corrections. 

XI. BRAIDING IN THE ANTI-PFAFFIAN AND BONDERSON-SLINGERLAND HIERARCHY STATES 



(214) 



In Refs. 3"> 30. it was pointed out that the particle-hole conjugate of the MR Pfaffian state is a distinct state and 
that, in the absence of Landau level mixing, it must be equal in energy to the MR state. Upon inclusion of Landau 
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level mixing, the anti-Pfaffian state appears to be lower in energ y 104 ' 105 . A candidate wavefunction for this state can 
be obtained by particle-hole conjugation 1 — of the MR wavefunction. 

In Ref.[6^, it was shown that one can construct hierarchical states over the MR and anti-Pfaffian v =1/2 states, and 
that these Bonderson-Slingerland (BS) hierarchy states provide candidates for all the (remaining) observed second 
Landau level plateaus. Moreover, it was numerically demonstrated^ - that the BS candidate for v = 12/5 is a 
competitive state, along with the RR and HH candidates. 

Here we will show that the wavefunction orthogonality results obtained in this paper for the MR wavefunctions 
imply the same orthogonality for the anti-Pfaffian and BS hierarchy wavefunctions obtained from them. Hence, the 
braiding statistics of these states are similarly obtained through analytic continuation of the wavefunctions. 

We begin by demonstrating that orthogonality of any two wavefunctions implies the orthogonality of their particle- 
hole conjugate wavefunction. For a general wavefunction ^{•q^z/) with quasiholes, one generates its particle-hole 
conjugate wavefunctioni ^ by filling one Landau level, introducing N\ holes at coordinates £ a , and projecting the 
holes onto the wavefunctions ^ by multiplying by ^{fj^] £ a ) and integrating over the holes' coordinates: 
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where v and S are respectively the filling fraction and shift of ^, while is the number of fluxes associated with 

the quasiholes of ^. 

Consider the overlap for two arbitrary particle-hole conjugate wavefunctions with quasiparticles (at matching po- 
sitions): 
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Now we use Eq. (|211l) to re- write this as 

JVj 
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~ JVi 



(219) 



where G l a iQ (?7 /J , is the overlap of the original two wavefunctions ^ p and ^ a (with N% electrons). 

It is now trivial to apply Eq. (|219l) to the anti-Pfaffian state, using our previously obtained results for G a ^ of the 

MR state. (For the MR state at v = 1/2, one has 5 = 3 and N^ h ' = n/2 for n fundamental charge e/4 non-Abelian 
quasiholes.) 

A similar, but slightly more complicated, argument applies to hierarchical wavefunctions, such as those of the 
Haldane-Halperin (HH) state d 102 ' 103 and Bonderson-Slingerland (BS) states^ 9 -. In particular, the wavefunctions for 
these states can be constructed by projecting Abelian quasiparticles of the Laughlin, MR, or anti-Pfaffian state into 
a Laughlin state, in a manner similar to how holes of the filled Landau level were projected onto a wavefunction 
to generate its particle-hole conjugate. One can then similarly use the results from Section [X] on orthogonality of 
wavefunctions with quasiparticles in different positions and the orthogonality postulate of Eq. (|213[) . 
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For concreteness, let us first examine the v = 2/7 HH hierarchy state. Wavefunctions of this state with n charge 
-e/7 quasiparticles at positions rj^ are given by 



n 

life 


n JVi 

- n,f 7 n n fe - 


U a ) Y\_ ( S 




fi—l a—l 


a<b 


N x 


Ni N 


N 


life 

a<b 


-« b ) i/3 nnfe- 

a—1 i—1 


Zi) Y[ (z% " 

i<j 



' e e " =1 



c=l 



(220) 
(221) 



where W i (tt a ; Zj) is the v = 1/3 Laughlin wavefunction with iVi = ^(A" — n) + 1 charge e/3 quasiholes at coordinates 
u a , which are projected into a Laughlin-type wavefunction 



n Ni 



$ fe; u a ) = n fe - ^) 3/7 n n fe - «■) n fe - ^) 7/3 12 R]Ua]2 

lt<v fi=l a=l a<b 

with n quasiholes at rj^. Taking the inner product and using Eq. (|213p . we obtain 
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where G*^^, ry M ) is the overlap for the Laughlin-type state's wavefunctions $, which, by the plasma analogy argument, 
is a constant, up to corrections that are exponentially suppressed in the separations — rj v \ between quasiholes. 
Hence, the inner product is equal to a constant, up to exponentially suppressed corrections. From this, it follows 
that the Berry's connection for transporting the quasiparticles at ?7 M is trivial (up to Aharonov-Bohm terms) and 
the braiding statistics of the HH hierarchy states are similarly obtained from direct analytic continuation of the 
wavefunction. The braiding statistics phase of e _l3Tr / 7 for a counterclockwise exchange of a pair of the charge —e/7 
quasiparticles in the v = 2/7 HH state can easily be read off the wavefunction in Eq. ([2201) . 

Now, let us examine the slightly more complicated case of the v = 2/5 BS hierarchy stated formed by condensing 
Laughlin-type quasiparticles (i.e. Ii excitations) of the v = 1/2 MR state. Wavefunctions of this state with n charge 
e/5 non-Abelian quasiholes at positions rj^ are given by 
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are the v = 1/2 MR wavefunctions with n charge e/4 o\ quasiholes at 77^ and Ni = + 1 charge e/2 72 Laughlin 
quasiholes at u a [ V I'^ IR are the MR wavefunctions with n charge e/4 quasiparticles defined previously in this 

paper] , and the N\ Laughlin quasiholes are projected into the Laughlin-type wavefunction 
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with n quasiholes at r/^. We take the inner product and use Eq. (|213|) to obtain 
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where G^(fj^, rj^) is the overlap for the Laughlin-type state's wavefunctions Hence, the inner product is equal, up to 
an overall constant and exponentially suppressed corrections, to that of the MR state times that of the Laughlin-type 
state (which is simply constant). From this, it follows that the Berry's connection is trivial (up to Aharonov-Bohm 
terms) and the braiding statistics of the BS hierarchy states are similarly obtained from direct analytic continuation 
of the wavefunctions. 



XII. STATES BASED ON OTHER CFTS 



The Coulomb gas formalis m 66 ! 67 that we used to describe the Ising CFT in this paper can similarly be used to 
describe any of the CFT minimal models M(p,p'), so we can use the methods developed in our paper to analyze the 
braiding statistics of candidate quantum Hall states constructed from such CFTs. Unfortunately, this becomes more 
complicated in general, since the structure of general minimal models is more complex than that of the Ising CFT 
(which is A4(4, 3)) and requires both types of screening operators to describe quasiparticles, which can create obstacles 
to applying Mathur's procedure 6 ^ or result in rather complicated plasmas about which little is known. However, we 
can make a few preliminary statements regarding certain cases. 



A. The M(5,4) state 

One may attempt to construct a quantum Hall state from the A^(5,4) CFT, using the e" = ^(4,1) = ^(1,3) operator 
with scaling dimension h e n —3/2 for the particles (electrons or bosons) in the same way the Ising ip operator was used 

for the MR stated. Specifically, the s"(zi) operator from the VW(5,4) sector is paired with an e*V M / 2 ^( Zi ) operator 
from a U(l) charge sector, and then the correlation function of an even number N of such operators produces a 
"ground-state" wavefunction at v = 1/M. (We note that M > 3 is necessary for the resulting wavefunction not to 
diverge as z, — > zj.) 

In the Coulomb gas formulation of the A^(5,4) CFT, the e" operator corresponds to vertex operators with 04,1 = 
— |a_ = or ai, 3 = — a + = ~\J\ (and appropriate screening operators). Using the ot\^ representation, we 
can apply the plasma analogy for this ground-state exactly as we did for the MR ground-state, resulting in a 2D 
two-component plasma of N charge Q2 = ~\f^g particles at Z{ (corresponding to the electrons/bosons) and N charge 
— Q2 particles at w a (corresponding to a+ screening operators) at temperature T = g. This plasma has coupling 
constant T 2 = Q\jT — 5, which is greater than the critical value T C2 = 4 above which such plasmas are not in the 
screening phase (as previously discussed). 
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Including the U(l) charge sector, we expect the resulting combined plasma to be screening in the U(l) charge 
sector for sufficiently small M, but not screening in the A4(5,4) sector. This can be used to demonstrate that the 
pair correlation functions of such states exhibit long-ranged correlations, and hence do not describe gapped states. 
The braiding statistics of quasiparticle excitations for such states is ill-dchncd, so whether or not we can construct an 
analogous plasma for wavefunctions with quasiparticles is irrelevant. 



B. The M{p,p- 1) states 

One may attempt to construct a quantum Hall state from the unitary M.(p,p — 1) CFT for p > 5, using the Abelian 
0(p-i,i) — 0(i,p— 2) operator with scaling dimension /i( p _x,i) = (p— 2){jp — 3)/4 for the particles (electrons or bosons) in 
the same way the Ising -0 operator was used for the MR state. In this case, this operator from the M(p, p — 1) sector 
is paired with an eW ' M / 2 f( z i) operator from a U(l) charge sector, where M is an even or odd integer depending on 
whether is integer or half-integer. The correlation function of an even number N of such operators produces 

a "ground-state" wavefunction at v = 1/M. 

In the Coulomb gas formulation of the M(p,p — 1) CFT, the (/>( p _i,i) operator corresponds to vertex operators 

with a p _i,i = — ^—^—O!— = ^r^y^p or a i,p-2 = — E -r La + = — ^y^^/p-zj (and appropriate screening operators). We 

can apply the plasma analogy for this ground-state similar to how we did for the MR ground-state, except in this 
case, the plasma analogy construction for the ground-state necessarily results in a two-component plasma where the 
different particle types carry charges of unequal magnitudes. Specifically, using the ai jP _2 representation results in 

a 2D two-component plasma of N charge Q2 — — \J ' ^pl y-^ particles at (corresponding to the electrons/bosons) 

and ^-^-N charge — ^z%Q2 particles at w a (corresponding to a + screening operators) at temperature T — g. 

Little appears to be known about the screening properties of plasmas with particles carrying charges of unequal 
magnitudes. Examining these plasmas coupling constants, the larger charges have Q\/T = p(p — 3) 2 /(p — 1) and 

the smaller charges have Q%/T = 4p/(p— 1). Since both of these are greater than T C2 = 4 (the critical 

value of the coupling constant for the two-component plasma with charges of equal magnitude) for p > 4, we expect 
(by comparison to the case with charges of equal magnitude) that these plasmas are not in their screening phase at 
temperature T = g. Including the U(l) charge sector, we again expect the resulting combined plasma to be screening 
in the U(l) charge sector for sufficiently small M, but not screening in the A4(5,4) sector, and consequently that 
these wavefunction do not describe gapped states. 



C. The Gaffnian state 



One may also attempt to construct a quantum Hall state from the non- unitary M(5, 3) CFT, using the ip = 4>(4,i) = 
0(1. 2) operator with scaling dimension h$ = 3/4 for the particles (electrons or bosons) in a similar way, producing 
the so-called Gaffnian wavefunctions 107 . In this case, the tp(zi) operator from the A4(5,3) sector is paired with an 
e iy/(2M+i)/4ip(zi) p era t or f rom a U(l) charge sector, and then the correlation function of an even number N of such 
operators produces the ground-state wavefunction at v = 2/(2M + 1). 

In the Coulomb gas formulation of the Ai(5, 3) CFT, the ip operator corresponds to vertex operators with 0:4,1 = 

— |a_ = \J~^ or ol\.2 — —\oi + — —\J~^2 ( an d appropriate screening operators). In this case, a similar plasma analogy 
construction for the ground-state results in a two-component plasma where the different particle types carry charges 
of unequal magnitudes. Specifically, using the a.4.1 representation for the ip operator results in a 2D two-component 
plasma of N charge Q2 = y / 27g/5 particles at Zi (corresponding to the electrons/bosons) and 37V/2 charge — 2Q 2 /3 
particles at w a (screening operators) at temperature T — g. Alternatively, using the representation for the ip 
operator results in a 2D two-component plasma of N charge Q2 = — y5g/3 particles at Zi (corresponding to the 
electrons/bosons) and N/2 charge — 2Q2 particles at w a (screening operators) at temperature T — g. We can also 
construct the analogous plasma for Gaffnian wavefunctions with up to two fundamental quasiholes. These quasiholes 

carry non-Abelian charge a, which corresponds to vertex operators with «2,i = = yjj m the Coulomb 

gas formulation. These map to test particles with charge Q2/2 in the analogous plasma resulting from the 0:4.1 
representation of ip- 

Examining these two plasmas coupling constants for their larger charges we have Q%/T = 27/5 and (2Q2) 2 /T = 
20/3, which are greater than T C2 — 4 (the critical value of the coupling constant for the two-component plasma with 
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charges of equal magnitude), whereas the smaller charges of these plasmas give (2Q2/3) 2 /T = 12/5 and Q\jT — 5/3, 
which are less than T C2 — 4. Hence, it is difficult to make a prediction by comparing to the two-component plasma 
with charges of equal magnitude. Whether or not these plasmas are in their screening phase is an interesting question 
that will be addressed in future research^. 

If these plasmas do not screen, then we expect that including the U(l) charge sector will result in a combined 
plasma that screens in the U(l) charge sector for sufficiently small M, but not in the A4(5,3) sector. This could 
then be used to demonstrate that pair correlation functions of such states have long-ranged correlations, and hence 
do not describe gapped states. However, if these plasmas are in the screening phase for parameters corresponding 
to the Gaffnian wave functions, then something must go awry elsewhere, as the braiding statistics resulting from 
analytic continuation of the wavefunctions with quasiparticles would result in a non-unitary anyon model, which 
cannot describe a topological phase. 



D. The Read-Rezayi States 



The next step in extending our work would be an adaptation of our plasma analogy construction to the Read-Rezayi 
Zfc-parafermion states^. To accomplish this, we need a Coulomb gas formulation of the Zfe-parafermion CFTs. Such 
a construction, based on the related Coulomb gas construction for the SXJ(N)k Wess-Zumino-Witten models, is known 
rather well^Sr^is, but it has not yet been developed to a point where an explicit representation of the overlaps of the 
desired wavefunctions (with excitations) can be written as a partition function of a 2D plasma. We suggest this as a 
subject for future work. 

XIII. DISCUSSION 

In this paper, we have constructed a new representation for the matrix elements between different MR Pfaffian 
four-quasihole and six-quasihole states. This representation allowed us to conclude that the Berry's matrix is trivial 
in the basis given by conformal blocks in ac = 1/2 + 1 CFT. This result implies that this CFT encapsulates the 
topological properties of this quantum Hall state - in other words, that this CFT can be used to compute braiding, 
fusion, etc. for quasiparticles in the MR Pfaffian stat o 34 ' 42 ! 43 and, by a straightforward extension, the anti-Pfaffian 
stat o 35 ' 36 and BS hierarchy states built on these^. This was, of course, the hope right from the beginning^ 4 -, but 
there was no proof, although there is strong evidence coming from a variety of arguments and numerical calculations, 
as we review below. Our paper provides a proof. 

We now consider the relation of our proof to previous results. In Ref. 43, a method was proposed to calculate the 
desired matrix elements, in terms of correlation functions in a perturbed CFT. The basic idea is that the integrals 
in Eq. (fTJ) can be written as the correlation functions of the operators corresponding to quasiholes if the operators 
which correspond to electrons are put into the action of the CFT and treated as a perturbation. For instance, TrG 
for n — 4 quasiholes, defined in Eq. (|126[) . is given by the product of the Ising model correlation function in Eq. 1)1070 
and the charge sector correlation function in Eq. (|12ip : 



r N 

TrG = / ]Jd 2 z k (^rhKOT^^Vfe^K^ 
^ fc=i 

x e ( Zl ,z 1 ) e 'Vf^^)... e ^ iZ - M )e'^»^)e- , i^^ 2z * (2 '^ (228) 

Here, we have made the abbreviation 4>(z, z) = (p(z) + (p(z). This can be re- written in the form 

r N a - 

TrG = / l[d 2 z k (0 9 (r ?1 ,r? 1 )0 9 (r ?2 ,?7 2 )0 g (r ? 3, V^O^, fj±) O^,^) . . .O e {z N ,z N )e~ l ^s S d2z ^ ] ) (229) 

^ k=l 

where 

O q (r], fj) = a(r), rfie'^k***",*) and O e (z, z) = e(z, s) e V¥*(*'*> (230) 
This can now be rewritten in the following form [Eq. (8.3) of Ref. |43] |: 

TrG = hm^(0,(r ?l! f? 1 )0 ? (77 2 ,f? 2 )0 9 (r ? 3 ! %)0 g (r 7 4,%) e*/^* 0.(^)^5^/^ 
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= limj^ {o q (rn,fh)o q {m, V2)o q (m, m)o q ( m ,f} 4 ) ^^ mJd2z4 ' izrz) ) c ^ c+to ( 231 ) 

This equation implies that we should view t J d 2 z O e {z, z) as a perturbation of the Lagrangian of the Ising model + 
a non-chiral boson and compute the correlation function of four O q operators in this perturbed theory. By taking 
N derivatives of the correlation function with respect to the coupling constant t for this perturbation, a correlation 
function with N electron operators is obtained and, therefore, matrix elements for iV-electron wavefunctions. The 
integrals in Eq. ([1} are exponentially-decaying at long distances if the renormalization group (RG) flow of the perturbed 
CFT is to a massive fixed point. Although this approach is highly suggestive, it is unclear how to show that the 
perturbed theory indeed flows to a massive fixed point. In Ref. 43, it was suggested that the perturbation has 
negative scaling dimension since correlation functions of O e increase with distance as a result of the background 
charge. However, it was pointed out in Ref. that O e should be viewed as an operator of positive scaling dimension 
(which is, in fact, an irrelevant perturbation for M > 2), with the background charge merely shifting the charge 
neutrality condition. Thus, it is not clear how to show that the perturbation in Eq. (|231j) is relevant. This is even a 
problem if this method is applied to the Laughlin states. Since all quantum Hall wavefunctions have a (Laughlin-like) 
bosonic charge sector, they all suffer the same complication. Furthermore, a leading-order RG calculation can show 
that the initial flow goes away from the conformal theory, but it is much harder to predict the fixed point at which 
the flow ends; it is necessary for the infrared fixed point to be massive in order for this approach to succeed. For these 
reasons, a Coulomb gas approach was initiated in Ref. Hi| this approach has been brought to fruition here. 

The method of Ref. |43j was recently discussed again by Read in an important paper—. It was shown there that 
this method can be applied more straightforwardly to the case of a chiral p-wave superconductor, which is described 
by the Ising CFT. This state has no charge sector and the perturbed Ising CFT is simply a free massive Majorana 
fermion. This result is important because it allows one to directly compute Rj a and , the two possible phases 
that result when two hc/2e vortices with fusion channel / and ip, respectively, are exchanged in a counter-clockwise 
fashion. This is a significant advance compared to approaches relying on the BCS wavefunction for this stat o 46 ' 47 i 49 ' 50 , 
which have only derived the ratio i?j <T /i?^ CT . Ref. [521 also presents an alternative calculation, based on a bosonization 
procedure related to the one which we use in Appendix [F] This calculation gives all of the needed matrix elements 
for two, four, and six quasiholes for a chiral p-wave superconductor (i.e. pure Ising CFT with no charge sector). 

A recent Monte Carlo evaluation^ of the matrix elements between MR Pfaffian two-quasihole states with even- 
and odd- numbers of electrons is consistent with R^ ai = 1 (which was already obtained in an earlier numerical 
calculation 4 ^) and R^" 1 — i for relatively large system sizes N ~ 150. Numerical diagonalization^ 4 of the 3-body 
Hamiltonian in the presence of pinning potentials for the quasiparticles is also consistent, for small systems (N ~ 16), 
with these R- matrices and also with the fusion rules of the c = 1/2 + 1 CFT. (The calculations of Refs. [H and are 
very similar, in principle. The difference - aside from the difference in system sizes - is that wavefunction overlaps are 
computed by Monte Carlo evaluation of overlap integrals in Ref. [H. In Ref. H3, different wavefunctions are computed 
in an orbital basis by exact diagonalization of the Hamiltonian with pinning potentials; the overlaps are then obtained 
from the inner products of the corresponding vectors in this basis.) 

It should be straightforward to extend the numerical calculations of Refs. 48.53.54 to compute the F-matrices as 
well and, thereby, fully determine the braiding properties of quasiholes using the logic of Section IIXB I It should also 
be possible to extend those numerical calculations to the six-quasihole case and, thereby, fully determine braiding 
without making any assumptions beyond the gap. 

Finally, quasiparticle braiding has recently been computed^ 1 - using coherent states in bases obtained from the 
"thin-torus" quasi-one-dimensional limit of the v = 1 bosonic version of the MR Pfaffian state. The key step in this 
derivation is changing between the basis obtained from the limit in which the torus is thin in one direction to the 
dual basis obtained from the limit in which it is thin in the other direction. The relation between these two bases 
(which one might recognize as the modular S-duality) is constrained by their properties under magnetic translation. 
The change of basis is not computed directed, but rather is determined by consistency. However, this only determines 
the change of basis and subsequent braiding relations up to an 8-fold degeneracy. This is precisely the same 8-fold 
degeneracy of anyon models that are consistent with the Ising fusion algebra (i.e. are solutions to the pentagon and 
hexagon equations)^, as described in Section IIXB I Thus, the results obtained by this method are equivalent to 
assuming no more than the fusion algebra (or, equivalently, the modular S-matrixa2£) and locality. 

It should be noted that an alternative approach to quasiparticle braiding in this state relies on a mapping between 
the MR Pfaffian state and a chiral p-wave superconductor. Wavefunction analytic continuation has been computed in 
the latter state directly from the BCS wavefunctio n 46 i 47 i 49 i 50 , and the Berry's matrix has been computed^ up to an 
overall phase (these studies find the ratio RT T /Rj' 7 — i but cannot determine itself; in this respect, our calculation 
gives more information, as does the calculation of Ref. l52l) . However, the mapping between the MR Pfaffian state and 
a chiral p-wave superconductor— is not exact. While the approximate mapping between the two is highly suggestive, 
the relationship between the two states is, strictly speaking, established by computing quasiparticle braiding in both 
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and comparing the result (i.e. comparing universal quantities in both states). This is done by comparing the results of 
this paper with the combined results of Refs. 46.47.49.50. Finally, the statistics of vortices in the non-Abelian phase 
of Kitaev's honeycomb lattice model^, which exhibits Ising topological order, has recently been computed through 
an explicit numerical computation of the Berry's matrix>ii£. 
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Appendix A: The free boson and the Coulomb gas 

This section sets the conventions and the normalizations for free bosons used throughout this paper. 
Consider a free boson field (compact such that (j> = <j> + 2ir) with the action 

S= JL [ d 2 z{Vcpf . (Al) 

47T J 

The value of g is often not particularly important. In the CFT Coulomb gas formalism convention that we follow, we 
set 5 = 1/4. The correlation function of <p is given by, up to an additive constant 

(4>{z,z )<f>(w,w)) = ! l -. (A2) 

9 

The electric operators are given by e i9fc< ^ Zfc:Zfc ) . Their correlation function is 

\zk~zA, (A3) 



h I V 9 k<l 



where one must have charge neutrality ^ fc = 0. 

We can also consider magnetic operators. They can be thought of as vortices in the field cj>, such that <f> changes 
by 2-Km when going around a magnetic charge m. Let us denote magnetic operators by O m . For electric operators 
located at positions z k and magnetic operators at positions w a , the correlation function isii£ 



A- 



exp j - flfcg; In 
\ g k<l 



W ~ zi\ + g ^2m a m b In \w a - w b \ + q k m a sag (z k - wi) \ , (A4) 



a<b k,a 



where J^k Ik = J2a m <> = °- 

Throughout the paper, we also use the holomorphic field f>{z), such that 

(<p(zMw)) = - ln( 7 W) - (A5) 
We can write 4>{z, z) = tp(z) + <p(z) and treat <fr(z, z) as the sum of two independent fields ip(z) and <p(z). As a result, 

(]le ia '- vM )=]l(z k -zi) S ^ L . (A6) 

\ k I k<l 
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We notice that, just as exp [iq<j>(z, z)\ = cxp [iq (<p(z) + (p(z))], we also have O m (z, z) = exp [img ((p(z) — <p{z))]. Hence, 
we can write 

eiav {z) = e i$<K*,*)o^(z,z) (A7) 

e ia<p(z) = e i%<K>,*)0-jL(z,z). (A8) 

In other words, the holomorphic/antiholomorphic vertex operator with coefficient a corresponds to an operator 
carrying electric charge q = a/2 and magnetic charge m = ±a/2g. 

For further reference, we also show how to calculate the important correlation function 

/nc fa ^)e^/'^w\=n(* fc -^)^e*? afcP/ ^ tog( **-* ) l (A9) 

\ k I k<l 

where now we can have a^ 7^ 0, since there is a background charge density p that can be used to maintain 
charge neutrality through p J d 2 z = — J^k a k- The integral over the logarithm can be calculated under the following 
assumptions. First of all, we take the domain over which z is integrated to be a disk. This means the imaginary part 
of the logarithm integrates to zero for symmetry reasons. As for its real part, we observe that 

V 2 log|z| = 2tt5 2 {z,z), (A10) 

where V 2 is the Laplacian and 5 2 (z,z) is the 2D delta function. This allows us to calculate this integral by solving 
the corresponding Laplace equation 

V 2 W J d 2 z\og\w- z\ =2tt, (All) 

which gives 

J d 2 z\og\w - z\ = ^M 2 - (A12) 

Taken together, this gives 

/ J] e ia^(*») e i, J \ = Y[( Zk Zl )^^ akW? . (A13) 

\ k I k<l 

For the charge sector of quantum Hall states, we had electrons with a = Qi = \J2Mg and a neutralizing background 
charge density p\ = — ^7 that gives rise to the Gaussian factors. 



Appendix B: Mathur's Procedure for Relating Products of Contour Integrals to 2D Integrals 

The purpose of this appendix is to review Mathur's trick which expresses sums of products of conformal blocks in 
terms of a 2D integral (i.e. as the classical Boltzmann weight for a plasma). Consider a 2D integral of the form 

/ <PwJ2f a (w)Q af ,ff,{w) (Bl) 

We will assume that f a {w)Q a pfp(w) is single- valued in D so that this integral is well-defined. Let us suppose that 
D is a simply-connected region with no singularities. Then, we can re-write the integral over the interior of D in the 
form 

J D d 2 wf a (w)Q aP .Uw) = Jj 2 w j a {w)Q a p^- Qf dw'ffiiw'^j (B2) 

where P is any point in the interior of D. Since there are no singularities in D, the integral from P to w is independent 
of the path. Then, since f a (w) depends only on w and not w, 



f d 2 w f a (w)Q a pf {w) = J d2w '^( K J dw ' ' fa{w)Q a pfp{w') 
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- f dw dw' f a {w)Q a pfp(w') 

1 J dD JP 

% - dw J a {w)Q a[3 (J dw'fp(w') ) ( B8 ) 



The penultimate step involves an integration by parts. Since x — (w + w)/2, y = (w — w)/2i, the complex derivative 
is defined by dg = dg/dw = (d x g — id y g)/2. Thus, 

/ d 2 wdg— / dxdy (d x g — id y g)/2 — — (p dw g (B4) 
J d J d 2 J dD 



We are interested in integrals of the form Eq. (|B1|) in which f a (w), f a (w) are conformal blocks. Thus, we expect 
them to depend on the coordinates z\, Z2, ■ ■ ■ , z m of all of the other fields besides the one at w, and it will be singular 
when w approaches any zu- In the expressions that concern us, w is the coordinate of a screening charge and the Zk are 
the coordinates of the other screening charges, the electrons, and the quasiholes. In order to avoid these singularities, 
we split the complex plane into annuli with inner and outer radii 1| and \z^\. We additionally define the 
points zq = and z m+ i = oo, so that the annuli cover the entire complex plane. Each of these annuli contains no 
singularities. However, they are not simply-connected. Therefore, we cut the annulus Ak open along a line from Zk-i 
and Zk , as shown in Fig. [2] We then have a simply-connected region bounded by the union of the curves Ck, Sk, 
— Ck-i, —S' k - The circular contour — Ck-i runs from P% to P£, while Ck runs from to Pi. Then according to 
Eq. (fim 

d 2 w f a (w)Q al 3f(s(w) = - f dwf a (w)Q a fi dw' fp(w') 

l JdA k J Pi 

- I dwf a (w)Q a fi dw'fp(w') + - dwf a (w)Q a p dw' ' fp(w') 
2 Jc k J Pi 2 Js k J Pi 

if _-_ f w i /if - - - f w i i 

dw f a (w)Q a p dw'fp(w) + - / dw f a (w)Q a p dw fp(w') (B5) 



Let us now define 



2 >-C k -i J Pi 2 JS' JP, 



pP' pit) 

J c = V / dwf a (w)Q a p dw'fp(w') (B6) 
(Ic)a = j p dwf a (w) (B7) 



49 



where C is a contour from P to P' Then, we can re-write the four terms on the right-hand side of the second equality 
in Eq. (|B5p in the form: 



y~] / dw f a (w)Q a p l dw'fp(w') = J Ck 

a,P J °" Jp l 

/ dw j a {w)Q a fi / dw'fp(w') = \^ dwf a (w)Q a pi / dw'fp(w') + / dw'fp(w') 
a Js k J Pi ^- / Js k \Jpi J Pi 



(B8) 



a,0 ' 

J2( ! s k ) a QMlc k )p + J Sk (B9) 



V* / dw f a (w)Q a p I dw'fp(w') = - V" / dw f a (w)Q a p ( / dw'fp(w') + / dw' fp(w') ) 

= -E( J ^-x)«^( 7 ^)^--/c fc _ a (BIO) 



a,j9 



dwf a (w)Q a J dw'fpiw') = -J2 dwf a (w)Q a J dw'ff 3 (w') = -J s > (Bll) 
-si Jp,; „ a J St J Pi 



a, p " °k 



Thus, we have: 



a,0 Ak ~ " a,p 



SUP 



When wc sum over the different annuli Ak, the Jc k — Jc k - ± terms will cancel. Now consider the terms in the second 
set of square brackets on the right-hand-side of Eq. (|B12|) . If f a (w) is taken in a counterclockwise direction from a 
point w on S k to the corresponding point on Sk, then it is transformed by the monodromy matrix M: 

[/"HU = X> Q/3 [/„(«,)]«,, (B13) 

P 

This monodromy is unitary (at least for the unitary CFTs), so M _1 = M'. Since f a {w)Q a pfp(w) is single-valued, 
M satisfies 

ApQM = Q. (B14) 

Here, we have used matrix notation and suppressed the indices. Consequently, Js k — Js' ■ This leaves only the terms 
in the third set of square brackets on the right-hand-side of Eq. (IB12[) . to which we now turn. 
Since there are no singularities in Ak, 

(I Ck ) a + (I Sk )a ~ (Ic^Ja ~ (^)a = (B15) 

Meanwhile, 

(Is k )a = ^2M a p(Ig it )p (B16) 

P 

Combining these two equations, we have 

I s , = (l-M)- 1 (/ Cfc -/o Jl _ 1 ) 

I Sk = {M^ -l)-\l Ck -I Ck ^) (B17) 

Substituting these expressions into the third set of square brackets on the right-hand side of Eq. (|B12|) . we find that 
the cross terms cancel so that we are left with 



d 2 wy2fa( w )Qapfp(w) = / d 2 W f a (w)Q a/3 ff j (w) 

J a,P k,a,P^ Ak 

= i-lc k QQ- - Mu)- 1 ^,, + Icu-Ml - Mfc)- 1 /^.,] (B18) 

k 

= Hc k Q(l - Mk)- 1 ^ + IcuQO- ~ Mk+^Ic,] ■ (B19) 

k 
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We have added a subscript k to M in this equation to emphasize that this is the monodromy matrix which results 
from deforming' S' k to Sk, and regrouped terms to obtain the form in the last line. This expression is Eqs. (2.4) and 
(2.19) of Ref. Ifia, where there is a slight typo in Mathur's Eq. (2.19), in that the index i—1 there should actually be 
i + Note that, by construction, both the left- and right-hand sides of this expression are single- valued. By applying 
Eqs. (|B18I) . (|B19[) repeatedly to each screening charge integral, we obtain the desired sum of products of contour 
integrals. 

We now specialize to the case of the Ising CFT in order to demonstrate some of the potential technical subtleties of 
applying Mathur's procedure. Consider the squared norm of the correlation function of two ip fields (or, equivalently, 
the correlation function of two energy operators). The generalization to N tp fields (where N is even) gives the square 
of the Pfaffian, which is a factor in the square of the ground state wavefunction. The basic structure is already 
apparent for just two tp fields, however, so we will begin with this simple case. Let us put the two tp operators at 0, z. 
We will call the screening operator coordinates Wi, w 2 . There is a single conformal block so, with appropriate choice 
of normalization, we can simply take Q = 1. When a screening operator e la - ip is taken around a i/)'b vertex operator 
or around another screening operator, it changes by —1. Thus, the monodromy Mk is ±1, depending 



l<X3l<p — „-ia_ip 



on whether an even or odd number of operators is contained within Ck- 
two tp fields 



Consider the left-hand side of Eq. (|9Tj> for 



d wi d W2 



— ia_ifi[z) ia_<p(w 1 ) e ia_ip(w 2 ) g — ia_y(Q) 



r 2 

/ d 2 wi d 2 w 2 \z\ 3 \wi - w 2 \ 3 Y[ (l 



z — wA 3 \w,\ 3 



(B20) 

We can show that it is equal to the right-hand side of Eq. (p?T]) by applying Eq. (|B19|) . First, we use Eq. (|B19[) to 
reduce the w 2 integral for fixed w\. There are two cases, < \w\\ < \z\ and \wi\ > \z\. For \w\\ < \z\, Eq. (|B19I) tells 
us that: 



d w 2 



e -ia_ip(z) e ia_tp(wx) e iu_ip(w 2 ) e ~ia_ip(0) 



\ Wl \<\z\ 



[Q--M z )- 1 -(l-M 1 )- 1 ]l Cwi Ic vl 

I [(1 - A'^)- 1 - (1 - Af,)- 1 ] I c Jc z (B21) 



Here, C Wl is the circle \w\ = \wi\ — e and C z is the circle |io| = \z\ — e. M\ is the monodromy of the holomorphic 
conformal block 



_¥>(*) 



(B22) 



when w 2 encircles the origin with < |iOa| < |wi | ■ Since such a circle encloses the ip field at the origin's vertex 
operator, e -1 "-^ -*, we have Mi = —1. M z is the monodromy when w 2 encircles the origin with \wi\ < \w 2 \ < \z\. 
Since such a circle encloses both the tp field at the origin's vertex operator, e _lQ -^( ) j and the screening operator 
e ia_tp(wi) ^ ^.jjjg m0 nodromy is M z — 1. Finally, is the monodromy when w 2 encircles the origin with \w 2 \ > \z\; 
since such a circle encloses the the tp fields' vertex operators, e~ la - v ^ and e~ la - v ^ z \ and the screening operator 
e la ~^ w ^\ this monodromy is I m = -1. 

Clearly, there will be divergent terms due to M z = 1. However, these divergence terms cancel, and so can be 
regulated (e.g. by adding small values to the vertex operators' charges that are taken to zero at the end). We will do 



this by letting q 
have 



T and q = 1 — S where S will be taken to zero at the end. Hence, for \wi\ < \z\, we 



d w 2 



e ~ia_ip(z) e ia_tp(wx) e ia_ip(w 2 ) e -ia_ip(0) X j 



\wi\<\z\ 



[S- 1 - (I - q)- 1 ] 




dwo e~ ia - ip ( z ' ) e ia - i p( w ^) e K*_v(w>2) e - ia -v(o) 



dw 2 e~ ia - v(z ^ e ia -V( w 2) e ia_<p(w!) e -ia_ v (0) 



(B23) 
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For \wi\ > \z\, we obtain a similar expression, but with M z = = — 1 and Mi = 1 — 5: 

2 

K|>|*| 



d 2 w 2 



e —ia_tp(z) e ia_<p(wi) e %a_<p(w2) e ~ia_tp(0) 



= ~ 2 [(1 - 9 2 )" 1 - (1 - l)- 1 ] 



+^[(l-g)- 1 -(l- 9 2 )- 1 ] 



dwo e ia - v( - Wl ^ e~ ia - v( - z ^ e ia -f( w 2) e -ia_(p(0) 



C z 



dw 2 e ia - v< - Wl ^ e ia -<P( W2 ) e - ia -v( z ) e - ia -v{o) 



We can now perform the w\ integral. For instance, 



(B24) 



d 2 w\ 



\wi\<\z\ 



dwo e -iQ - v ( z ) e ia ~f{t"i) e ia_ip(w 2 ) e -iu_ip(0) 



dwi f dw 2 e- ia - v(z) e la - v(wi) e la - vi - W2) e~ la - v{0) 
c z Jc z 



J-terms. 



(B25) 



Since w 2 is always enclosed when w\ encircles the origin (by definition, C\ is the circle at radius \w\ \ — e), the winding 
of w\ around w 2 contributes to the monodromy, in addition to the winding of wi around and w 2 around 0, giving the 
total combined monodromy of M = q. When w\ lies on the contour C z , the contour C Wl becomes the same contour 
(but point split, so that it is at infinitesimally smaller radius.) The J-terms cancel off, and so can be neglected. 
Similarly, we get 



/ d 2 wi ( f L 

J\wi\>\z\ \JC m 



dw 2 e ia - v ^ Wl ^ e iQ -^(™ 2 ) e ~ ia -v(2) e -»«_¥>(o) 



d 2 w\ 



KI<M 



d Wl f dw 2 e- ta - v(z) e m - viwi) e m - viw2) e-™-^ 
c z Jc z 

2 



J-terms 



(B26) 



dw 2 e~ ia - V< ' Z ' > e ia -V( w i) e ia_ip(w 1 ) e -iQ_ip(0) 



c z 



^[-(l-ff)- 1 ] 



d 2 w\ 



\wi\>\z\ 



dwi f dw 2 e-'"-*' e^-v^i) e ^-H>{w 3 ) e -««_<p(o) 
c z Jc z 

2 



+ J-terms (B27) 



dw 2 e ia -V{ w i) e -ia_ip(z) e ia_ip(w 2 ) e ~ia_tf>(0) 



Putting this all together, we find 
d 2 W\ d 2 w 2 



dwi f dw 2 e ~ ta - v{z) e ta - v{wi) e la - v{w2) e - ra -^°) 
c z Jc z 



J-terms. 



(B28) 



e -ia_tp(z) e ia_tp(wi) e ia_<p(w 2 ) e ~ia_ip(a) 



dw! f dw 2 e~ la - v(z) e la ~ v{wi) e la - v(w2) e- 101 -^ 
c z Jc z 



1 

Tg 



(v 3 t>)v&°(o)>| 



(B29) 
(B30) 



where we used q = — 1 and the definition of screened vertex operators to obtain the last line (notice that the divergent 
terms canceled prior to taking 6 — > 0). This is precisely the two-?/> version of Eq. (1^7)) . We note that we could have 
derived the same result with the "0 fields at two arbitrary points z\ and z 2 . In that case, one would produce screening 
charge contours at both radii \z\\ and \z 2 \ using this procedure. In the end, the terms with screening charge contours 
at the smaller radius will vanish due to Felder's rules, so one will again be left with both screening charges attached 
to the ip field's vertex operator that is at the larger radius. The derivation here allowed us to avoid some extra steps 
by situating one of the ip fields at the origin. 



r>2 



From the preceding derivation, we can now see how the general N-ip case works. Since the monodromy for taking a 
screening operator around a ip operator is —1, there will only be a non-trivial contribution from the contour associated 
with every second tp operator. Consequently, as in Eq. (|B29[) . there will be two screening operator contour integrals 
attached to every second ip operator, precisely as on the right-hand side of Eq. ([97)1 . 

Turning to conformal blocks with a fields, we consider first the combination that gives us the trace of the overlap 
matrix, namely the correlation function of 17(77, fj) operators and N energy operators. For this combination, Q a p — 5 a p. 
The monodromy matrices are again diagonal: the monodromy for a screening charge to go around a ip operator is 
— 1; to go around a a operator, it is ±i. In a similar manner to the steps that led to Eq. (|B29|) . the diagonality of 
Q and M simplifies matters and leads to Eq. (|115[) . Correlation functions with both order and disorder operators, 
which give us the off-diagonal elements and the difference between the diagonal elements of the overlap matrix, are a 
little more complicated because Q is no longer diagonal. They are considered in detail in the next appendix. 



Appendix C: Correlation function of two order and two disorder operators in the Ising model 

As explained in Section IViTll we can take advantage of Eqs. ([T4"g|) , (fT4"4")) . (|TI5|) . (|H5]I . (fH7|) . (HE]), and (H"4Uf to 

represent the correlation function of two order and two disorder operators in the Ising model in terms of the following 
Coulomb gas correlator 

{^{m,m)v{mJn)p{mim)o{m,m)) = (ci) 

Evaluating the correlation function results in the following expression, which is the particular case of Eq. (I156P with 
N = 0: 

_ _ _ _ — — _ 3 3 _ _ 1 _1 3 

(^1,^1)^2, m)cr(ri3,m)<j(m,v4)) = {visvu) 8 (rinvu) 8 (V23V24) 8 (77237724) 8 m -ml 4 \m~m\ 4 

/f 1 3 3 3 "l 

d 2 w - 771) (w - f\2)Y \{w-r) 2 ){w-f)i)}~ 1 - \w-rj 3 \^ |w-?74r % }+ c.c. (C2) 

Our goal is to calculate the integral over w and show that the result of integration indeed coincides with this correlation 
function, as given, for example, in Ref. l58l . 

The first step of the calculation, as standard in practical calculations of four point correlation functions, is to use 
global conformal invariance to assign specific values to three of the variables: 

771 = 00, 772 = 1, 773=0;, 774 = 0. (C3) 

x remains arbitrary and is effectively the only free parameter upon which the correlation function depends. Without 
loss of generality, we take a; to be a real variable (which can be analytically continued to the complex plane later if 
needed) satisfying 

< x < 1. (C4) 

When taking the limit r\\ — \ oo, we multiply the correlation function Eq. (|CJ2|) by (771771)^ to ensure a finite result, 
since the conformal dimension of the order and disorder operators is 1/16. This results in the following expression for 
the correlation function 

(/j(oo,oo)/j(1, l)a(x,x)a(0,0)) = (1 - x)% (1 - x)~% \x\* J d 2 w (1 - (1 - w)^ \x - w\~% \w\~? 

x)*(l- x)~*\x\* J d 2 w(l - w)~% (l - wy \x - w\~% \w\~% (C5) 

We now introduce the following convenient notation 

/1 = (1 — x)^xi(l — w)~i (x — w)~iw~i (C6) 

f 2 = (l - x)~ix^(l - w)*(x - 7«)~3 w ~4 } (C7) 

as well as the following matrix 

«-(;;)■ <c8) 
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Jx 1 




FIG. 3: The integration contours in Eq. (|C11|) . 



Then Eq. (jC5P can be rewritten in the following compact way using matrix notation 

(H(oo,oo)ii(l,i)a(x,x)tr(0,0)) = [ d 2 wfQf. 



(C9) 



We are now in a position to use the techniques developed by Mathur in Ref. |68| specifically to compute integrals of 
this sort. In that paper, it was shown that an expression in the form of Eq. (|C9| can be rewritten in the form of 
Eq. (|HT9]> . Hence, Eq. JCS]) becomes 



>(oo,oo)/i(l, l)a(x, x)a(0, 0) ; 



-IMQ(1 ~ Mx)- 1 ^ + 7 (1) Q(1 - Af2) _1 / (1) - I (2) Q (1 - M 2 ) _1 1^ + J< 2 >Q (1 - M 3 ) _1 1^ . (CIO) 
Here 1^ and 1^ are contour integrals defined according to 



I (1 



dw f a 



j(2) 



dwf a . 



(Cll) 



|«;|=1 



The contours of integration are shown in Fig. [3] Mi, Mi and M3 are the analytic continuation matrices of the 
functions /, defined in the following way: Take the function / just above the real axis and analytically continue it 
over w along a big circle centered at zero to the values just below the real axis. The new value fc is then given by 



fc = M 1 f, 0<w<x; fc = M 2 f, x < w < 1; f c = M 3 f, x > 1. 
In this case, these matrices are diagonal and can be calculated in a straightforward fashion, to give 



Mi 



i 0' 
,0 *, 



M 2 



-1 
-1, 



M, 



-i ' 
-i, 



(C12) 



(C13) 



All that remains is to compute the integrals of Eq. (|C11|) and substitute them into Eq. (|C10|) to find the answer. 
The appropriate integrals are computed by deforming the contours from \w\ — x to < w < x, and from \w\ — 1 to 
1 < w < 00. The resulting integrals are then standard and can be expressed in terms of hypergeometric functions— 



3 3/ 3 3 3 

(1-i)«i 5 i dw (1 ~ w)~^(x — w)~iw~i 

J \w\—x 

3 3/ 3 \ f X 3 3 

(1 — x)* x* ( — 1 + e~ % ~%~j J dw (1 - w)~ 4 (x - w)~*w~ 
Jo 

nzi'iiL,^ ,i -ir(i)r(i) ,^/3 i i 

V2e**{l-x)'x 8 r(l) (,4'4 ; 2 ;a: 

1 (.2/ 



(1 — x) s x' s ® dw (1 — w) 4 (x — io ) 4 w 4 



- ] / dw (1 — w) 4 (x — io ) 4 w 



j_ 3 / 

(1 - X) _ SXS ( -1 + 

_ pi 

y/2e iS *(l-x)-*x-* / dC C~ S (1 - C) _l (l - 
Jo 



V2e l * (1 — x) <•.!■ 



_ + r(i)r(i)^/ i i i 



e" 4 



r(i) 

* r (iTO ( i_ a . ) -* a .-* > /i + v T3 



r(|) 



3 _ 3 

4 w 4 



(l-x)sxs ® dw (1 — w) 4 (x — tu) 

>/|u)| = l 

(1 - x)^x^ ^ 

\/2e _< *(l -x)txi e - 4 T L f d((i(l-()-i(l-x() 

Jo 



1 — e y dw(l-w) 4 (x — w) 4 w 



V^e* 4 (1 - x) «i« 



r(!)r(i)^/3 5 3 



r(|) V4'4'2 



r ' A ' r> ' x 



e« (1 _ x) -i, 4 ^r 



(1 — x) 8 i ! i (1 — to) 4 (x — to) 4 w 4 

J\w\ = l 

(1 - x)~^x^l - e~ 1 ^ / dw (1 - w)i(x - w)~%w~ 



V2e- l2 i(l-x)-^xie~^ f d( C i (1 - ()i {I - x(Y 

Jo 

r ^ i 3r(i)r(|) /3 1 3 \ 
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We define the conformal blocks 



Jo ee e '¥ T{ ^ ] V2(l - x)-ix-ifl + VT^ = = v^, (C18) 

V 2 / 

r (l) r (|) A*,, N-i -i.A ,/i :_./or(2) 



This is the conventional definition of the Ising CFT's conformal blocks, up to the overall phases, which are unimportant 
here since they cancel with their conjugates, and the overall constant p 1 v2, which causes no problems since it 
enters as a common factor to both blocks. Employing vector notation, we can now write 



/W = y/2 K , I {2) = V2[ >1 . (C20) 




Substituting Eqs. ([UBj) . (1UT31) . (|C20)) into Eq. (fCTU|) . we finally find 

(Kvi,m)K r i2,fj2)v(v3,v3)(j(Ti4,m)) = FoFq -?vF\- (C21) 



This is indeed the right answer for this correlation function. Not only does it match the expected form in Eqs. (1139[) . 
(I140p , and (|141[) , it also coincides with the explicit expression for this correlation function worked out in Ref . l5o . 



Appendix D: Screening in a two-dimensional two-component plasma 



In this secti on, w e examine the screening properties of a 2D two-component plasma using field theoretical analysis 
(see, e.g. Ref. Ill8l ). The interaction energy of a number of electric charges q k located at positions z k in 2D is given 
by 



$ = - Yl ln 1 



(Dl) 



l<3 



Suppose we have a 2D plasma of a large but equal number of electric charges of magnitude q and —q at temperature 
T . The partition function of this plasma can be written as 



Z= Y[d-z, , 



a*,.*-*/ T 



IF 



T 52 Uli Mzi-Zj\ 



(D2) 



where qi are either q or —q. 

With the help of Eq. (|A3[) . we observe that this same partition function can be rewritten as 



,iqi<f>(zi,zi) 



(D3) 



while making the identification 



T = .g. 



(D4) 



This identification of T with g is not particularly necessary. We could have multiplied all the electric charges by an 
arbitrary factor, and simultaneously multiplied the temperature by the square of this factor. This would keep the 
partition function in Eq. (|D3[) exactly the same. However, T = g is a convenient choice that we use in this paper. 

The partition function of such a plasma is most easily computed in the grand canonical ensemble, using methods 
developed originally in the context of the Kosterlitz-Thouless transition^. For completeness, we present the derivation 
here. The partition function in the grand canonical ensemble looks like 



00 \ f n I 

Z = H E ^/IF 2 M ex p 

n=\ qi=±q ' fe=l \ 



■^2qi(/)(zi,zi) 



1=1 



(D5) 
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Here A is the fugacity, the parameter whose logarithm gives the chemical potential of the charges. We can now sum 
over qk = ±q and over n to findii^ 

J U(f> e~-k / d 2 z(V<P) 2 +2\} d 2 zcos(q4>) 

Z= JV^e-^I^^ 2 ' (D6) 

The behavior of this plasma depends crucially on whether A is a relevant or irrelevant perturbation. Since the 
correlation function 

( e lq ^ z '^ e -^K™) } = \z- w\~ q2/a , (D7) 

the dimension of this perturbation is 

AA = 2-f. (D8) 

The perturbation is relevant if A\ > or, equivalently, q < y/Ag. Otherwise, it is irrelevant. 

If the perturbation is irrelevant, then the plasma does not screen. If it is relevant, the plasma screens. In the latter 
case, the correlation function of two electric operators is a constant if they are farther away from each other than the 
correlation (screening) length £ 2 of the plasma 

250 e i 9i^( z i^i)g i '320(z2,22) e -£r/d 2 z(V0) 2 +2A/d 2 zcos(g0) ^ _j_ Q( e ~\^-^\/^\ 

where C 2 is a constant independent of Z\ and z 2 - To see this, we simply observe that if A is a relevant perturbation, 
typical values of <f> are restricted to the minima of cos (q(f>), and so </> fluctuates very little about this value. 

Indeed, upon identifying g with the temperature of the plasma, we see that it is at high temperature that the 
plasma screens. At low temperature, charges of opposite magnitude bind strongly pairwise to form neutral dipoles 
and hence the plasma does not screen. The transition temperature between these two phases is seen to be T C2 = q 2 /A. 

Additionally, one can consider the behavior of the magnetic operators in the screening phase. Their correlator goes 
to zero at large distances 

V^O m (z 1 ,z 1 )0- m (z 2 ,z 2 )e-^^ d2z ^^ 2+2 ^ d2zcos ^ ~ 0(e-l^-^l/^) . (D10) 

This corresponds to the well known fact that magnetic charges are confined in the electrically-screening phase. To 
see why this is so, recall that the operator O m {z, z) creates a vortex in the field 0, so that <j) winds by 2irm in going 
around (z, z). As was pointed out earlier, <fi is restricted to the minima of cos (q(f>) in the screening phase. Therefore, 
the two magnetic charges m and — m are necessarily connected by a region across which (j) changes by 27rm. Since (j> 
must leave the minima of cos (q<fi) to do this, it is energetically favorable for this region to be as small as possible. 
Consequently, this region will take the form of a string connecting the magnetic charges, with (f> rapidly changing by 
2-Km as it crosses the string perpendicularly. The energy cost of such a configuration is proportional to the length of 
the string and the magnetic charges are thus confined. In order that the integrand of the partition function for the 
plasma (of electric charges ±q) be single-valued when there are test particles with magnetic charges ±m included, 
the charges must satisfy a Dirac quantization condition: qm € Z. (This condition is satisfied for the plasmas with 
magnetic charges that arise in our plasma analogy mapping for the MR state.) This makes it clear that <fi can pass 
from one minima of cos (q(f>) to another as (f> changes by 27rm. If one does not satisfy the Dirac quantization condition, 
one must introduce a branch cut or deal with the multi-valued integrand in some other manner. It is not clear 
whether this can provide something physically meaningful, but, if so, the confinement argument still applies to these 
(fractionally quantized) magnetic charges in the plasma. 



Appendix E: Debye Screening Length of General Plasmas 

The considerations of the previous Appendix tell us whether or not a two-component plasma screens. When a 
plasma does screen (either one or two component, or even more general plasmas with multiple types of Coulomb 
interactions), we would like to know what its screening length is. Deep within the screening phase, we can use the 
Debye-Hiickel theoryi^S. We review this theory here, and generalize it to plasmas with multiple types of Coulomb 
interactions. We consider the case where there are m different types of Coulomb interactions and S different particle 
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species. The jth species particles have density rij (~r^) and carry the feth type of Coulomb charge <Tj ■ We start with 
the Poisson equation 



V 2 0« = --pW (^) 



(El) 



for the fcth electric potential (f>^ ( r ) and charge density ( r ). We take the convention in which the Coulomb 
energy in D spatial dimensions between two point charges q\ and q2 at r \ and fa, respectively, takes the form 



which corresponds to 



$ = < 



<7i<7 2 log |~r\ -^2 1 forD = 2 

for £> > 3 



(E2) 



eo = 



_i_ 

2tt 



4tt d / 2 



for L> = 2 



for D > 3 



(E3) 



Next, we assume that the system is in thermal equilibrium, so the particle densities are given by a Boltzmann 
distribution with respect to the Coulomb energies 



( r ) = ny exp 



•r- 



1 m 



(E4) 



where 7V*- is the homogeneous density of the jth species of particles (far away from test particles). 

Charge neutrality for the fcth type of Coulomb charge can be obtained either by balancing charge among the 



different species of particles to sum to zero (i.e. X^i=i 1 



0) or through a uniform neutralizing background 



charge density [i.e. p 



(fc) 

neutralizing 



(*)„(0)l 



In either case, charge neutrality allows us to write 



3=1 



(E5) 



giving zero charge density where the potential vanishes. 

Combining Eqs. (|E1[) . (|E4[) . and (|E5I) . we obtain the generalized Poisson-Boltzmann equation 



s ( m 

V 2 0W C?) = -J2lfnf 1 - exp -^E^ 

£ ° 3=1 I L i=l 



(E6) 



This differential equation for <^( fe ) is obviously non-linear, but is approximately linear when and where the plasma is 
weakly coupled (i.e. the Coulomb energies are small compared to the temperature) so that Y^iiLx , lf l <t>^ C?) *C T. 
In this regime, we can expand the exponential 



exp 



z=l 



(0, 



(E7) 



to obtain the linear approximation of Eq. f|E6l) 

m 

v 2 0«(^) ~ EW'H?) 



;=i 



A - 1 Y^J°U fe U') 
Afei - —fl^ n 3 1j 1] ' 

U 3=1 



(E8) 
(E9) 
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which generalizes the Debye-Hiickel equation. A is a symmetric, real, positive-definite matrix, so all of its eigenvalues 
A a (where a = 1, . . . , m) arc positive. It is now straightforward to solve this differential equation by changing to a 
basis in which A is diagonal, giving m independent diffusion equations 

V 2 ^ Q ) (?) ~ A a ^ Q > (?) , (EIO) 

where (f>^ — J2k Sak4>^ k \ f° r S the similarity transformation that diagonalizes A, i.e. SAS^ 1 = diag[Ai, . . . , A m ]. 
Requiring the potentials 4>^ to go to zero at infinity, we know the solutions must generally have an exponentially 

— 1/2 

decaying behavior, with decay lengths £ a = A a . We define the longest decay length to be the Debye screening 
length of the plasma 

£ D = m^x{\- 1 / 2 }. (Ell) 
We now consider several examples relevant to this paper: 



1. The One-Component Plasma 



For a one-component plasma, there is one species of particles with charge Q and a neutralizing background, which 
gives 

n^O 2 

For the v = 1/M Laughlin states, the corresponding plasma is a 2D one-component plasma with T = g and 
Q = ^/TM~g. This gives 

A = i (E13) 

l B 

where = v/2-k1 2 b is the electron density of the quantum Hall fluid. Thus, the Debye screening length for a 
Laughlin state is l\ = Is I V%- 



2. The Two-Component Plasma 



For a two-component plasma, there are two species of particles with charge Q and — Q, respectively, which gives 

A-*!?* (E14) 

For the p-wave superconductor, the corresponding plasma is a 2D two-component plasma with T — g and Q = a/%. 
This gives 

A = 127m(°) (E15) 

1/2 

where is the fermion density. Thus, the Debye screening length is £2 = [l27rn^°^l 



3. The Moore-Read Pfafflan States' Plasma 



For the v = 1/M MR Pfaffian states, the corresponding 2D plasma has temperature T = g, two types of Coulomb 
interactions, and two particle species: the first with charge = y/2Mg and = \/3g, the second with charge 
— Q( 2 ) = ~^/ig. There is also a neutralizing background for Coulomb charge of type 1. This gives 



A = 



27rn(°) 



(Q 



1 



2 a/67M 
a/67m 6/M 



(E16) 
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where — v/2n£ 2 B is the electron density of the quantum Hall fluid, which must also be the density of the screening 
charges. The two eigenvalues of this A are A± = (M + 3 ± \J M 2 + 9) /M£ B . Thus, the Debye screening length is 

f M \ 1/2 

Id = ; t B . (E17) 

U + 3- VM 2 + 9J V ' 

For M = 2 this is i D 1.2 t B . 



Appendix F: Conformal Blocks of n a fields and N ip fields: a basis of n quasihole wavefunctions in which 

braiding properties are manifest 

In this section, we compute the conformal blocks which correspond to n-quasihole (and TV-electron) wavefunc- 
tionsi^. The four-quasihole case was computed in Ref. [U here we extend this result to arbitrary even n . From the 
previous discussion, it is clear that we only need these conformal blocks in the two-, four-, and six-quasihole cases. 
However, for completeness, we compute them for arbitrary numbers of quasiholes. These wavefunctions have the nice 
property of furnishing, through thei r explicit analytic continuation, representations of the n-quasihole braid groups. 

The basic strateg y is to use Refs. Il2ll)l22l to compute 



(cricr 2 CT1CT2 • ■ • 0102 ipilpi ...ipt) 



(Fl) 



Here, we have two chiral Majorana fermions rpi, 4*2 with their two spin fields 01, 02. Since <j\ and 02 are completely 
independent, this is equal to the product of 



(01 Ol ...<Tt1pl1pi...1pi) 



(F2) 



and 



(02 ■■■02) 



(F3) 



The bosonization formulas derived in Refs. H2I123 allow us to compute all of the con formal b locks of Eqs. (|F1|) and 
(|F3p . thereby giving us the desired conformal blocks of Eq. (|F2|) . According to Refs. ll2lUl22l if we have two chiral 
Majorana fermions, ip±, ij>2, we can combine them into a single Dirac fermion which can be bosonized: 



(F4) 



so that ipi = cos ip. Bosonizing chiral spin fields is trickier. The individual spin fields do not have a simple expression, 
but the product of two factors of 0102 can be written in the form: 



01 



("i) 2 ("i) 01M02M =e^)l 2 e- l ^ 2)/2 ± e ~Mvi)/2 ^Am)/2 



(F5) 



In a conformal block in which a given set of fields 0102 ■ 0102 fuse to 1 • 1 we take the + sign; if they fuse to ip\ • "02, 
we take the — sign. 

Thus, we can compute the square of a conformal block of n fields by computing: 



'(PuPa.— .P n / 2 ) 



r.=0,l 



(-1) 



' P/ e i(-l) ri (Vi-V2)/2 e i(~l) r 2( l p 3 - ipi )/2 



c ,i(-lf"/=( v ,„_ 1 - y „)/2 



(F6) 



We have employed the shorthand cp^ = ^(n M ) and (— l) r p = (— l)^3 r 3 p i . The subscript (p 1 ,p 2 , ■ ■ ■ ,p n / 2 ) on the 
left-hand-side is used to specify the conformal block of this correlation function which we are computing: p. t = 0, 1 
denotes that the (2i — l) th and 2i th fields fuse to I or ip, respectively. There is an overall parity constraint 
Y!,iPi = 0(mod2) for Eq. (|F6I) . since there are no additional i\) field insertions in this correlation function. The two 
values r i = 0, 1 correspond to whether we have used the first or second term on the right-hand-side of Eq. (|F5I) in 
the bosonic correlation function on the right-hand-side of Eq. (|F6[) (which is equivalent to the usage of r i = 0, 1 in 
Section llXT) . Thus, we have 



(00 . . . 0, 



(pi,p 2 > 



I! %i-l,2j'-l ^24, 2j 
i<j 



n/2 

UV2 

i,3 



n/2 



e (- i ) r - p n 



J k,l 



(F7) 



k<l 
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where 



Xk,l 



V2k-l,2l 7 hl-l,2k 
^2fe— l,2i— 1 %fc,22 



(F8) 



This generalizes the formulas for the four-er conformal blocks in Ref. l58l which has only two terms in curly brackets. 
Now consider the correlation function 



(cri<T 2 0"lCT 2 • • • <J\<J2 ipilpi ...ipi) 



(F9) 



with TV Majorana fermion fields ipi = cos ip. We will initially consider the even N electron number case and briefly 
mention the odd N case at the end of this section. This can be computed in the same way as before by choosing 
half of the fermions to be e lip , the other half to be e~ lip , and then summing over all permutations, which gives the 
Pfaffian. Hence, we obtain: 



((T1CT2 <J\(T2 ■ ■ ■ CTlO-2 1pli>l ■■■Ipl) 



/n V 

11 "fei-l, 2j-l "iifij 
i<j 



CPl.Pa.-.P n/a ) 



n/2 

n %»-i. 

1,j 



2.1 



n/2 



J/ 2 ,T, 



E (-ir P lKr'"" *(l + n,3 + r 2) ...)(2- ri) 4-, 21 ...) 



r i= 0,l 



(F10) 



where 



(l+r 1 ,3+r 2 ,...)(2-r 1 ,4-r 2 ,...) 



Pf < 



Zi - Zj \rj 2 - Zi Tji - Zj 



% - Zi 7]4 - Zj 



+ 3) (Fll) 



Including the charge sector of the v = 1 /M MR Pfaffian wavefunctions and dividing by (aa . . . a) 
finally obtain: 



( Pl ,p 2 ,..., Pn/2 ) 



, we 



(n X 

11 72i-l,2j-l %i,2j" 



( Pl ,p 2 ,..., Pi ) 



n/2 

n 



1/2 



-1/2 



E (-ir p Il 



J k,l 



J 



k<l 



■i/2 



(1+^,3+^,.. 0(2-^,4-^,. ..) 



(F12) 



r i= 0,l 



fc<Z 



where 



n N 



N 



,1/2 



*(l+r 1 ,3+r a ,...)(2-r 1 ,4-r 3 ,...) = ^(1+^,3+^, ...)(2-r 1 ,4-r 2 ,...) x ]J II ^ 

/i— 1 2—1 Z<J 

^ ('7l + r 1 - Z 0('73 + r 2 - Z 0---(';„-l + r„ ~ Z (»?2 - r j - Z j ) (Vi - r 2 ~ Z 3 ) • • • (.Vn - r „ ~ 2 3 ) + (^j) ^ N 



n, ,M 3 £ ;l z »l 2 



Pf 



n( 



:,-z,) M e ^^ |Zl|2 (F13) 



are electron wavefunctions with normalizations that do not explicitly contain the quasiholes' braiding statistics. 
This equation expresses the 2^ _1 basis vectors with TV even in terms of the 2? _1 basis vectors 

\Pl 1P2 ''"'Pn/2' 

^ (l+r ,3+r 2 ,...)(2-r ,4-r ,...)■ The basis vectors *&(x+r 3+r ...)(2-r ,4-r ,...) are intuitive and easy to write down (and 
were, therefore, written down in Ref. 1431) . However, they are not orthonormal and their braiding properties are com- 
plicated. Meanwhile, according to the result shown in this paper that Berry's matrices are trivial in the conformal 
block basis, the basis vectors "J^ p p j have simple braiding properties given by the branch cuts in their definition 

Eq. (lFT2j) . 
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In order to compute wavefunctions for an odd number TV of electrons, we can compute the wavefunction for N — 1 
electrons and n + 2 quasiholes as before, and then obtain the desired wavefunction by taking rj n+ i — > r/ Il+2 = z±, 
dividing by the appropriate power of r} n +i — Vn+2, and correcting the charge sector terms so that z\ corresponds to 
an electron coordinate: 

n N 

*(p 2 --.p n/2 )^ 1 ' ■■■,Vn,z 1 ,...,z N )=l[ ( Vfl - Zl )^ J] ( Zl - z.f- 1 e -^l^l 2 

|U=1 i=2 

x lime"(™-+l) ; 1} (?7i, . . .,r]n,Vn+i = zi +e,i] n+2 = z 1 ;z 2 ,.- .,z N ) (F14) 

The piS which index the odd N electron number wavefunction on the left-hand-side satisfy the parity constraint 
ErirP, = l(mod2). 



Appendix G: Incomplete Direct Approach to Plasma Analogy for n Quasiparticles 

In this section, we present a more direct approach to constructing the plasma analogy for n quasiparticle wave- 
functions. Unfortunately, the argument is incomplete, but presenting the argument serves to clarify the obstacle in 
proceeding in this manner, and why we needed to use the methods involving disorder operators. 

We consider the conformal blocks with n a operators (which correspond to n fundamental quasiholes). There are 
2'?~ 1 such conformal blocks, which we denote by: 

J^aiVn'i z i) = (c(J7i)o-(»72) • • • a{q n -i)a{r) n )il){zi) . . . ip(z N )) a (Gl) 

where a = (771, . . . , 7r n / 2 ), for Wj = 0, 1 with the overall constraint 7r„/ 2 = for N even and 7r„/ 2 = 1 for TV odd. We 
can represent these in the following way (defining 7r =0): 

/ n/2 y 

^ 1 ,...^ /2) (^; ^) = ^ J] tfr'-^mj-jv^fa) v 3 f( Zl )v™(z2) ■ . ■ v2?{z N ^)v™{z N ) J (G2) 

n/2 

= II r f 1 ^ U2 r 1 ■ f duN-i f dun 

j=i Jc Tj Jc zi Jc zi Jc Z3 Jc Z3 Jc ZN1 Jc^N-i 

I n/2 

( II V2 ?(mj-i)e iil ~* j - l)a - Hwj) V2™(mj)e iWja - Hwi+l) 

x V^{z 1 )e ia -^ ) e ia -^ ) V^{z2) . . . V^{z N - 1 )e ia -^ UN - l) e ia -^ UN ] V^{z N )} (G3) 

n/2 

= TT f dwj f dui f du 2 f du 3 f du 4 . . . f du N -i f du N f^ 1 ^..^ n ){w a ;u i ;n ll \Zi)(GA) 
i=\ J c Tj Jc zl Jc zi Jc Z3 Jcz 3 Jc ZN _ t Jc.^ 

where C x is the contour at radius \x\ centered on the origin, rj = r]2j-i- 7rj _ 1 , and we have defined 

/n/2 

\j=i 

x V^{z 1 )e ia -^e ia -^V^{z 2 ) . . . l / 3 TKi)e ,a -^- l) e" 1 ^ ( "" ^T^iv)) • (G5) 

There are iV screening operators with coordinate Ui for the ip fields and n/2 screening operators with coordinates w a 
for the a fields, and the conformal blocks are determined by the placement of the w a screening charge contours. 
Specifically, ttj = 0, 1 indicates that the contour for Wj is at radius | ??2j— 1— 7Tj_i |- Strictly speaking, f a does not encode 
the fusion channel without the knowledge of these contours, but we will nonetheless use the subscript label to remind 
us of the contour placements. 



G2 



We want to compute the overlap 

N 



G a ,p (v^v'v) = / Yl d2 z k J 7 a (ri^z l )J 7 p('q f _ l ;z. l ) 

J k=l 

IT d 2 z k IT d> duij f dw'j f du[ (f> du' 2 j> du' 3 (p 

k=l j=l Jc ri J Cr> } JC Z1 JC Z1 JC Z3 JCz 



du' 4 . . . d> du N _ 1 <b du' N 



x® dux f du 2 f du 3 i> du 4 . . . i> du N -x f du N f a (w a -,u i ;T] fl ,;Zi)f0(w' a ;u' i ;r] pi \Zi) (G6) 
'CW Jc zi Jc Z3 Jc Z3 Jcz N _ 1 Jc^^ 

71,2 1 I 

FT* dw 3 f dw' j T a ,p{w a ,w' a ]r] IJ ,,r)tj) (G7) 
f = \Jc rj Jc r , 



where r'j = °r\2j-\-ix'._ 1 corresponds to j3 — (ir[, . . . , 7r^ 2 ), and 



TV 

r a ^(w a ,w' a ]f]^,r)^)= T\d 2 z k (f> du[ f du' 2 f du 3 f du' A ...f du' N _ 1 f 
J k=1 Jc zl Jc zi Jc Z3 Jc» 3 Jcz N _ 1 Jc n 



du' N 



x j> dui f du 2 f du 3 d> du4...d> du N -i f du N f a (w a ; Uf, i]^; z { ) fp (w' a ; u-; r^; Zi) . (G8) 
J c*i Jc zi Jc Z3 Jc Z3 JCz N _ 1 Jcz N _ 1 

Now let us see if we can convert the pairs of u' i: Ui contour integrals into d 2 Ui integrals. If we can do so, then we 
can deduce the properties of T at p and G a ^ via plasma analogy. We follow Mathur's steps, starting from the integral 

d 2 ui . .. J d 2 u N f a (wjiUiifj^, fp {w'^u^-q^ z 4 ) . (G9) 

In this case, we are not considering screening operators of non-Abelian fields, and thus monodromies will not take one 
to a different conformal block, so we do not need the bilinear form notation. Also, we are not restricting our attention 
to diagonal components, so the plasma potential will not be monodromy invariant (i.e. not single- valued) , but this is 
not necessarily a problem as long as we keep track of monodromies. We will show that this expression decomposes 
into 



d 2 ui . . . J d 2 UNf a (w a ] Ui',Vn> z i) fp ( w a'> u i')V^ z i) = (terms with contours all at radii \ zi\) 

+(terms with at least one contour at radius \f]^\) + J-terms at branch cuts. (G10) 

The terms with the contours all at the radii \zi\ are equal to (u) a , w' a ; fj^, r)^) (up to some constant). This is 
because the contributions to this term from integrals with contours on the wrong zi will vanish, leaving only r a; ^. 
This is known from Ref. l67l . since contours in the unallowed configurations give rise to overall multiplicative terms 
with canceling phases. The terms with at least one contour at r/^ will vanish after taking the §c §c i dw j 

contour integrations, because this results in unallowed configurations of screening contours (i.e. too many screening 
charges on the V21 operators). Most of the J-terms will cancel each other, just as in Mathur's construction, except 
for the ones on either side of a branch cut. The branch cuts occur when w a ^ w' a . (We note that w\ and w[ can 
actually be treated the same as the Ui screening charges, since its contour placement is uniquely specified). We want 
to show that, in the end, the J-terms vanish or cancel, at least after performing the dw' a and dw a integrations. 
If we could make the J-terms vanish, then Eq. (|G7|) would become 



G a ,/j (^>i) =11 f dw J ? dw 'j^o ( ,p{w a ,w' a \'q IJ ,,r] IJL ) (Gil) 



where one can now apply the plasma analogy to 

N „ N 



f a ,p{w a ,w' a - 1 -q ll ,-q ll )= j Y\ d2z kd 2 u k f a (w a ]Ui;r]^ Zi) fp (w'^u^r)^; Zi) = / Y\_d 2 z k d 2 u k e * /T = e F . (G12) 

** I 1 •* 1 1 



fc=l fc=l 
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Here, T = g and $ describes the 2D Coulomb interaction between N charge Q = y/3g particles at Zi, N charge — Q 
particles at m, n charge Q/2 particles at ?7 M , ^ particles with electric charge —Q/2 and magnetic charge m = yJi/Ag 
at w a , and ^ particles with electric charge —Q/2 and magnetic charge —m at w' a . Hence, F is the free energy of a 
classical 2D two-component plasma at temperature T of N charge Q particles and N charge —Q particles, with n 
charge Q/2 test particles at 77 M , ^ test particles with electric charge —Q/2 and magnetic charge m at w„, and ^ test 
particles with electric charge —Q/2 and magnetic charge — m at w' a . By confinement of magnetic charge in a screening 
plasma, we know that Eq. (|G12[) will vanish unless w a — w' a , which shows that the result is proportional to S a p, as 
desired, since a ^ (3 requires w a ^ w' a for at least one a. 

We now proceed by seeing what happens for a single u screening charge. We partition the plane into a number of 
annular regions Di such that the positions of the other coordinates where there are potentially singularities or branch 
cuts are left outside of Di. Leaving the other coordinates implicit, we have for a region D 

d 2 uf a (u) fp («) = J / duf a (u) fp (u) (G13) 

where 



d 2 jq d 



fp (u) = / du'fp («') . (G14) 



p 



We define 







J a/3 — 


/; 


T c - 


/; 



duf a (u) fp (u) (G15) 

dufp (u) (G16) 
for C a contour running from P to P'. Now taking the same steps as Mathur, we get 

dVa (u) i> («) = | [jg + + J* - J^J| 2 - - Jf|] (G17) 

(G18) 



D 

i 
2 



tCi tC 2 , tSi tS 2 iu r lrCird n , f r l FC 2 7-C2 1 M " ~ M /3 TC 2 rCi 

•'a/fl _ J Q0 + -Vi - ^ ~ I 1 - M a) i a 1 p ~ U - JW /3j J q J /3 + (1 _ j^Q (1 - Mp) 13 



where we used Jg 1 + /g — 7g 2 — I^ 2 =0 and Ig 1 = Mpl^ 2 , and now M a and are not equal for the regions D; 
between to w'j +1 (i.e. from radius 1 772 j | — e to 1 772^+1 1 — e) when 7Tj 7^ 7^, but are otherwise equal. This is not a 
problem, because (as previously mentioned) this only gives us an extra term with a contour at r]2j and one at 772.7+1 , 
but these give a vanishing result when one evaluates the dwj+\ and dw'j +1 contour integrals. 

The J-terms are however a more difficult problem. We know that these will cancel as long as Wj = w'j , since then 

J^p = J^p and J^p from region Di is equal to J^ 2 . from region Di+\. However, when wj 7^ w'p there is a branch cut 
running between Wj+i and w'j + i and the J-term on the two sides of the cut do not cancel each other. Thus, we are 
stuck with a J-term integrated around these branch cuts, and no obvious way to cancel them out. 

Appendix H: Explicit Examples of Orthogonality for Unmatched Quasiparticles 

In this appendix, we provide the derivation of the overlap Eq. (|211l) for wavefunctions describing the v — 1 filled 
Landau level with n holes and of the overlap Eq. (|212[) for wavefunctions describing an arbitrary quantum Hall state 
(that has a plasma analogy) with one Laughlin-type quasihole. 

1. v = 1 Integer Quantum Hall State 

We first consider the v = 1 filled Landau level state, which can be solved exactly. The wavefunction for one filled 
Landau level of TV electrons with n holes is 

*i (*/>; *o = n fa* - ^) fi n - *o n ^ - *>) e 1 * ^ ]zA ■ ^ 

p.<v n=l i=l i<j 



G4 



Taking the inner product of two such wavefunctions with holes not necessarily at the same positions, one finds 

N 



J k=l 

/N n n N N 1 A /, ,2,1 / | 2N \ 1 y> 1 12 

if** n - ^ « - ^)] n n ^ - *) « - *)] n i* - ^ 2 e - A M j ~ 5 £ 

fc— 1 /^<^ /x=l i— 1 «<j 

n 

= d V (-i) w TT e -Ki^(rti 2 +i<i 2 -2^ (f .)<) 

7reS„ M=l 

n 

= C 1 (2tt)" J] (-^"Il <*lll (^) - <) . (H2) 

where Ci is the (unspecified) normalization constant. For this, we note that one can think of wavefunction with N 
electrons and n holes as one filled Landau level of N + n particles, which is a Slater determinant state. Thus, the 
integral over the N electron coordinates z\ gives the 2n-point particle correlation function, which factorizes into a 
product of 2-point functions summed over permutations, as per Wick's theorem. Finally, we used the fact that the 
(normalized) two-point functior>i2£ 

is the lowest Landau level projection of the delta-function, in the sense that 

d 2 z St Lh (z - z') f (z) e-^l 2 = / (/) e-^'l 2 (H4) 
for any holomorphic function f(z). We also note that 

e -i(\z\ 2 + W\ 2 ~2zz') = e - i \z-z'\\ i (zz'-zz') = e -±\z-z>\ 2 e >±l m [zz'} ^ ^ 

so the two-point function has Gaussian decay with the distance between z and z' . Hence, we find that 

||*! (77 M ;^)I| 2 = Ci + O (e-l""-"-l 2 /<) . (H6) 

This is a somewhat stronger result than given by the plasma analogy, which nominally involves 0(e _ ' I?A '~' ? ' / '^ 1 ) 
corrections. 

One could also arrive at the result of Eq. (|H2[) by noting that 

N n N N _i « 2 

r = / n^** n n - *) k - «)] n i* - ^ e 2 1 m 

is holomorphic in 77^, antiholomorphic in 77^, and 

r (*7m>^) = r = r (^O)'^) = r ( H8 ) 

for any 7r g S n . By the plasma analogy, we know that G (fj^ , 77^ ) = Ci + O (e I* 7 ** ' ' ) and thus can uniquely 
obtain the result of Eq. (|H2[) . i.e. that 



n — — ( 1 i 2 +l ' I 1 

n n 
M< 1/ 7rGS„ /j=1 

where C\ can now be identified as the undetermined constant from the plasma analogy, by using only the plasma 
analogy and the analytic properties of T. 
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2. Laughlin-type Quasihole 



For the v = 1 /M Laughlin states, Laughlin demonstrated such an orthogonality for the single quasihole wavefunction 
using the plasma analogy and analyticity££. Specifically, he showed that 



N 

J fe=l 



N 



JY 



k=l i=l 
(7 ie -I^(l')| 2 + l'7T-277 I 7') i 



I |2M 4M 



(H10) 



where Ci 



1 VI 2 ), the 



/ 1 2 

^ i (??; Zj) • For this, he noted that, except for the Gaussian factors exp {—jjj \r}\ - 

inner product is holomorphic in rj and antiholomorphic in 77, and the plasma analogy indicates that G(fj, rj) = C\ 
These properties uniquely determine the result of the inner product. One similarly has that 



_J_ e -^(\v\ 2 +\n'\ 2 -mv') = 5 2 ( _ ,s 
2nM LLLm w ; 



(Hll) 



is a projection of the delta- function into a lowest Landau level with a re-scaled magnetic length of = V~M£b, be. 



d 2 V 5t LhM (v - rf) f fo) e-^ 2 = f (tf) e ™IVP 



(H12) 



for any holomorphic function f(r)). 

The same argument applies to wavefunctions with one Laughlin-type quasihole for any state with a plasma analogy 
(e.g. the I2 excitation in the MR state). Specifically, one has 



G(fj,v')= f 11^% fa*)** 
•* fe=i 

r N N j 1 1 2 \ 

j 1 — 1 i — i 



fe=l i=l 



*i (^) 



(H13) 



where (z^) is the ground-state wavefunction and (?y; z,) the wavefunction with one Laughlin-type quasihole 
at r/. 

It is difficult to generalize these methods of obtaining explicit overlap results that go beyond the qualitative behavior 
obtained in Section [X] for cases that involve multiple Laughlin-type quasiparticles or different types of quasiparticles. 
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The self-energy for the neutralizing background charge density are not included in the Coulomb potential energy <3?i . This 
can be done safely without altering any of the subsequent arguments, since it simply contribute a constant to this energy. 
More generally, in the standard basis, one should label every fusion channel, i.e. also have a label for the collective topological 
charge of the 1 st through (2j + l) th fields. However, for the Ising a fields, these fusion label are always equal to a as a 
consequence of the simple fusion algebra, so we do not bother to label them. 

For example, tensor product (fusion) of two adjoint representations of SU(3) give 8<8>8 = 27©1O0 1O08©8©1, which 
therefore has iVf 8 = 2. 

In fact, these are the only modular non-Abelian anyon models with 3 topological charge types. 

We note that when the topological S- matrix, S a b = (5H C dc) ^ J2 C Nabdcj^, is unitary, the theory is called "modular" 
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and obeys the relation N£ b = J2 X S ax Sb x S* x /Sox- 

As we were completing this paper, we learned that the conformal blocks and associated electron wavefunctions computed 
in this Appendix were obtained recently by a slightly different method by E. Ardonne and G. Sierra—^. 



